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Ductile  phase  toughening  of  MoSi2  matrix  reinforced  with  Nb  filaments  or 
laminae  for  high  temperature  structural  applications  has  been  investigated.  To  improve 
the  high  temperature  stability  of  the  composites,  AI2O3  coating  between  the  matrix  and 
reinforcements  was  identified  as  an  effective  diffusion  barrier  and  thermomechanically 
compatible  with  both  the  matrix  and  the  reinforcements.  The  techniques  investigated 
for  producing  the  AI2O3  coating  on  the  surfaces  of  niobium  reinforcements  included 
sol-gel  processing,  physical  vapor  deposition,  and  hot  dipping  in  a molten  aluminum 
followed  by  anodizing  to  form  alumina.  The  optimization  of  these  coating  techniques 
was  conducted  and  the  results  indicated  that  the  best  coating  was  provided  with  physical 
vapor  deposition,  followed  by  the  sol-gel  technique. 

In  order  to  optimize  the  ductile  phase  toughening,  the  interactions  between 
cracks  and  ductile  reinforcements  were  examined.  The  effects  of  various  material 
parameters,  such  as  size  of  the  ductile  phase,  bonding  strength  of  the 
matrix/reinforcement  interface,  and  orientation  of  the  ductile  laminae  on  the  toughening 
were  also  explored.  The  experimental  results  indicated  that  the  ductile  phase  toughening 


in  the  present  system  was  large-scale  bridging  in  nature;  the  crack  growth  resistance 
increased  with  size  of  ductile  phases  for  large  initial  cracks;  ductile-laminae-reinforced 
composites  provided  a two-dimensional  toughening;  and  the  matrix/reinforcement 
interface  exhibited  a dual-effect  on  the  toughness,  that  is,  a strong  interface  would  give 
rise  to  a high  damage  tolerance  and  low  work  of  fracture,  while  a weak  interface  would 
provide  a high  work  of  fracture  and  low  damage  tolerance. 

Based  on  the  experimental  observation,  a rigorous  analytical  model  was 
proposed  to  describe  the  large- scale-bridging  (LSB)  of  the  ductile  phase  toughening. 
The  LSB  model  allowed  for  the  prediction  of  the  typical  R-curve  behavior  of  the 
composites,  encompassing  the  effects  of  the  bonding  strength  of  the  matrix/ 
reinforcement  interface,  the  size  of  ductile  phases,  and  the  intrinsic  mechanical 
properties  of  the  reinforcement  and  matrix.  The  LSB  model,  which  was  found  to  be 
consistent  with  the  experimental  results,  provides  a guideline  for  the  design  of  ductile- 
phase-toughened  composites. 
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CHAPTER  1 
INTRODUCTION 


The  need  to  increase  the  operating  temperature  of  heat  engines  to  improve  their 
efficiency  has  provided  a large  impetus  in  recent  years  for  the  development  of  high 
temperature  structural  materials  such  as  carbides,  nitrides,  aluminides  and  silicides. 
These  ceramics  and  intermetallics  have  very  high  melting  temperatures,  high  elastic 
moduli,  low  densities,  and  can  withstand  hostile  environments.  Their  structural 
application,  however,  is  currently  limited  by  their  tendency  to  fail  catastrophically  from 
small  surface  and/or  internal  flaws,  and  by  thermal  shock. 

Accordingly,  much  effort  has  been  made  to  improve  the  toughness  of  these 
materials.  Two  general  approaches  have  received  a great  deal  of  attention:  (1) 
increasing  toughness  through  composite  approaches  (e.g.,  introduction  of  a second 
phase  in  the  form  of  particles,  flakes,  whiskers,  fibers,  or  even  voids  to  produce  a 
ceramic  matrix  composite  or  intermetallic  matrix  composite);  (2)  reducing  flaw  size 
and  population  through  improved  processing  (e.g.,  elimination  of  agglomerates, 
pressure  sintering,  etc.).  Although  significant  improvements  in  minimizing  strength- 
degrading defects  have  been  achieved  through  improved  processing  in  recent  years,  this 
approach  is  reaching  its  limitation  since  the  critical-sized  flaws  cannot  be  eliminated 
entirely.  In  contrast,  composite  approaches  have  proved  to  be  effective  in  improving 
the  toughness  and  have  increased  the  apparent  fracture  toughness  of  brittle  materials  to 
useful  levels  (about  20  MPa.m1/2).  However,  further  advances  are  needed  to  develop 
ceramics  and  intermetallics  with  the  desired  level  of  reliability  for  structural 
applications. 
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It  has  been  established  that  the  mechanisms  responsible  for  composite 
toughening  consist  of  crack  deflection,  crack  bowing,  crack  bridging,  crack  blunting 
and  crack  shielding  (i.e.,  phase  transformation  toughening  and  microcracking 
toughening).  Among  these  mechanisms,  crack  bridging  has  been  shown  to  be  the  most 
effective  toughening  mechanism  in  the  sense  that  it  could  increase  the  toughness  of 
composites  to  a useful  level  and  provide  such  toughening  at  ambient  temperatures  as 
well  as  at  elevated  temperatures.  Crack  bridging  is  achieved  generally  by  the 
incorporation  of  ceramic  whiskers  or  fibers,  or  ductile  fibers  or  particles.  Much 
exploration  has  been  made  in  the  area  of  ceramic  fiber-  or  whisker-reinforced 
composites.  For  example,  SiC  whiskers  have  been  used  to  reinforce  alumina  [1], 
mullite[2],  silicon  nitride  [3-6],  glass  [7,8],  magnesia-alumina  spinel  [9],  zirconia 
[10,11],  and  mullite-zirconia  [12,13].  Continuous  SiC  fibers  have  also  been  used  to 
reinforce  silicon  carbide  [14-17],  silicon  nitride  [18]  and  glass  [19],  One  of  the  most 
important  findings  in  the  brittle-fiber-  or  whisker-reinforced  brittle  matrix  composites  is 
that  a brittle  matrix  bonded  to  a brittle  reinforcement  by  a brittle  interface  results  in  a 
tough  material,  a result  defying  intuition.  Another  finding  is  that  the  toughness  of  fiber- 
or  whisker-reinforced  composites  is  flaw-size  independent.  In  contrast,  investigation  on 
ductile  phase  toughening  has  received  relatively  less  attention  than  ceramic  fiber  or 
whisker  toughening.  To  obtain  an  optimal  ductile  phase  toughening,  it  is  necessary  to 
consider  the  following  effects  on  the  toughness:  (1)  matrix/reinforcement  interface,  (2) 
morphology  of  ductile  phases,  (3)  size  of  ductile  phases,  and  (4)  intrinsic  mechanical 
properties  of  ductile  phases.  Besides,  most  of  the  theories  of  ductile  phase  toughening 
[20-26]  are  based  on  fracture  mechanics  of  small-scale  bridging  which  is  usually  not 
representative  of  ductile  phase  toughening.  Therefore,  more  accurate  models  (i.e., 
models  of  large-scale  bridging)  are  needed  to  predict  the  overall  mechanical 
performance  of  the  composite  from  the  properties  of  the  constituents  (matrix, 
reinforcement,  interface).  Finally,  there  are  often  deleterious  interactions  between 
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ductile  reinforcements  and  intermetallic  matrices  during  the  processing  and/or  at 
elevated  service  temperature.  Thus,  it  is  necessary  to  establish  a coating  technique  to 
prevent  the  interactions  and  maintain  the  integrity  of  the  reinforcement. 

It  is  the  goal  of  this  research  to  explore  ductile  phase  toughening  in  brittle 
matrices  for  structural  applications  above  1200^C.  Owing  to  the  projected  use 
temperatures,  the  composites  must  exhibit  adequate  mechanical  properties— high 
specific  strength,  toughness,  good  thermal  shock  resistance,  creep  resistance,  fatigue 
resistance-as  well  as  stability  in  oxidizing  atmospheres.  Thus,  molybdenum  disilicide 
(MoSi2)  has  been  chosen  as  the  matrix  material  in  the  present  study,  because  of  its 
excellent  oxidation  resistance,  high  melting  point  (2020^0)  and  relatively  low  density 
(6.31  g/cm3).  Similarly,  the  suitable  ductile  reinforcements  are  limited  to  the  refractory 
metals,  notably  W,  Ta,  Mo  and  Nb.  Nb  is  of  particular  interest  because  its  thermal 
expansion  coefficient  is  closest  to  that  of  MoSi2,  reducing  the  potential  for  matrix 
cracking  upon  thermal  cycling  during  processing  and  in  service.  The  work  presented  in 
this  dissertation  consists  of  four  major  elements: 

(1)  establishment  of  a processing  technique  and  the  process  control  parameters 
for  producing  a diffusion  barrier  between  the  MoSi2  matrix  and  Nb  reinforcement; 

(2)  study  of  the  flow  behavior  of  the  constrained  Nb  and  the  effects  on  the 
toughening; 

(3)  investigation  of  the  effects  of  various  material  parameters,  e.g.,  size  of  the 
ductile  phase,  bonding  strength  of  the  matrix/reinforcement  interface,  and  orientation 
of  the  ductile  laminae  on  the  toughening; 

(4)  development  of  models  of  the  bridging  micromechanics  and  guidelines  for 
the  design  of  ductile-phase-reinforced  composites. 

For  controlling  the  matrix/reinforcement  interaction,  chemical  compatibilities  of 
various  potential  coating  materials  (mullite,  AI2O3,  and  Z1O2)  with  MoSi2  and  Nb 
were  conducted  via  microstructural  and  compositional  analyses.  Based  on  these 
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analyses,  AI2O3  was  chosen  as  the  final  diffusion  barrier  coating  because  it  was 
chemically  compatible  with  both  MoSi2  and  Nb  plus  it  also  had  a thermal  expansion 
coefficient  close  to  those  of  both  the  matrix  and  reinforcement.  The  techniques  used  for 
producing  the  AI2O3  coating  on  the  Nb  surfaces  included  (1)  sol-gel  processing;  (2) 
physical  vapor  deposition;  (3)  hot  dipping  the  reinforcement  in  molten  aluminum, 
followed  by  anodizing  to  form  AI2O3.  The  process  control  parameters  for  the 
techniques  and  the  effectiveness  of  the  coatings  as  a diffusion  barrier  were  evaluated. 

Study  of  the  flow  behavior  of  the  constrained  Nb  was  conducted  using  a tensile 
test  on  a single  Nb  lamina  imbedded  in  an  MoSi2  matrix.  The  effects  of  the 
matrix/reinforcement  interface,  the  mechanical  properties  (yield  strength  and  work 
hardening  rate)  and  size  of  the  ductile  phase  on  the  flow  behavior  were  evaluated. 
Based  on  the  observations  from  the  tensile  tests,  an  approximate  model  was  proposed 
which  gave  insight  into  the  influence  of  yield  strength,  work  hardening, 
matrix/reinforcement  interfacial  bonding  strength  and  size  of  the  ductile  reinforcement 
on  the  flow  behavior.  The  stress-displacement  function  developed  was  used  to  calculate 
the  bridging  stress  across  the  crack  surface  in  the  large-scale-bridging  model  proposed 
in  this  study. 

The  role  of  the  matrix/reinforcement  interface  in  the  fracture  toughness  of 
MoSi2/Nb  composites  was  assessed  via  four-point  bend  tests  on  chevron  notched 
specimens.  Variation  of  the  interfacial  bonding  was  obtained  by  depositing  different 
oxide  coatings  (AI2O3  and  ZrC>2)  or  by  the  development  of  a reaction  product  layer 
between  the  matrix  and  reinforcement.  Measurement  of  fracture  energy  (bond 
toughness)  of  the  interfaces  was  carried  out  on  a chevron  notched  short  bar  specimen.  It 
was  established  that  whether  or  not  a strong  interfacial  bonding  was  conducive  to 
toughness  depended  on  the  criterion  used  to  describe  the  toughness  of  the  composites. 
Four-point  bend  tests  were  also  used  to  evaluate  the  effects  of  the  size  and  orientation 
of  ductile  laminae  on  the  toughness  of  the  composites.  The  results  indicated  that  the 
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toughness  of  the  composites  increased  with  increasing  size  of  the  Nb  laminae  and 
ductile  laminae  offered  two  dimensional  toughening. 

All  the  toughness  evaluations  and  examinations  of  the  interaction  of  cracks  with 
ductile  reinforcements  indicated  that  (1)  ductile  phase  toughening  involved  large-scale 
bridging,  i.e.,  the  bridging-length  was  at  the  same  magnitude  as  the  crack  length,  the 
specimen  size  or  the  distance  from  the  crack  to  the  specimen  boundaries;  and  (2)  the 
toughening  of  the  brittle  MoSi2  matrix  with  Nb  could  be  described  by  the  equilibrium 
stress  distribution  across  the  crack  surface.  According  to  these  observations,  a general 
approach  was  proposed  to  compute  the  ductile  phase  toughening  in  the  case  of  large- 
scale  bridging.  The  approach,  which  allowed  for  the  prediction  of  the  typical  R-curve 
behaviors  of  ductile-phase- toughened  composites,  encompassed  the  effects  of  the 
bonding  strength  of  the  matrix/reinforcement  interface,  the  size  of  ductile  phases  and 
the  intrinsic  mechanical  properties  of  the  ductile  phases  and  matrix. 


CHAPTER  2 
LITERATURE  SURVEY 

2.1  Theoretical  Background 

The  basis  for  toughening  brittle  materials  can  be  best  presented  using  a fracture 
mechanics  approach  for  crack  propagation  within  the  materials.  Using  this  approach, 
the  stress  state  near  a sharp  crack  is  related  to  the  size  of  the  crack  and  the  applied 
load,  and  is  described  in  terms  of  a stress  intensity  factor.  The  relation  between  these 
three  parameters  is  given  in  the  following  form  [27] 

K = payfjra  (2.1) 

where  K is  the  stress  intensity  factor  (MPa.m1/2),  /3  = (3  is  a dimensionless 

number  depending  on  the  geometry  of  the  cracked  part,  a is  the  characteristic 
dimension  of  a crack  (m),  L is  a unified  length  dimension  describing  the  geometry  of 
the  cracked  part  (m),  and  a is  the  applied  load  (MPa). 

A cracked  part  will  fracture  when  the  remotely  applied  load  increases  to  a 
critical  level  (Of)  at  which  the  stress  intensity  factor,  K,  reaches  a critical  value  (KJ. 
The  critical  stress  intensity  factor,  Kc,  is  defined  from  eq.  (2.1)  as  the  value  of  K at 
which  a failure  stress,  Of,  causes  a crack  of  size  2a  to  propagate,  leading  to  failure  of 
the  cracked  part.  According  to  this  definition,  a material  fractures  when  the  following 
condition  is  satisfied  [27] 

K > Kc  (2.2) 

1^  is  a material  property  and  is  commonly  referred  to  as  the  material's  fracture 
toughness.  The  fracture  toughness  of  a material  is  directly  related  to  the  fracture 
surface  energy.  Based  on  the  conservation  of  energy,  an  equation  describing  the 
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relation  between  failure  stress,  Of,  and  the  fracture  surface  energy,  y,  for  brittle 
materials  can  be  derived  from  Griffith's  theory  which  takes  the  form  of  [28,29] 


where  E and  v are  Young's  modulus  and  Poisson's  ratio  of  the  material,  respectively. 
Comparing  eqs.  (2.1)  and  (2.3),  it  is  evident  that  fracture  toughness  increases  with 
increasing  fracture  surface  energy.  Fracture  toughness  is  also  related  to  the  plasticity  of 
the  material.  Plastic  deformation  at  the  crack  tip  blunts  the  crack,  cutting  off  the  elastic 
singular  peak  of  the  stress  distribution  near  the  crack  tip  [30,31].  Therefore,  much 
higher  external  stress  has  to  be  applied  to  fracture  the  material. 

The  J-integral  proposed  by  Rice  [32,33]  provides  a convenient  way  to  handle 
fracture  of  materials  with  plasticity.  The  J-integral  is  an  energy  line  integral  defined  for 
two-dimensional  problems  (Fig.  1)  and  is  given  by  [32,33] 


where  x is  any  contour  surrounding  the  crack  tip  and  traversing  in  the  counterclockwise 
sense  as  shown  in  the  figure,  S is  the  arc  length  of  the  curve,  0 is  the  strain  energy 
density,  Tj  is  the  traction  vectors  defined  by  the  outward  normal  n along  x and  equals 
°ij  n,  Uj  is  the  displacement  vectors,  and  xj  is  the  axes  of  the  coordinate  as  defined  in 
Fig.  2.1. 

The  relation  between  the  J-integral  and  the  critical  stress  intensity  factor  for 
mode  I crack  propagation  is  [34] 


where  C = 1 and  (1  - d)  for  plane  stress  and  plane  strain,  respectively.  Because  the  J- 
integral  is  independent  of  the  chosen  path  [32,33],  it  can  be  evaluated  along  a path 
chosen  to  give  the  greatest  computational  advantage.  Thus,  either  the  stress  analysis  or 
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Fig.  2.1  Flat  surface  notch  in  two-dimensional  deformation 
field,  x denotes  any  curve  surrounding  the  tip  [32]. 
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energy  approach  can  be  used  to  analyze  fracture  of  materials,  depending  on  the 
problem  in  hand. 

2.2  Toughening  Mechanisms  in  Brittle  Matrix  Composites 

From  eqs.  (2.1)  and  (2.2),  it  is  clear  that  in  order  to  increase  the  failure  stress 
of  materials,  several  approaches  may  be  taken:  reduce  the  size  of  the  existing  critical- 
sized flaws,  or  increase  the  fracture  toughness  of  the  material,  or  introduce  a 
compressive  stress  on  the  surfaces  of  the  material.  In  practical  terms,  the  first  approach 
has  limitations.  Even  if  machining  flaws  and  inclusions  are  eliminated,  many  flaws  can 
be  introduced  on  the  surface  during  handling  and  in  service.  The  latters  can  be  achieved 
by  microalloying  the  inter  metal  lies  to  improve  their  plasticity,  or  by  tempering  and  ion- 
exchange  strengthening,  or  by  compositing  approaches. 

Any  composite  approach  will  utilize  one  or  more  of  the  following  means  to 
improve  the  toughness  of  brittle  materials: 

1 . Increase  the  fracture  surface  area. 

2.  Create  an  irreversible  process  zone  near  the  crack  tip  to  enhance  energy 
absorption. 

3.  Change  the  fracture  mode,  i.e.,  mode  I to  mode  II  or  III. 

4.  Reduce  the  stress  intensity  at  the  crack  tip. 

In  the  following  few  pages  the  mechanisms  utilizing  the  above  means  to 
improve  the  fracture  toughness  of  brittle  materials  will  be  reviewed.  The  mechanisms 
usually  fall  into  several  overlapping  categories  and  may  be  highly  interactive. 
However,  the  discussion  of  these  mechanisms  will  be  carried  out  in  separate  categories 
for  convenience. 

Crack  Deflection.  Crack  deflection  involves  out  of  plane  bending  (tilting  and 
twisting)  of  a crack  front.  The  process  is  pictured  schematically  in  Fig.  2.2.  The 
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Fig.  2.2  Caption  overleaf. 
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Fig:  2.2  Mechanisms  for  producing  crack  deflection,  (a)  Residual  stresses 

from  thermal  expansion  mismatches,  (b)  Microcracking,  leading  to 
crack  branching  and  deflection,  (c)  Elastic  mismatches  between  phases, 
leading  to  induced  stresses  and  crack  deflection,  (d)  Weak  interfaces, 
which  deflect  crack. 
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approaches  for  ensuring  crack-microstructure  interaction  and  thus  producing  crack 
deflection  include  the  introduction  of  residual  stresses  from  thermal  expansion 
anisotropy  and/or  thermal  expansion  mismatches,  Fig.  2.2(a);  stress  induced 
microcracking  at  second-phase  particles,  Fig.  2.2(b);  induced  stresses  from  elastic 
mismatches,  Fig.  2.2(c),  and  weak  interface  between  phases,  Fig.  2.2(d).  The 
toughness  increase  due  to  crack  deflection  comes  not  only  from  the  increase  in  the 
fracture  surface  area,  but  also  from  the  change  of  the  fracture  mode  from  pure  mode  I 
to  a mixture  of  mode  I and  mode  II  (tilted  crack)  or  mode  I and  mode  III  (twisted 
crack).  As  the  fracture  mode  changes  from  mode  I to  a mixed  mode,  the  primary 
driving  force  (i.e.,  the  Kj)  decreases,  and  the  crack  must  also  be  driven  by  the  shear 
components  of  the  stress  field.  In  addition,  the  typically  larger  values  of  Knc  and  Knic 
relative  to  Kic  are  also  contributed  to  the  increased  toughness. 

Effects  of  crack  deflection  on  material  toughness  have  been  examined 
theoretically  by  Faber  and  Evans  [35]  for  composites  containing  dispersions  of  spheres, 
discs  and  rods.  In  their  treatment,  the  local  stress  intensity  factors  at  crack  front,  Ki, 
K2  and  K3,  as  a function  of  deflection  angle  were  assessed  and  the  driving  force  for  the 
crack  advance  was  assumed  governed  by  the  average  strain  energy  release  rate  across 
the  crack  front  which  contained  planar  portion  as  well  as  tilted  and  twisted  portions.  Of 
particular  interest  in  their  results  is  the  prediction  that  toughening  is  invariant  with 
particle  size,  and  depends  only  on  particle  shape,  distribution  and  volume  fraction.  The 
most  effective  morphology  for  toughening  is  the  rod  with  a high  aspect  ratio.  As  high 
as  four-fold  increases  in  fracture  toughness  has  been  predicted  for  rods  with  an  aspect 
ratio  of  12.  Additional  important  results  from  their  analysis  are  that  toughening  comes 
almost  entirely  from  crack  twist  and  is  volume  fraction  invariant  above  ~ 0.2  for  all 
three  morphologies  studied. 

Experimentally,  several  investigators  have  observed  evidence  of  crack  deflection 
and  toughening  by  this  mechanism.  Faber  and  Evans  [36],  in  a companion  study  to  the 
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one  discussed  above,  verified  the  independence  of  deflection  toughening  on  particle  size 
and  the  morphology  effects  predicted  by  their  model.  They  also  observed  [37]  the 
existence  of  crack  deflection  toughening  in  SiC/Al203  and  SiC/ZrC>2  composites. 
About  60%  and  100%  increases  in  toughness  over  the  sintered  SiC  matrix  were 
reported  for  the  SiC/Al203  and  SiC/Zr02  composites,  respectively.  Swearengen  et  al. 
[38]  observed  deflection  toughening  in  composites  of  a low  expansion  borosilicate  glass 
containing  alumina  spheres  and  as  high  as  a three-fold  increase  in  toughness  was 
achieved.  Carter  and  Hurley  [39]  reported  that  crack  deflection  was  an  important 
toughening  mechanism  in  SiC-whisker-reinforced  MoSi2,  and  a toughness  increase  of 
50%  over  the  unreinforced  matrix  was  obtained.  All  in  all,  increases  of  2 - 4 times  in 
fracture  toughness  over  that  of  the  matrix  are  possible  by  the  crack  deflection 
mechanism. 

Crack  Bowing.  Crack  bowing  is  a process  of  in-plane  bending  of  a crack  front 
whereby  a propagating  crack  is  momentarily  pinned  at  discrete  locations  by  appropriate 
dispersed  second-phase  inhomogeneities.  Figure  2.3  is  a schematic  of  this  process, 
which  shows  that  a crack  front  interacts  with  an  array  of  second-phase  inclusions,  i.e., 
the  crack  front  bows  out  between  the  inclusions.  This  phenomenon  could  be  achieved 
by,  for  example,  dispersion  of  second  phase  particles,  fibers  or  voids. 

Direct  observation  of  crack  bowing  has  been  made  using  ultrasonic  fractography 
techniques  [38,40-42],  which  leaves  a permanent  record  of  crack  front  advancement  on 
the  fracture  surface  in  the  form  of  rib  mark.  Other  investigators  also  reported  crack 
bowing  produced  by  voids  in  NaCl  [43,44],  MgO  [45]  and  AI2O3  [46],  by  whiskers  or 
fibers  in  MgO  [47]  and  in  glasses  [48],  by  brittle  particles  in  a brittle  epoxy  [49]  and 
glass  [50]  and  by  ductile  particles  in  glasses  [41,42,51,52], 

Toughening  by  crack  bowing  is  attributed  to  a line  tension  effect  at  the  crack 
front,  a concept  first  postulated  by  Lange  [53].  In  his  and  other  subsequent  analyses 
[54-56],  it  was  assumed  that  the  increment  of  energy  absorbed  for  an  increment  of 
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Secondary  crack 


Fig.  2.3  Schematic  of  crack  bowing. 
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crack  extension  consisted  of  two  parts,  i.e.,  one  part  associated  with  the  energy  to  form 
new  surface  area  and  the  other  part  associated  with  the  energy  to  form  the  increased 
length  of  crack  front.  The  localized  pinning  of  cracks  would  result  in  an  increase  in  the 
energy  associated  with  the  crack  front  (line  tension  effect)  since  the  length  of  the  crack 
front  was  increased  due  to  crack  bowing.  Their  analyses  [54-56]  suggested  that  a 
substantial  increases  in  toughness  and  strength  could  be  achieved  by  a high  volume 
fraction  of  a dispersed  phase  of  small  size.  However,  a high  volume  fraction  of  small 
size  particles  is  difficult  to  achieve  because  of  uniformity  of  dispersion  and  chemical 
reactivity  problems.  In  addition  to  these  difficulties,  there  must  be  some  effective 
mechanisms  for  ensuring  that  the  crack  front  intersects  the  particles,  rather  than 
bypassing  them.  As  noted  in  the  discussion  of  crack  deflection,  this  would  require 
significant  residual  stresses.  With  this  condition  satisfied,  the  composites  would  also 
exhibit  crack  deflection  toughening  and  toughening  by  other  mechanisms  (for  instance, 
crack  bridging  and/or  crack  shielding),  and  the  contribution  of  the  crack  bowing  would 
probably  be  negligible.  Due  to  the  aforementioned  limitations,  only  about  20%  - 100% 
increase  by  crack  bowing  in  toughness  over  the  basic  matrix  have  been  reported. 

Crack  Shielding.  Crack  shielding,  as  it  is  usually  considered,  is  a toughening 
mechanism  involving  a process  zone  in  which  the  microstructure  changes  due  to  the 
stress  field  in  the  near  vicinity  of  crack  tip.  The  process  zone  shields  the  crack  from  the 
applied  stresses,  hence  reducing  the  stress  intensity  at  the  crack  tip.  Approaches  to 
achieve  crack  shielding  consist  of  (1)  phase  transformation  in  the  process  zone,  Fig. 
2.4(a),  such  as  martensitic  transformation  of  zirconia  in  alumina-zirconia  [57-59]  and 
mullite-zirconia  composites  [60-64];  and  (2)  stress-induced  microcracking  in  the 
process  zone,  Fig.  2.4(b),  as  observed  in  rock  [65],  AI2O3/Z1O2  composites  [66]  and 
ceramic-glass  materials  [67,68]. 

Transformation  toughening  in  ZrC>2  containing  materials  is  related  to  the 
diffusionless  transformation  of  metastable  tetragonal  ZrC>2  into  a stable  monoclinic 
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Fig.  2.4  Microstructural  features  of  ceramics  toughened  by  crack  shielding, 
(a)  crack  shielding  arising  from  phase  transformation;  (b)  crack 
shielding  resulting  from  microcracking. 
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phase  with  a higher  specific  volume.  This  increase  in  specific  volume  effectively 
reduces  the  stress  intensity  at  the  crack  tip  by  load  transfer  from  the  crack  tip  to  the 
near-by  region,  as  occurs  in  the  yield  region  at  the  front  of  a crack  tip  in  a metal.  Thus, 
additional  load  must  be  applied  to  cause  the  crack  to  propagate.  The  contribution  of  the 
transformation  to  toughness  has  been  treated  theoretically  by  several  investigators  [69- 
73].  In  all  analyses,  the  fracture  toughness  of  the  composites  was  treated  as  the  sum  of 
the  matrix  toughness  and  the  transformation  toughening  contribution.  The  contribution 
from  the  transformation  toughening  was  related  to  the  chemical  free-energy  change  for 
the  transformation  of  tetragonal  ZrC>2  to  monoclinic  ZrC>2,  to  the  change  in  strain 
energy  and  the  excess  surface  energy  associated  with  the  transformation,  and  to  the  size 
of  the  transformation  zone  and  volume  fraction  of  the  metastable  phase  in  the 
composites.  All  analyses  lead  to  the  following  conclusions:  toughness  can  be  increased 
by  (1)  decreasing  the  service  temperature,  (2)  decreasing  the  alloy  content  in  Z1O2,  (3) 
increasing  the  volume  fraction  of  the  metastable  phase,  and  (4)  increasing  the  size  of 
the  transformation  zone. 

Phase  transformation  toughening  has  been  observed  in  many  other  systems. 
Moya  and  Osendi  [60]  reported  an  increase  of  Kic  from  2 to  3 MPa.  m 1/2  for  a 
mullite/Zr02  composite.  For  Zr02-toughened  alumina,  Kic  values  of  ~ 12  [70]  and  16 
MPa. m 1/2  [71]  were  reported.  Transformation  toughening  was  also  found  in  single 
crystal  and  polycrystal  zirconia  materials  containing  metastable  precipitates  (partially 
stabilized  zirconia)  or  consisting  of  only  tetragonal  phase  (fully  stabilized  zirconia). 
Toughnesses  as  high  as  9 MPa.  m 1/2  have  been  reported  for  tetragonal  polycrystals 
[74],  while  ~ 12  MPa.m1/2  was  observed  in  MgO-partially-stabilized  zirconia  [75]. 
Transformation  toughening  to  date  is  one  of  the  most  successful  toughening  approaches 
used.  However,  it  is  highly  temperature  dependent.  As  temperature  increases, 
toughening  decreases,  going  to  zero  when  temperature  reaches  a point  where  the 
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decrease  of  Gibbs  free  energy  due  to  the  transformation  is  offset  by  the  increase  of 
strain  energy  associated  with  the  transformation. 

Microcracking  in  a process  zone  around  macrocrack  tips  is  a result  of  the 
combination  of  the  residual  stresses  and  the  stress  field  associated  with  the  crack. 
Opening  of  the  microcracks  requires  the  contribution  of  energy  from  the  applied  loads. 
Thus,  additional  load  must  be  applied  to  make  the  crack  propagate.  For  microcracking 
to  occur,  there  should  be  property  mismatches.  Mismatches  in  thermal  expansion  are 
the  most  common  source  [76],  but  mismatches  due  to  phase  transformation  [77-79]  and 
differences  in  elastic  modulus  [78]  can  also  be  important.  It  is  the  combination  of  such 
residual  stresses  with  the  applied  stress  field  that  causes  microcracking. 

The  contribution  of  microcracking  to  toughness  has  been  analyzed  theoretically 
by  Hoagland  and  Embury  [80],  Evans  and  Cannon  [71],  Evans  and  Faber  [81],  and 
Hutchinson  [82].  In  all  treatments,  microcrack  toughening  was  assessed  in  terms  of  the 
change  of  elastic  energy  associated  with  the  release  of  residual  stresses  on  the 
microcrack  plane  and  reduction  in  the  elastic  moduli  of  the  process  zone,  and  the 
increase  of  surface  energy  arising  from  the  formation  of  the  microcracks.  Their  results 
showed  that  the  net  increase  of  toughness  was  a function  of  the  microcrack  density 
(number/unit  volume),  the  size  of  the  process  zone,  the  contributions  from  the  residual 
stress  field  and  from  the  increase  of  the  compliance  of  the  process  zone,  and  the  surface 
area  of  the  microcracks.  Microcrack  toughening  is  strongly  particle  size  dependent 
[68,78,81,83].  At  relatively  small  particle  sizes,  particles  are  stable  even  in  the  near-tip 
field  and  therefore,  microcracking  does  not  occur.  As  the  particle  size  increases,  a zone 
of  microcracks  is  formed,  resulting  in  an  increase  in  toughness.  At  some  critical 
particle  size,  spontaneous  microcracking  occurs  on  cooling,  leading  to  no  toughening 
for  the  materials.  Thus,  toughening  only  occurs  with  median  size  of  particles. 

Experimentally,  microcrack  toughening  has  been  observed  in  many  systems. 
Claussen  et  al.  [84]  observed  an  increase  of  Kjc  from  ~ 5 to  10  MPa.m172  for  ZrC>2 
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toughened  AI2O3  composite.  Fu  and  Evans  [83]  reported  a toughness  increase  in  three 
single  phase  polycrystal  materials,  i.e.,  AI2O3,  BaTi03  and  Nb20s.  Although  it  is 
possible  to  utilize  microcracking  to  toughen  brittle  materials,  the  toughening  increments 
due  to  microcracking  are  relatively  small,  compared  with  transformation  toughening, 
because  transformation-related  dilatation  is  substantially  larger  than  microcracking 
dilatation.  Another  problem  with  microcrack  toughening  is  the  concern  of  the 
coalescence  of  microcracks  to  form  strength-limiting  macrocracks.  To  date,  there  is 
little  experimental  evidence  in  the  way  of  microcracking  effects  on  strength. 

An  additional  mechanism  which  may  be  classified  under  the  stress-induced 
microcracking  is  crack  branching.  A few  ceramics  have  been  demonstrated  to  have 
intricate  crack  bifurcations  [85].  The  crack  opening  associated  with  the  branched  cracks 
provides  volume  expansion  in  the  process  zone  and  should  result  in  toughening  [86]. 

Crack  Blunting.  Crack  blunting  occurs  when  the  sharp  crack  tip  is  blunted  via 
limited  plasticity  or  viscous  flow  at  the  crack  tip  or  when  the  sharp  crack  intercepts  a 
void  or  a ductile  phase  inclusion.  Blunting  of  the  crack  reduces  crack  tip  stress 
intensity,  and  therefore  improves  the  fracture  toughness. 

Viscous  flow  would  normally  be  associated  with  glassy  phases  and  high 
temperatures.  Thus,  the  toughening  effect  via  viscous  flow  is  very  limited  because  such 
flow  may  indeed  be  the  source  of  failure  through  gross  viscous  flow.  Similarly,  the 
effects  of  crack  tip  blunting  through  microplasticity  are  very  limited  because  of  the  high 
Peierl's  stresses  and  limited  slip  systems  available  in  typical  ceramics  and 
intermetallics.  Introduction  of  ductile  second  phases  should  provide  a seemingly 
effective  blunting  mechanism.  However,  other  toughening  mechanisms  (crack  bridging, 
crack  bowing  and  crack  deflection)  would  operate  with  the  introduction  of  the  ductile 
phases  and  the  contribution  of  the  crack  tip  blunting  process  would  probably  become 
negligible. 
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Incorporation  of  porosity  has  been  shown  to  improve  toughness.  An  example 
was  lead  zirconate  titanate  material  containing  an  array  of  ordered  voids,  which  showed 
1-2  fold  increases  in  the  work  of  fracture  and  100%  increase  in  strength  of  the  ceramic 
[87].  The  toughening  was  attributed  to  crack  bowing,  blunting  and  shielding,  the  latter 
through  the  local  changes  in  compliance  associated  with  the  voids.  Another  example  is 
provided  by  Passmore  et  al.  [50]  with  porosity  in  AI2O3. 

Crack  Bridging.  Crack  bridging  refers  to  a phenomenon  where  unbroken 
ligaments  behind  a crack's  wake  "bridge"  across  the  crack  faces  and  hold  them  together 
to  some  extent,  thus  reducing  the  stress  intensity  at  the  crack  tip.  These  ligaments  can 
be  coarse  grains  in  noncubic  monolithic  materials  such  as  alumina  [88,89],  a ductile 
metal  [20,21,90-95],  strong  whiskers  [96-101],  and/or  fibers  [102-104], 

The  effects  of  crack  bridging  on  material's  toughness  have  been  the  focus  of 
several  theoretical  treatments  [20-26,92,93,105-107],  In  these  analyses,  the  action  of 
the  bridging  ligaments  is  treated  in  terms  of  a zone  of  tractive  forces  in  the  crack  wake 
which  tend  to  hold  the  crack  faces  together  and  reduce  the  stress  intensity  seen  at  the 
crack  tip,  as  shown  in  Fig.  2.5.  Thus,  the  total  stress  intensity  acting  at  the  crack  tip  is 
assumed  to  be  a superposition  of  an  opening  mode  stress  intensity  caused  by  the 
external  stress  and  a stress  intensity  of  opposite  sign  caused  by  the  clamping  forces  of 
the  bridging  ligaments.  The  key  to  analyzing  the  toughening  effect  of  the  bridging  lies 
in  the  estimation  of  the  clamping  forces  which  in  turn  depends  on  the  nature  of  the 
ligaments,  i.  e.,  ductile  phases  or  ceramic  fibers,  and  the  characteristics  of  the 
matrix/reinforcement  interface,  viz.,  frictional  bridging  or  strong  particle  bridging.  For 
ceramic  reinforcements,  the  magnitude  of  the  clamping  forces  is  treated  according  to 
the  characteristics  of  the  interface  in  terms  of  frictional  bridging  [20,108]  or  strong 
bonding  bridging  [20],  or  simplified  spring  bridging  [25,106,107],  These  analyses  lead 
to  two  important  conclusions:  (1)  the  toughness  increases  as  the  strength  of  the 
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Fig.  2.5  Crack  bridging  by  ductile  ligaments.  Note  that  the  clamping  forces  of  die 
bridging  ligaments  are  a function  of  the  distance  away  from  the  crack  tip. 
The  length  of  the  arrows  indicates  the  magnitude  of  the  clamping  forces. 


Table  2.1  Ceramics  with  Enhanced  Toughness  [112] 


Toughening  Mechanism 

Material 

Maxium  Toughness 
MPa.m1/2 

LAS/SiC 

~ 20 

Fiber  reinforced 

Glass/C 

~ 20 

SiC/SiC 

~ 20 

Whisker  reinforced 

Al203/SiC  (0.2) 

10 

Si3N4/SiC  (0.2) 

14 

AI2Q3/AI  (0.2) 

12 

Ductile  network 

B4C/AI  (0.2) 

14 

WC/Co  (0.2) 

20 

PSZ 

18 

Transformation  toughened 

TZP 

16 

ZTA 

10 

Microcrack  toughened 

ZTA 

7 

Si3N4/SiC 

7 
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reinforcement  increases,  and  (2)  the  stress  required  to  propagate  a matrix  crack  is 
almost  independent  of  the  crack  length,  as  contrasted  with  the  behavior  of  cracks  in 
unreinforced  brittle  materials,  for  which  the  strength  decreases  with  a '1/2. 

Evidence  of  crack  bridging  by  ceramic  whiskers  or  fibers  were  reported  in 
several  composites.  A Kic  of  25  MPa.m1/2  was  achieved  in  SiC-fiber-reinforced 
ZrTi04  [109],  and  an  increase  of  Kjc  from  6 to  12  MPa.m1/2  was  observed  in  SiC- 
whisker-reinforced  zirconia  [110].  MoSi2  was  also  toughened  with  SiC  whiskers  and  an 
increase  of  50%  in  toughness  was  reported  [111].  An  important  conclusion  from  these 
experiments  is  that  a weak  interface  leads  to  a higher  toughness  than  a strong  interface. 
The  debonding  at  the  matrix/reinforcement  interface  guarantees  that  the  matrix  crack 
bypasses  the  fiber  or  whisker,  leaving  them  intact.  Otherwise,  the  crack  propagates 
through  the  reinforcement,  and  there  would  be  litter  improvement  in  toughness. 

Bridging  by  ductile  phases  also  provide  an  effective  means  to  improving 
toughness.  This  will  be  detailed  separately  in  the  next  section  because  the  topic  is  the 
focus  of  the  present  study. 

Finally,  it  would  be  useful  to  finish  this  section  using  a table  (Table  2.1)  from 
Evans  [112],  which  summarizes  the  maximum  toughness  achieved  using  different 
toughening  mechanisms.  From  the  table  and  the  above  discussion,  it  is  clear  that  crack 
bridging  and  phase  transformation  to  date  are  the  most  effective  mechanisms.  However, 
phase  transformation  toughening  is  limited  to  relatively  low  temperatures.  In  contrast, 
crack  bridging  can  provide  toughening  at  ambient  temperatures  as  well  as  at  elevated 
temperatures.  Thus,  it  appears  that  crack  bridging  is  the  most  important  toughening 
mechanism  at  elevated  temperatures. 

2.3  Ductile  Phase  Toughening 

The  primary  mechanism  responsible  for  ductile  phase  toughening  is  the 
bridging  of  ductile  ligaments  behind  the  advancing  crack  tip,  although  other  effects 
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such  as  crack  deflection  [46,113,114],  crack  bowing  [44,45,54,55]  and  crack  shielding 
[23]  also  contribute.  The  physical  mechanism  of  the  bridging  is  shown  in  Fig.  2.5. 
The  ductile  ligaments  spanning  the  advancing  crack  must  stretch  as  the  crack  opens 
until  they  fracture  or  decohere.  The  work  of  stretching  and  fracturing  the  bridging 
ligaments  contributes  to  the  toughness  of  the  composite.  When  the  size  of  the  bridging 
zone  in  the  wake  of  the  crack  tip  is  small  relative  to  the  crack  length  and  the  specimen 
dimension  ( i.e.,  small-scale  bridging  ),  the  steady-state  toughness  of  the  composite, 
Kcom,  can  be  related  to  the  fracture  energy  increase  due  to  the  bridging,  AG,  with  the 
aid  of  the  equation  (2.4)  by  taking  the  contour  t as  shown  in  Fig.  2.6.  The  final 
relation  takes  the  form  of  [25,1 15] 

( 1 - v 2 ) K 2 

v com  / com 


+ a g 


(2.6) 


where  v,  E and  K are  Poisson's  ratio,  the  elastic  modulus  and  fracture  toughness, 
respectively.  The  subscripts  'com'  and  'm'  refer  to  the  composite  and  matrix, 
respectively.  The  fracture  energy  increase,  AG,  is  related  to  the  work  of  stretching  and 
fracturing  ductile  ligaments  as  follows  [20,23] 

AG  = Vf  IJc  (u)  du  (2.7) 

where  a(u)  is  the  nominal  stress  carried  by  the  constrained  ductile  reinforcement  for  a 
given  crack  opening  u,  Vf  is  the  volume  fraction  of  the  reinforcement  and  u*  is  the 
crack  opening  at  the  point  when  the  ductile  reinforcement  fails.  The  integral  in  eq. 
(2.7)  is  commonly  referred  to  as  the  work  of  rupture  of  the  ductile  reinforcement. 

From  eqs.  (2.6)  and  (2.7),  it  is  evident  that  the  key  to  predict  the  increased 
fracture  toughness  is  to  find  the  stress-displacement  function,  o(u),  of  the  constrained 
ductile  phase.  Since  a(u)  is  different  from  that  measured  in  a simple  tensile  test, 
several  investigators  have  attempted  to  relate  o(u)  to  the  uniaxial  stress-strain 
properties  of  the  ductile  phase.  The  methods  used  included  a slip  line  field  analysis 
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Fig.  2.6  J -integral  path  for  bridging  toughening.  The  shaded  area  indicates  the 

bridging  zone.  Contour  T +^2+  ^3+  ^4  >where^l  *s^ar 
from  the  crack  tip, ^3  is  rigth  at  the  tip,  T 2 & ^4  are  4ie  upper  and 
lower  surfaces  of  the  bridging  zone. 
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[20,23],  finite  element  methods  [23,24],  spring  models  [25,106,107]  and  geometric 
models  [21,23,24,26].  These  analyses  conclude  that  the  work  of  rupture  increases 
with  increasing  strength  of  the  ductile  reinforcement  and  decohesion  at  the 
matrix/reinforcement  interface. 

Experimental  determination  of  the  o(u)  has  also  been  initiated  by  Ashby  et  al 
[26].  Their  test  specimen  consisted  of  one  ductile  reinforcement  (Pb)  imbedded  in  a 
precracked  brittle  matrix  (glass),  and  a tensile  test  is  conducted  to  evaluate  the  stress- 
displacement  function,  o(u).  They  found  that  the  work  of  rupture  was  enhanced  by 
limited  decohesion  at  the  interface.  They  also  found  that  the  work  of  rupture 
normalized  with  the  yield  strength  and  size  of  the  ductile  phase  was  constant  for 
different  sizes  of  lead  in  glass.  Other  investigators  have  also  conducted  similar  tests  on 
other  systems  such  as  Nb  in  TiAl  and  TiNb  in  TiAl  [116,  117],  and  Nb  in  MoSi2  [118, 
119-121].  Working  with  TiAl/TiNb  systems,  Deve  et  al.  [117]  found  that  whether  or 
not  extensive  decohesion  was  desirable  for  a high  work  of  rupture  depended  on  the 
work  hardening  capability  of  the  ductile  phase.  For  MoS^/Nb  system  [120,121],  it  was 
found  that  the  normalized  work  of  rupture  depended  on  the  size  of  the  niobium  rather 
than  being  a constant,  and  varied  from  5.3  to  2.8  when  the  half  thickness,  aQ,  of  the 
niobium  changed  from  0.5  to  0.125  mm.  These  results  indicated  that  the  effects  of 
decohesion  and  size  of  ductile  phase  on  the  toughness  were  complicated  and 
interdependent  with  other  material  properties. 

Evidence  of  ductile  phase  toughening  has  been  presented  by  several 
investigators.  Khaund  and  Nicholson  [46]  observed  a maximum  composite  toughness  of 
greater  than  fifty  times  that  of  the  glass  matrix  for  composites  containing  20  vol.  % of 
aluminum.  Krstic  et  al.  [122]  reported  a similar  increase  in  fracture  toughness  in  a 
similar  composite.  A KIC  of  about  5 MPa.  m 1/2  was  achieved  in  borosilicate  toughened 
with  Fe-Ni-Co  alloy  particles  [113,114].  Improvement  of  toughness  was  also  found  in 
composites  reinforced  with  continuous  ductile  phase  network  [95,123-125],  An  increase 
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of  Kic  from  - 1.4  to  4.5  MPa.m1/2  was  observed  for  FeO/Fe  network  [125]  and 
Zr02/Zr0  network  composites  [124].  In  WC/Co  network  composites,  a Kic  of  20 
MPa.m1/2  has  been  achieved  [95].  Toughening  of  TiAl  with  pancake-shaped  Nb  has 
raised  Kic  from  8 to  20  MPa.  m 1/2  [126]. 

Parameters  affecting  ductile  phase  toughening  include  (1)  the  morphology  of  the 
ductile  phase,  (2)  the  nature  of  the  matrix/reinforcement  interface,  (3)  size  of  the 
ductile  phase,  and  (4)  intrinsic  mechanical  properties  of  the  ductile  phase.  Thus  far,  the 
morphologies  studied  include  particles  [46,113,114,122],  continuous  boundary 
networks  [95,123-125],  fibers  [91,127],  and  pancakes  [126].  For  particle-reinforced 
composites,  some  systems  have  shown  an  increase  in  toughness  [46,122],  some  have 
not  [46,122,128].  A common  observation  in  these  studies  was  that  whenever  matrix 
cracks  were  attracted  to  ductile  particles,  a large  improvement  in  toughness  was 
obtained.  Otherwise,  little  improvement  was  gained.  In  contrast,  toughening  effect  was 
observed  in  all  continuous  network,  fiber,  and  pancake  reinforced  composites.  The 
reason  for  this  is  that  the  special  geometric  shapes  of  continuous  network,  fibers,  and 
pancake  guarantee  the  crack/ductile  phase  interactions,  therefore  the  inherent  toughness 
of  ductile  phase  is  utilized. 

Theoretical  analyses  of  the  stress  field  [122,129]  indicated  that  crack  could  be 
attracted  to  or  pass  ductile  particle  indiscriminately  only  when  the  following  criteria  are 
satisfied  simultaneously:  (1)  Young's  moduli  of  matrix  Em  and  ductile  particle  Ep  are 
the  same  or  Em  > Ep,  and  (2)  thermal  expansion  coefficients  of  matrix  am  and 
particle  (Xp  are  the  same  or  ctp  is  slightly  lower  than  am.  This  prediction,  which  has 
been  substantiated  by  experiments  [46,122,128],  provides  a guidance  for  the  designing 
of  particulate  reinforced  composites. 

It  seems  that  different  shapes  of  ductile  phases  also  require  different 
matrix/reinforcement  interfaces.  As  mentioned  before,  weak  interfaces  enhance 
decohesion  at  the  interface,  and  therefore  enhance  the  work  of  rupture  of  ductile  phases 
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if  the  crack  bridging  mechanism  operates  [23,24,26].  For  particulate  toughened 
composites,  the  scenario  may  be  different,  i.e.,  cracks  tends  to  travel  along  the  weak 
interface  instead  of  utilizing  the  inherent  toughness  of  the  ductile  particles,  and 
therefore  strong  interfaces  aid  the  operation  of  the  bridging  mechanism,  leading  to  the 
utilization  of  the  toughness  of  the  ductile  particles  [46,122]. 

Effects  of  the  size  of  ductile  phases  have  been  examined  experimentally  by 
several  investigators.  The  systems  examined  were  A^CtyNi  laminate  [130],  WC-Co/W 
particle  [131],  and  NiAl/Nb  fibers  [132],  In  the  size  ranges  tested,  all  of  the  results 
exhibited  that  large  size  ductile  phases  had  better  toughening  effects  than  smaller  ones. 
These  results  are  in  agreement  with  the  theoretical  analysis  [26]. 

So  far,  the  effect  of  the  intrinsic  properties  of  ductile  phases  has  not  been 
studied  experimentally  and  systematically.  However,  theoretical  analyses  indicate  that 
toughness  increases  with  increasing  strength  and  work  hardening  rate  of  ductile  phases 
[23,24,26,120,121]. 

Finally,  it  should  be  pointed  out  that  the  modeling  work  mentioned  above 
[20,21,23,24,26]  analyzed  the  effects  of  ductile  phase  constitutive  behavior  on 
toughness  quantitatively  based  on  the  assumptions  of  small  scale  bridging  of  an  infinite 
crack  in  an  infinite  specimen.  However,  in  cases  of  potential  toughening,  the  bridging 
zone  is  often  a significant  fraction  of  crack  length  [91,116,133]  and  toughening  is 
likely  to  be  sensitive  to  the  specimen  geometry  and  crack  length.  Recognizing  this 
situation,  Ravichandran  [134,135],  Cox  [136],  and  Cox  and  Low  [137]  have  recently 
proposed  models  describing  large  scale  bridging  to  take  the  effects  of  the  specimen 
geometry  and  crack  length  into  account.  In  their  treatments  [134,135],  it  is  assumed 
that  the  effective  stress  intensity  acting  at  the  crack  tip,  K^p,  is  a superimposition  of  an 
opening  mode  stress  intensity  caused  by  the  external  stress,  Kjnf,  and  a stress  intensity 
of  opposite  sign  caused  by  the  clamping  forces  of  the  bridging  ligaments,  R^dge*  i.e., 

Kfip  = Kjnf  - Kforidgg  (2.8) 
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By  using  appropriate  Green's  function,  Ktjp  can  be  related  to  the  external  stress,  Ojnf, 
and  the  bridging  stress,  ^bridge,  the  following  equation  [134,135] 

K up  = J*(  - ^bridge (x, a))  h(x,a)  dx  (2.9) 

where  a is  the  length  of  the  crack,  x is  the  position  along  the  crack  face,  and  h(x,a)  is 
the  Green's  function  which  is  a function  of  specimen  and  crack  configuration  and  can 
be  found  in  handbooks  for  some  standard  specimen  geometries.  Bridging  stress  is 
arrived  at  iteratively  [138].  However,  a knowledge  of  stress-displacement  function, 
a(u),  of  the  constrained  ductile  phase  is  necessary  for  the  iteration,  which  imposes  a 
demand  for  determining  o(u),  a situation  similar  to  the  small- scale-bridging  analysis.. 

2.4  Physical  and  Chemical  Properties  of  MoSio  matrix 
Molybdenum  disilicide  is  an  extremely  promising  intermetallic  compound  for 
high  temperature  structural  applications  in  oxidizing  atmospheres.  The  phase  diagram 
(Fig.  7)  [139]  indicates  that  MoSi2  melts  congruently  at  a high  temperature  of 
2020®C,  has  a tetragonal  structure  below  1900^C  and  a hexagonal  structure  above 
this  temperature.  It  has  a moderate  density  (6.31  g/cm^),  and  an  excellent  electrical 
and  thermal  conductivity.  It  is  extremely  brittle  with  a fracture  toughness  of  about  3 - 5 
MPa.  m 1/2  [91,111,140]  and  has  low  impact  strength  at  low  temperatures,  but 
exhibiting  a brittle  to  ductile  transition  at  approximately  925®C  [141,142].  Above  this 
temperature,  MoSi2  behaves  more  like  a metal  showing  yielding  and  stress-relieving 
characteristics.  Since  it  undergoes  creep  and  plastic  deformation  at  temperatures  above 
transition  temperature,  the  strength  at  elevated  temperatures  is  relatively  low. 
Therefore,  to  make  MoSi2  as  a viable  structural  material,  it  is  necessary  to  improve  the 
room-temperature  fracture  toughness,  elevated-temperature  strength  and  creep 


resistance. 
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Fig.  2.7  Mo- Si  Phase  Diagram  [139]. 
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MoSi2  has  excellent  high-temperature  oxidation  resistance.  Its  use  as  a heating 
element  in  furnaces  that  operate  in  air  at  temperatures  up  to  1700®C  is  a clear 
demonstration  of  its  excellent  high-temperature  properties.  Its  outstanding  oxidation 
resistance  is  due  to  the  formation  of  a glassy  silica  (SiC>2)  layer  which  acts  as  a 
protective  film  at  high  temperature  [143].  However,  MoSi2  can  disintegrate 
catastrophically  (pest)  during  oxidation  at  temperatures  about  400  - 600®C  [143-146]. 
At  these  temperatures,  a mixed  oxide  layer  containing  amorphous  SiC>2,  monoclinic 
M09O26,  and  M0O3  platelets  is  formed  on  the  surface  of  MoSi2,  as  contrasted  with 
the  formation  of  only  SiC>2  layer  at  higher  temperature  oxidation  [143].  The  mixed 
oxide  layer  is  voluminous,  poorly  adherent  and  even  flakes  off,  resulting  in  a complete 
disintegration  of  MoSi2  specimens  [144,146].  However,  a recent  research  [143] 
indicated  that  MoSi2-pest  could  be  prevented  by  ensuring  the  material  is  > 95  pet 
dense  or  by  high-temperature  preoxidation  to  produce  a protective  SiC>2  layer.  Thus, 
MoSi2-pest  is  not  expected  to  prevent  its  application  as  a monolithic  or  composite 
matrix  material  for  high-temperature  applications. 


CHAPTER  3 

EXPERIMENTAL  PROCEDURE 

3.1  Microstructural  Characterization  and  Raw  Materials 
The  microstructure,  fracture  surfaces,  and  interfaces  of  the  composites  were 
characterized  with  optical  microscope,  scanning  electron  microscope,  and  electron 
microprobe.  X-ray  diffractometry  (XRD)  was  used  to  assist  the  determination  of  the 
phases  presented  in  the  commercially  pure  MoSi2  powder.  The  particle  size  of  the 
MoSi2  powder  was  determined  based  on  the  principle  of  the  centrifugal 
photosedimentation  using  a particle  size  analyzer  (Model  CAPA-700).  The  grain  size, 
d,  of  the  hot-pressed  MoSi2  matrix  and  Nb  reinforcements  was  assessed  by  the  mean 
linear  intercept  method  using  the  following  formula: 


where  Lj  is  the  length  of  the  measuring  line  and  Pl  is  the  number  of  boundary 
intercepted  by  the  line. 

The  matrix  of  the  composites  was  commercially  pure  MoSi2  powder  of  -325 
mesh  (supplied  by  Johnson  Mathey  Inc.).  The  powder  had  an  irregular  shape  as  shown 
in  Figure  3.1.  Depending  on  the  lot  number  of  the  powder,  the  average  particle  size 
ranged  from  1.2  to  about  4.8  pm  diameter  determined  by  the  method  of  centrifugal 
photosedimentation.  The  composition  provided  by  the  vendor  was  Mo  > 61.0%,  Si 
36.6%,  C 0.06%,  N 0.03%,  O 1.22%  and  metallic  impurities  < 0.5%  (wt%).  Note 
that  there  is  a high  oxygen  content  in  the  as-received  powder.  Analysis  of  the  x-ray 
diffractometry  (XRD)  showed  only  MoSi2  peaks  (20  angles  from  5®  to  100^), 
indicating  other  phases  if  present  were  below  the  volume  fraction  detectable  by  XRD, 
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Table  3.1.  Mechanical  properties  of  as-received  niobium  * 


Thickness  of  Nb 
foil  (mm) 

1.0 

0.5 

0.25 

Microhardness 
(HV,  kg/mm2) 

110 

79.4 

82.3 

Yield  strength 
(oo,  MPa) 

204 

121 

123 

Tensile  strength 
(cu,  MPa) 

277 

195 

204 

Elongation  (8) 

0.403 

0.523 

0.468 

Strength  coefficient 
(MPa) 

454 

353 

370 

Strength  hardening 
coefficient 

0.173 

0.245 

0.243 

ao/HV 

1.85 

1.52 

1.49 

au/HV 

2.52 

2.45 

2.48 

* Only  one  sample  for  each  condition  was  tested. 


Fig.  3.1  Commercially  pure  MoSi2  powder  (supplied  by  Johnson  Mathey  Inc.) 
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or  in  the  form  of  amorphous  film  on  the  surface  of  MoSi2  powder  such  as  a thin 
amorphous  SiC>2  layer. 

Different  sizes  of  niobium  in  the  form  of  foil  or  filament  were  used  as  ductile 
phases  in  the  present  study  depending  on  the  experimental  objectives.  The  niobium  had 
a purity  of  99.8%  based  on  metallic  constituents  (supplied  by  Johnson  Mathey  Inc.). 
Tensile  and  hardness  tests  on  the  as-received  Nb  indicated  that  the  mechanical 
properties  of  Nb  depended  on  the  size  of  the  niobium  (see  Table  3.1). 

3.2  Fiber  and  Foil  Coating  Techniques 

The  coatings  applied  to  the  Nb  surfaces  included  ZrC>2  or  AI2O3  . The  former 
was  applied  by  dipping  the  metal  into  a coating  solution  consisting  of  97%  ZrC>2  and 
3%  AI2O3  (supplied  by  ZYP  Coatings,  Inc.),  while  the  latter  was  produced  via  sol-gel 
processing,  physical  vapor  deposition  or  hot  dipping  in  the  molten  aluminum  followed 
by  anodizing.  The  techniques  for  producing  AI2O3  coating  are  detailed  below. 

Sol-gel  coating  technique.  The  sol  solution  was  an  aluminum-alkoxide-derived 
sol  (using  aluminum-sec-butoxide  (ASB)),  hydrolyzed  in  excess  water  and  peptized 
with  aluminum  nitrate,  A1(NC>3)3.  The  procedure  for  preparing  the  sol  solution  has 
been  described  in  detail  by  Clark  et  al.  [147],  The  ratio  of  aluminum  nitrate  to  ASB  in 
the  sol  was  from  0.04  :1  to  0.07:1,  a range  which  has  been  reported  to  produce  the 
minimum  volume  gel  with  the  highest  solid  concentration  [148].  Solutions  with  ratios 
of  water  to  ASB  ranging  from  100:1  to  25:1  were  prepared.  The  compositions  of  the 
sols  used  are  summarized  in  Table  3.2.  Glycerol  was  added  to  some  sol  solutions  to 
improve  the  flexibility  of  the  dried  gel  [149]. 

The  sol-gel  coatings  were  produced  by  dipping  Nb  foils  into  the  sol  solution  for 
5 seconds  or  by  electrophoretic  deposition  onto  the  Nb  surface.  For  the  latter 
technique,  the  Nb  foil  was  kept  parallel  to  the  anode  with  a separation  distance  of  1 
cm.  A constant  DC  voltage  ranging  from  0.5  to  2 volts  was  applied  across  the 
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Table  3.2  Composition  of  the  sol  solution 


Designation 
of  the  sol 

A1(N03)3 

ASB 

H20  : ASB 

h2o 

glycerol 

PH 

Viscosity 

(centipoise) 

SG  1 

0.04  : 1 

100  : 1 

100  : 2 

3.35 

2.13 

SG2 

0.07  : 1 

100  : 1 

100  : 2 

3.05 

2.27 

SG3 

0.04  : 1 

50  : 1 

100  : 2 

3.05 

4.41 

SG4 

0.07  : 1 

50  : 1 

100  : 2 

2.85 

4.52 

SG5 

0.04  : 1 

25  : 1 

50  : 2 

2.80 

SG6 

0.04  : 1 

100  : 1 

nil 

3.35 

_ 

electrodes  and  the  coating  time  was  kept  at  30  minutes  for  all  the  experiments.  All 
coatings  were  conducted  within  24  hours  after  the  preparation  of  the  sol  to  avoid  aging 
of  the  sol.  Once  a coating  was  applied,  the  sample  was  withdrawn  from  the  sol  at  a rate 
of  10  cm/sec  for  both  dip  coating  and  electrophoretic  coating.  After  withdrawing  from 
the  sol,  the  Nb  foils  were  suspended  vertically  for  48  hours  at  the  ambient  temperature 
before  hot  pressing  with  MoSi2  powder.  Some  air  dried  Nb  foils  were  further  dried  at 
500®C  for  1 hr  with  a heating  rate  of  3®C/min  before  the  hot  pressing. 

Physical  vapor  deposition.  Physical  vapor  deposition  of  AI2O3  on  Nb  was 
conducted  using  a Thermionics  100-0030  system  with  an  electron  beam  heated  source, 
operating  in  2 x 10*6  torr  vacuum.  The  emission  current  was  200  mA  with  a potential 

difference  between  the  cathode  and  the  anode  being  3 KV,  resulting  in  a deposition  rate 

0 

of  approximately  10  A /sec.  By  changing  the  deposition  time,  three  different  deposit 
thicknesses  were  obtained,  i.e.,  1,  2 and  3 pm. 

Hot  dipping  and  anodizing.  Hot  dipping  of  Nb  foils  into  the  molten  aluminum 
bath  was  conducted  in  an  alumina  crucible  which  was  kept  inside  a box-type  resistance 
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furnace.  Aluminum  used  to  prepare  the  molten  bath  had  a nominal  purity  of  99.999 
wt%  based  on  metallic  constituents.  Two  different  bath  temperatures,  930®C  and 
870^C,  were  chosen  to  evaluate  the  effects  of  the  molten  metal  temperature  on  the 
thickness  and  uniformity  of  the  coating.  Dipping  time  varied  from  2 to  5 min 
depending  on  the  bath  temperature  used.  All  samples  were  withdrawn  from  the  bath  at 
a rate  of  10  cm  /sec.  Since  the  molten  bath  was  exposed  to  the  air  directly,  an  oxide 
layer  developed  on  the  surface  of  the  molten  bath,  which  caused  non-uniform 
aluminum  coating  on  the  first  few  Nb  foils.  These  were  excluded  from  the  subsequent 
processing. 

Anodizing  of  the  hot  dipped  Nb  foils  was  conducted  in  an  electrolyte  containing 
5 wt%  sulfuric  acid.  Such  low  acid  concentrations  were  used  because  high 
concentrations  tended  to  attack  the  alumina  formed,  and  the  film  produced  had  been 
reported  to  be  more  dense  using  a low  acid  concentration  [150].  The  process  was 
conducted  at  ambient  temperature,  and  a constant  DC  voltage  ranging  from  10  to  20 
volts  was  applied  to  produce  current  densities  ranging  from  5 to  25  mA/cm^. 
Anodizing  time  was  30  minutes  for  most  operating  conditions. 

3.3  Fabrication  of  MoSb  Matrix  Composites 

For  the  Nb  filament  toughened  composites,  the  filaments  were  pre-aligned  in  a 
plane  and  set  with  acrylic  resins.  The  Nb  filament  sheets  were  then  stacked  together 
with  MoSi2  powders  in  an  appropriate  thickness  ratio  and  vacuum  hot  pressed  at 
1700^C  for  40  min  under  a pressure  of  30  MPa.  The  diameter  of  the  filaments  was 
0.25  mm  and  the  volume  fraction  of  Nb  in  the  composites  was  10  or  20%  . 

For  the  Nb  foil  toughened  composites,  the  coated  or  uncoated  Nb  foils  were 
stacked  with  the  MoSi2  powders  and  vacuum  hot  pressed  at  1400^C  for  1 hr  with  a 
pressure  of  40  MPa.  The  thickness  of  the  Nb  foils  ranged  from  0.127  to  1.0  mm  and 
the  volume  fraction  of  the  niobium  was  20%  in  all  the  foil  toughened  composites.  In 
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order  to  minimize  residual  thermal  stresses,  the  hot  pressed  discs  were  held  in  the  hot 
pressing  chamber  at  80()0c  for  1 hour  before  cooling  down  to  room  temperature.  The 
residual  thermal  stresses  induced  with  such  processing  have  been  estimated  to  be  small 
with  a tensile  stress  of  19  MPa  in  the  matrix  and  a compressive  stress  of  79  MPa  in  the 
Nb  on  the  lamina  plane  [151]. 

3.4  Measurement  of  Intrinsic  Properties  of  the  Ductile  Phase 

Since  the  flow  behavior  of  constrained  ductile  phases  is  different  from  that 
measured  from  a simple  tensile  test,  the  experimental  procedures  proposed  by  Ashby  et 
al.  [26]  were  taken  to  measure  the  flow  behavior  of  constrained  Nb.  For  these 
experiments,  laminated  specimens  were  prepared  by  stacking  a Nb  foil  with  two  layers 
of  MoSi2  powder  at  an  appropriate  ratio,  and  then  vacuum  hot  pressing  at  1400^C  for 
1 hour  under  a pressure  of  40  MPa.  Three  different  thicknesses  of  Nb  foils  (1.0,  0.5 
and  0.25  mm)  were  used  to  prepare  the  laminated  composites.  The  hot  pressed  discs 
were  cut  into  rectangular  tensile  test  bars  with  dimensions  of  5.0  x 3.9  x 30.0  mm,  as 
shown  in  Fig.  3.2.  Straight-through  notches  in  the  MoSi2  matrix  were  introduced  using 
a diamond  wafering  blade  with  a thickness  of  0.15  mm.  The  distance  from  the  notch  tip 
to  the  laminate  interface  was  selected  as  300  pm  for  1.0  mm  thick  Nb  lamina  and  200 
pm  for  0.5  and  0.25  mm  thick  Nb  laminae.  Such  distances  have  been  found  to  be  small 
enough  not  to  distort  the  stress-displacement  curve  of  the  constrained  ductile 
reinforcement  and  at  the  same  time  maintain  an  appropriate  degree  of  elastic  constraint 
from  the  matrix  [119]. 

Displacement  controlled  tensile  tests  were  conducted  using  an  INSTRON 
machine  with  a displacement  speed  of  0.005  in/min.  Wedge  type  grips  were  used  to 
clamp  the  specimens  and  load  train  alignment  was  achieved  by  coupling  the  top  trip  to 
a universal  joint.  The  load  was  recorded  as  a function  of  cross-head  displacement. 
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Nb  MoSi2 


Fig.  3.2  Schematic  of  a composite  laminate 
tensile  test  specimen. 


Fig.  3.3  Short  bar  specimen  geometry. 
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3.5  Measurement  of  Fracture  Energy  of  Matrix/Reinforcement  Interface 
The  measurement  of  the  interfacial  fracture  energy  was  conducted  on  chevron 
notched  short  bars  using  the  procedures  recommended  in  Reference  152.  In  the  present 
study,  however,  a modified  specimen  geometry,  as  shown  in  Fig  3.3,  was  used.  An 
advantage  of  the  present  geometry  is  that  there  is  no  need  for  compliance  calibration 
due  to  symmetry  of  compliance  of  the  specimen  with  respect  to  the  interface.  The 
technique  involves  determining  the  critical  stress  intensity  factor  from  the  peak  load  of 
the  notched  bars  and  then  converting  the  critical  stress  intensity  factor  to  the  fracture 
energy  of  the  interface.  The  equation  used  to  calculate  the  fracture  energy  of  the 
interface,  Gmax>  is  [152] 

K 2 

r'  _ *MCSR  /O  -in 

Gmax  ~ „ (3-2) 

E 

where  E is  the  elastic  modulus  of  MoSi2  and  is  the  critical  stress  intensity 

factor  determined  from  the  peak  load  of  the  short  bar  test.  The  test  specimens  were  hot 
pressed  with  a coated  or  uncoated  Nb  foil  (0. 127  mm  thick)  sandwiched  in  between  two 
layers  of  MoSi2.  The  hot  pressed  discs  were  then  cut  into  rectangular  short  bar  with 
dimensions  of  8.28  x 9.25  x 14.29  mm  (Fig.  3.3).  The  notch  in  the  short  bar  was  cut 
parallel  to  the  foil  using  a diamond  wafering  blade  with  a thickness  of  0.4  mm.  The 
thickness  of  the  notch  was  chosen  to  be  slightly  larger  than  that  of  the  foil  to  insure  that 
the  initiation  and  propagation  of  the  crack  is  at  the  weakest  position  among  various 
interfaces  between  matrix/coating/reinforcement.  The  specimen  was  loaded  at  a 
constant  test  machine  crosshead  speed  of  0.005  in/min  using  a hydro- servo  controlled 
MTS.  The  load  was  recorded  as  a function  of  the  crosshead  displacement. 

3.6  Measurement  of  Toughness  of  the  Composites 
The  toughnesses  of  the  MoSi2  matrix  and  the  laminated  composites  were 
measured  by  4-point  bending  of  chevron  notched  specimens,  with  an  inner  and  outer 
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span  of  10  and  20  mm,  respectively,  using  a hydro-servo  controlled  MTS  with  a cross 
head  speed  of  4 x 10"4  mm/sec.  A schematic  of  the  bend  test  setup  is  shown  in  Fig. 
3.4.  The  notch  on  each  sample  was  cut  perpendicular  or  transverse  to  the  laminate 
plane  depending  on  the  purpose  of  the  experiments. 

The  peak  load  of  the  bending  tests  was  used  to  calculate  the  toughness  of  the 
composites,  Kmax,  with  the  aid  of  the  following  equation  [153] 


K 


max 


max  y 

B VW  ”” 


(3.3) 


where  Pmax  is  the  maximum  test  load,  B and  W are  the  width  and  height  of  the 
bending  bar,  respectively,  and  Ymjn*  is  the  minimum  value  of  the  dimensionless  stress 
intensity  factor  coefficient  as  a function  of  relative  crack  length  for  the  particular 
specimen  used.  It  should  be  noted  that  because  of  the  rising  crack-growth  resistance  for 
ductile-phase-toughened  composites,  the  maximum  load  and  the  minimum  value  of  the 
dimensionless  stress  intensity  factor  coefficient  do  not  occur  coincidentally  at  the  same 
crack  length,  and  therefore  Pmax  does  not  exactly  correspond  to  the  stress-intensity 
factor  at  failure,  but  a good  approximation  to  it  [154-157].  The  fact  that  Pmax  and 
Ymin  do  not  occur  coincidentally  at  the  same  crack  length  also  induces  a specimen 
size  effect  [154,157].  Thus,  the  dimensions  of  the  specimens  in  the  present  study  were 
kept  constant,  and  the  value  calculated  using  eq.  (3.2)  is  called  "damage  tolerance"  and 
designated  as  Kmax  in  the  paper,  rather  than  Kj£. 
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Fig.  3.4  Configuration  of  four  point  bend  test  on  chevron 
notched  specimens  (all  dimensions  in  millimeters). 


CHAPTER  4 

RESULTS  AND  DISCUSSION 
4.1  Microstructures  of  the  Composites 

4.1.1  Monolithic  MoSn 

Typical  microstructures  of  as-polished  MoSi2,  hot  pressed  at  temperatures 
between  1300-1700®C  under  a pressure  of  40  MPa  for  1 h,  are  shown  in  Fig.  4.1.  It  is 
evident  that  the  hot  pressed  monolithic  MoSi2  contains  second  phases  and/or  porosities 
as  indicated  by  the  dark  spots  in  Fig.  4.1.  Furthermore,  the  size  of  the  second  phases 
or  porosities  increase  with  increasing  processing  temperature.  High  magnification 
backscattered-electron  images  of  the  microstructures  processed  at  1400  and  1700^C 
(Fig.  4.2)  suggests  that  the  dark  spots  in  Fig.  4.1  fit  into  at  least  two  categories:  one  is 
second  phase  particles  which  are  solid  dark  areas  in  the  figure,  and  the  other  may  be 
porosities  or  second  phases  containing  high  Z elements  which  are  dark  areas  with  white 
points  in  them.  The  white  points  in  the  dark  areas  may  come  from  high  Z elements  or 
from  more  secondary  electron  signals  due  to  the  rough  surfaces  of  the  porosities.  Note 
that  the  MoSi2  hot  pressed  at  1 400^C  has  fewer  pores  or  second  phases  containing  high 
Z elements.  Furthermore,  it  is  noted  that  most  of  the  second  phases  and  porosities  are 
located  at  the  grain  boundaries  and  triple  points.  For  particle  sizes  larger  than  1 pm, 
electron  microprobe  analysis  indicated  that  the  majority  of  these  second  phase  particles 
were  SiC>2  and  others  were  Mo5Si3.  The  presence  of  Mo5Si3  and  SiC>2  in  the  hot 
pressed  monolithic  MoSi2  has  also  been  confirmed  using  TEM  technique  by  other 
investigators  [158-161]. 

If  the  oxygen  impurity  in  the  commercial  MoSi2  powder  is  assumed  to  be  in  the 
silica  particles,  then  the  volume  fraction  of  silica  can  be  estimated  to  account  for  about 
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• 9% 


42 


Fig.  4.1  Typical  optical  microstructures  of  as-polished  MoSi2  hot  pressed  at  (a)  1700°C,  (b)  1600°C,  (c)  1400°C 

and  (d)  1300°C  with  a pressure  of  40  MPa  for  1 hour. 
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Fig.  4.2  Backscattered-electron  images  of  MoSi2  hot  pressed  at 
(a)  1400°C  and  (b)  1700°C. 


44 


Table  4.1  Densities  of  hot  pressed  monolithic  MoSi2  measured 
using  Archimedes'  principle 


Hot  pressing  temperature, 

Density 

Theoretical  density 

<°C> 

(g/cm3) 

(%)* 

1300 

5.92 

97.21 

1400 

6.09 

100.00 

1600 

5.89 

96.72 

1700 

5.84 

95.90 

* Assuming  the  presence  of  6 vol  % SiC>2. 


6 vol%.  Thus,  by  assuming  the  presence  of  6 vol%  SiC>2  the  densities  of  the  hot 
pressed  MoSi2  measured  using  Archimedes'  principle  may  be  calculated  and  the  results 
are  summarized  in  Table  4.1.  The  results  suggest  that  the  highest  density  is  achieved  by 
hot  pressing  at  1400^C,  consistent  with  the  microstructural  observations.  The  reason 
for  decreasing  density  as  processing  temperature  increases  is  not  yet  clear.  The  grain 
size  of  hot  pressed  MoSi2  does  not  seem  to  depend  strongly  on  the  processing 
temperature  in  the  temperature  range  studied  as  shown  in  Fig.  4.3.  Instead,  the  second 
phases  become  coarser  as  processing  temperature  increases.  Insensitivity  of  the  grain 
size  to  the  processing  temperature  is  believed  to  be  related  to  the  presence  of  the  second 
phases  located  at  the  grain  boundaries  and  triple  points  and  provide  a pinning  force  for 
grain  growth.  Since  processing  at  140(Pc  provides  the  highest  density,  most  of  the 
composites  in  the  present  study  were  consolidated  at  this  temperature. 

4.1.2  MoSio/Uncoated  Nb  Interface 

The  microstructure  and  compositional  profiles  of  a MoS^/uncoated  Nb 
interface  hot  pressed  at  1400^C  are  shown  in  Fig.  4.4.  As  seen  in  the  figure,  Nb  and 
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Fig.  4.3  Secondary-electron  images  of  MoSi2  hot  pressed  at  (a)  1400°C, 
(b)  1600°C  and  (c)  1700°C.  The  corresponding  density  of 
each  sample  is  listed  in  Table  4. 1 
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Fig.  4.4  Microstructure  and  composition  profile  across  the  interphase  in  the 
MoSi2/uncoated  Nb  composite  hot  pressed  at  1400°C. 


Fig.  4.5  Microstructure  of  (Mo,Nb)5Si3  at  the  MoSi2/(Nb,Mo)Si2  interface  in  the 
MoSi2/uncoated  Nb  composite  hot  pressed  at  1400°C. 
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MoSi2  reacted  extensively  during  the  processing,  forming  three  new  intermetallic 
compounds  at  the  interface.  The  resulting  interaction  region  consisted  of 
Nb  | |Nb5Si3|  |(Nb,Mo)Si2  | |(Mo,Nb)5Si3  1 1 MoSi2 
However,  it  was  found  that  (Mo,Nb)5Si3  existed  in  a discrete  form  rather  than  a 
continuous  layer.  Most  of  the  (Mo,Nb)5Si3  phase  was  formed  around  the  dark  areas  in 
Fig.  4.4  which  were  porosity  or  Si02  . A high  magnification  of  (Mo,Nb)5Si3 
microstructure  at  the  MoSi2/(Nb,Mo)Si2  interface  is  shown  in  Fig.  4.5.  The  thickness 
of  the  interaction  zone  formed  was  measured  to  be  about  30  pm  for  hot  pressing  at 
1400^C  for  1 h,  and  the  nature  of  the  reaction  products  were  found  to  depend  on  the 
processing  temperature  and  time.  For  example,  Fig.  4.6  shows  the  interfacial 
microstructure  and  the  corresponding  composition  profile  of  a MoSi2/Nb  sample  which 
was  hot  pressed  at  HOO^C  for  30  min.  At  this  processing  temperature,  only  one 
(Nb,Mo)5Si3  interphase  was  formed  instead  of  the  three  intermetallics,  Nb5Si3, 
(Nb,Mo)Si2  and  (Mo,Nb)5Si3,  formed  at  1400^C. 

Such  temperature  dependence  of  interphase  formation  can  be  explained  in  terms 
of  Mo-Si-Nb  ternary  phase  diagram  and  temperature  dependence  of  the  diffusivity  of 
each  elemental  species.  The  available  isopleth  through  MoSi2/Nb  [162]  shows  that 
there  are  two  four  phase  equilibria,  at  1730  and  1740^C,  associated  with  MoSi2/Nb 
vertical  section.  The  isopleth  [162]  also  shows  that  the  relevant  phase  fields  remain 
essentially  unchanged  below  the  two  four  phase  equilibria.  Therefore,  the  composition 
path  in  the  present  study  can  be  analyzed  on  the  basis  of  the  800^0  isothermal  section 
[163],  which  has  been  slightly  modified  by  Lu  et  at.  [164]  to  eliminate  the  Nft^Si  phase 
[165]  and  is  shown  in  Fig.  4.7.  The  composition  paths  for  hot  pressing  at  1400^C  and 
1700^C  are  also  shown  in  Fig.  4.7  with  line  A and  B,  respectively.  The  non-planar 
interface  of  Nb5Si3/(Nb,Mo)Si2  in  Fig.  4.4  (hot  pressed  at  1400®C)  suggests  that  the 
composition  path  crosses  the  corresponding  two  phase  region  by  cutting  across  the  tie 
lines.  On  the  other  hand,  the  sharp  interface  of  Nb/Nb5Si3  indicates  that  the 
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Fig.  4.6  Microstructure  and  composition  profile  across  the  interphase  in  the 
MoSi2/uncoated  Nb  composite  hot  pressed  at  1700^C. 


Si 


Fig.  4.7  Isothermal  section  of  the  Mo-Si-Nb  ternary  diagram  [163]  at  80(PC,  showing 
the  composition  path  of  Nb/MoSi2  systems.  Line  A is  for  1400^C  and  line  B 
for  1700°C. 


49 


composition  path  runs  along  the  tie  line  in  the  corresponding  two-phase  region. 
Therefore,  the  composition  path  for  MOO^C  is  shown  in  the  following  route: 

MoSi2  — > (Mo,Nb)Si2  + (Mo,Nb)5Si3  (discrete  form)  ->  (Mo,Nb)Si2  + 
(Nb,Mo)Si2  ->  (Nb,Mo)Si2  -»  (Nb,Mo)Si2  + (Nb,Mo)5Si3  (cutting  across  the  tie 
lines)  ->  Nb5Si3  + Nb  (along  the  tie  line  of  the  binary  system)  ->  Nb. 

Similarly,  based  on  the  microstructure  examination  and  composition  analysis, 
the  composition  path  for  hot  pressing  at  1700^C  (Fig.  4.6)  is 

MoSi2  ->  (Mo,Nb)Si2  + (Mo,Nb)5Si3  (cutting  across  the  tie  lines)  — > 
(Mo,Nb)5Si3  ->  Nb5Si3  + Nb  (along  the  tie  line  of  the  binary  system)  — » Nb. 

Since  the  isopleth  of  the  Mo-Nb-Si  system  [162]  shows  that  the  relevant  phase 
fields  remain  essentially  unchanged  below  1730®C,  the  difference  in  the 
microstructures  at  different  processing  temperatures  below  1730^C  might  be  due  to  the 
changes  in  the  interface  reaction  rate  which  is  proportional  to  the  free  energy  of 
formation  of  these  compounds  or  due  to  changes  in  diffusion  kinetics  of  the  elemental 
species  to  the  reaction  sites.  According  to  the  phase  diagram  shown  in  Fig.  4.7,  the 
compounds  formed  in  the  MoS^/Nb  diffusion  couple  may  be  described,  depending  on 
whether  or  not  (Nb,Mo)Si2  forms,  by  the  following  formulae: 

MoSi2  + flNb->i  (Mox,  Nbj_x)5Si3  + c (Nbj.y,  MOy)Si2  (4.1) 

or 

MoSi2  + [(l-z)/z]  Nb  — >(l/5z)  (Moz,  Nb^S^  (4.2) 

where  the  coefficients  a,  b,  c,  x,  y and  z depend  on  the  composition  and  relative 
amounts  of  the  products.  For  the  above  equations,  it  is  assumed  that  1 mole  of  MoSi2 
reacts  completely  with  the  required  amount  of  niobium  to  form  the  compounds  listed  on 
the  right-hand  side  of  the  equations.  The  free  energy  for  reaction  (4.1)  can  be 
calculated  by  considering  several  intermediate  steps,  such  as  the  following: 

MoSi2  + Nb  = (4/7)  NbSi2  + (1/5)  Mo5Si3  + (3/35)  Nb5Si3  (4.3) 
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x Mo5Si3  + (1-x)  Nb5Si3  = (Mox,  Nbi_x)5Si3  (4.4) 

and 

y MoSi2  + (1-y)  NbSi2  = ( Nbi-y,  Moy)Si2  (4.5) 

Similarly,  equation  (4.2)  can  be  broken  down  into 

MoSi2  + (7/3)  Nb  = (7/15)  Nb5Si3  + (3/15)  Mo5Si3  (4.6) 

and 

(1-z)  Nb5Si3  + z Mo5Si3  = (Moz,  Nb^S^  (4.7) 


Reactions  (4.4),  (4.5)  and  (4.7)  correspond  to  the  formation  of  the  solid 
solutions  from  mixing  of  the  pure  components.  Based  on  the  available  thermodynamic 
data  [166],  free  energies  of  formation  for  consuming  1 mole  of  MoSi2  in  reactions 
(4.3)  and  (4.6)  have  been  calculated  for  two  different  temperatures,  1400^C  and 
1700°C,  and  the  results  are  listed  in  Table  4.2.  As  seen  in  the  table,  the  free  energies 
of  formation  for  reactions  (4.3)  and  (4.6)  do  not  vary  substantially  as  the  hot  pressing 
temperature  increases  from  1400  to  1 700^(2.  It  should  be  noted  that  depending  on  the 
amounts  of  MoSi2  and  niobium  participating  in  reaction  (4.3),  the  ratio  of  the  products 
NbSi2  to  Nb5Si3  may  not  be  fixed  as  indicated  in  the  equation.  Other  possibilities  such 
as  the  formation  of  much  more  Nb5Si3  than  NbSi2  or  a similar  amount  have  also  been 
considered.  However,  the  calculations  show  that  the  same  conclusion  mentioned  above 
(for  the  case  of  more  NbSi2  than  Nb5Si3)  still  holds  for  these  cases.  Since  no  data  on 
the  free  energy  of  mixing  for  reactions  (4.4)  and  (4.7)  are  available,  the  free  energy  of 
mixing  can  be  assumed  to  be  ideal.  This  may  not  be  an  unreasonable  approximation 
since  (MoxNbi_x,)5Si3  is  a continuous  solid  solution  as  shown  in  Fig.  4.7.  For  this 
case,  the  free  energy  of  mixing  is  proportional  to  temperature  for  a fixed  composition. 
Thus,  the  change  in  the  free  energy  of  mixing  for  reactions  (4.4)  and  (4.7)  from 
1400^C  to  1700°C  should  not  change  the  order  of  the  free  energy  of  formation  for 
reactions  (4.1)  and  (4.2).  A similar  argument  may  apply  to  reaction  (4.5),  since  the 
isopleth  of  the  Mo-Nb-Si  system  [162]  shows  that  the  relevant  phase  fields  remain 
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Table  4.2  The  standard  free  energy  of  formation  for  some  possible  reactions  in  the 

MoSi2/Nb  diffusion  couple 


Reaction 

AG  1400 

(kcal.mol-1) 

AG  1700 
(kcal.mol'1) 

4.3 

-44.46 

-49.76 

4.6 

-37.45 

-38.20 

essentially  unchanged  below  1730^C,  i.e.,  the  phase  relations  between  (Nb,Mo)Si2, 
(Mo,Nb)Si2  and  (Nb,Mo)5Si3  remain  unchanged  (Fig.  4.7),  implying  that  (Nb^.y, 
Moy)Si2  behaves  more  or  less  like  an  ideal  solid  solution.  Thus,  although  the  free 
energies  of  formation  for  reactions  (4.1)  and  (4.2)  change  with  temperature,  the  order 
of  the  free  energy  of  formation  for  reactions  (4.1)  and  (4.2)  does  not  vary  with 
temperature  in  the  range  of  interest.  Therefore,  the  free  energy  of  formation  of  these 
compounds  on  the  reaction  kinetics  may  not  explain  the  observed  microstructural 
change  from  1400®C  to  1700^C. 

In  contrast,  the  microstructural  changes  in  the  interaction  zone  at  different 
processing  temperatures  could  be  due  to  the  change  in  the  diffusion  kinetics  of  Si,  Nb 
and  Mo,  as  will  be  discussed  below.  The  composition  profile  in  Fig.  4.4  (1400^C) 
shows  clearly  that  Si  diffuses  far  deeper  than  Mo  into  the  Nb-rich  regions,  indicating 
that  Si  has  a higher  mobility  than  Mo.  This  is  consistent  with  the  available  literature 
data  [167,168],  which  show  that  diffusivities  of  Si  and  Mo  in  Nb  could  be  expressed  as 
follows. 

Dsi  =0.51xl0-2  exp(-48,000/RT) 
dMo  = 92  exp(-122,000/RT) 


700  - 1500°C 
1725  - 21820c 
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The  above  equations  indicate  that  the  diffusion  coefficients  of  silicon  in  niobium 
at  1400^C  and  1700^C  are  2.73  x 10"9  and  2.45  x 10"^  cm^s'l  respectively.  For 
molybdenum,  the  corresponding  values  are  1.56  x 10'^  and  2.81  x lO'^  cm^s'l. 
Because  of  the  low  mobility  of  Mo  relative  to  Si,  the  interphase  formation  at  1400®C  is 
dominated  by  the  diffusion  of  Si  into  Nb.  This  is  supported  by  the  formation  of  M>5Si3 
adjacent  to  Nb  instead  of  the  formation  of  (Nb,Mo)5Si3  However,  as  temperature 
increases,  diffusion  of  Mo  becomes  more  pronounced  relative  to  diffusion  of  Si  because 
of  its  higher  diffusion  activation  energy,  and  therefore  a stronger  temperature 
dependence  of  its  diffusivity.  Similarly,  Nb  diffusion  becomes  more  competitive  at 
high  temperatures,  which  combined  with  the  diffusion  of  Mo  leads  to  the  formation  of 
(Nb,Mo)5Si3  at  1700°C,  as  shown  in  Fig.  4.6.  The  limited  growth  of  (Mo,Nb)Si2 
interphase  at  1700^C  is  probably  due  to  the  faster  growth  of  (Mo,Nb)5Si3  caused  by 
the  more  pronounced  diffusion  of  Nb  relative  to  Si  at  high  temperature.  It  should  be 
noted  that  at  1400^C,  growth  of  Nb5Si3  depends  mainly  on  the  transport  of  Si  through 
the  (Nb,Mo)Si2  layer.  However,  at  1700^C,  as  Nb  diffusion  becomes  more 
pronounced,  the  growth  of  (Mo,Nb)5Si3  is  enhanced  not  only  in  the  direction  of  the 
Nb-rich  region  but  also  in  the  direction  of  the  (Nb,Mo)Si2  interphase.  Thus,  when  the 
diffusion  of  Nb  becomes  competitive  with  respect  to  the  diffusion  of  Si,  the  growth  of 
(Mo,Nb)5Si3  is  enhanced  and  that  of  (Nb,Mo)Si2  is  limited. 

Note  that  according  to  the  conservation  of  mass,  much  more  (Mo,Nb)5Si3 
should  exist  in  the  diffusion  couple.  However,  only  small  amount  of  (Mo,Nb)5Si3  has 
been  observed,  as  shown  in  Fig.  4.4.  This  is  believed  to  be  due  to  the  presence  of 
SiC>2  in  the  commercial  MoSi2  powder  as  shown  in  Fig.  4.4  and  reported  in 
References  158  - 161.  The  presence  of  SiC>2  makes  the  actual  composition  of  the 
starting  MoSi2  powder  shift  to  the  MoSi2*Si  two  phase  field,  instead  of  in  the  single 
phase  MoSi2  field.  The  presence  of  the  extra  silicon  is  also  reflected  in  the  composition 
profile.  As  shown  in  Fig  4.4,  Matano  interfaces  of  three  elements  are  not  at  the  same 
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position,  with  the  Si  Matano  interface  far  more  closer  to  the  Nb-rich  side  than  that  of 
Mo.  However,  the  difference  in  the  position  of  Matano  interfaces  becomes  smaller  at 
1700°C,  as  shown  in  Fig.  4.6.  The  results  suggest  that  the  presence  of  SiC>2 
contributes  more  to  the  diffusion  process  at  lower  temperatures,  whereas  at  high 
temperatures  contributions  from  MoSi2  and  Nb  become  more  pronounced. 

4.1.3  Chemical  Compatibility  of  Coating  Materials  with  Nb  and  MoSb 

The  chemical  compatibilities  of  the  potential  coating  materials  (AI2O3,  Z1O2 
and  mullite)  with  MoSi2  matrix  and  Nb  reinforcement  were  assessed  experimentally  via 
microstructural  and  compositional  analysis  of  the  samples  after  hot  pressing  MoSi2  and 
Nb  with  the  coating  materials.  Figure  4.8  shows  the  micrograph  and  compositional 
profiles  of  the  Al203/MoSi2  interface.  The  sample  was  produced  by  hot  pressing 
powder  mixtures  of  MoSi2  and  AI2O3  at  1700^C  for  40  minutes  followed  by  annealing 
at  1400°C  for  100  h.  As  seen  in  the  figure,  there  exist  two  very  distinctive  types  of 
interface:  a relatively  sharp  interface  with  a thickness  less  than  the  distance  between 
two  consecutive  scanning  points  of  the  electron  microprobe  (2  pm),  and  a relatively 
wide  interface  with  a gradual  change  in  the  composition  which  is  about  7 pm  thick. 
The  wide  interface  was  determined  to  be  rich  in  Si,  A1  and  O,  with  little  or  no  Mo,  as 
shown  in  Fig.  4.8.  As  discussed  in  Section  4.1.1,  the  commercial  MoSi2  contains 
SiC>2.  Consequently,  it  is  believed  that  the  interaction  layer  at  the  wide  interface  is  due 
to  a reaction  between  AI2O3  and  SiC>2  impurity  when  both  are  in  contact.  Figure  4.9 
which  shows  that  there  also  exist  sharp  interfaces  on  both  sides  of  an  AI2O3  particle 
further  supports  this  argument.  The  results  indicate  that  AI2O3  is  chemically 
compatible  with  MoSi2  and  that  the  presence  of  SiC>2  in  commercial  MoSi2  powder 
may  cause  the  degradation  of  the  AI2O3  coating  in  the  MoS^/Nb  composites. 

Similar  results  were  obtained  with  the  samples  of  MoSi2-Zr02  and  MoSi2- 
mullite.  For  instance,  the  interface  of  MoS ^/mullite  is  very  sharp  as  shown  in  Fig. 
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Fig.  4.8  Microstructure  and  composition  profile  across  an  AI2O3  particle  in  MoSi2  matrix 
hot  pressed  at  1700°C  for  40  min  and  annealed  at  1400°C  for  100  hours. 


Fig.  4.9  Composition  profile  across  an  AI2O3  particle  in  MoSi2  matrix  with  the  same 
processing  condition  as  the  sample  in  Fig.  4.8. 
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4.10  and  its  thickness  was  measured  to  be  around  lpm  (the  resolution  limit  of  the 
microprobe).  These  results  indicate  that  AI2O3,  ZrC>2,  and  mullite  are  chemically 
compatible  with  MoSi2.  The  microstructures  and  the  corresponding  composition  profile 
of  AI2O3/M)  interface  are  shown  in  Fig.  4.11.  A very  sharp  interface  (about  1 pm) 
indicates  that  Nb  is  chemically  stable  with  AI2O3  up  to  about  1700®C. 

The  above  findings  were  also  confirmed  by  thermodynamic  calculations.  Using 
the  thermodynamic  data  available  [169],  the  free  energy  changes  in  possible  reactions 
among  the  coating  materials  (AI2O3  and  ZrC>2),  MoSi2  and  Nb  were  calculated  at 
1400K  and  2000K.  The  reaction  components  considered  in  these  calculations  included 
Mo,  MoSi2,  Mo5Si3,  M03SP  MoO(g),  M0O2,  M0O3O),  Si,  SiO(g),  SiC>2,  Al, 
A1203,  AlO(g),  Al20(g),  Zr,  Zr02,  ZrO(g),  Nb,  NbO,  Nb02,  Nb205,  ZrSi2,  ZrSi, 
Zr5Si3,  Zr2Si,  and  02(g)  and  all  the  compounds  were  assumed  to  be  stoichiometric. 
Table  4.3  summarizes  some  reactions  which  have  lower  standard  free  energy  changes 
among  the  possible  reactions.  It  is  clear  from  Table  4.3  that  AI2O3  and  Zr02  are 
chemically  stable  with  MoS^/Nb  system  within  the  temperature  range  considered. 

However,  it  was  found  that  mullite  is  not  chemically  compatible  with  Nb,  as 
shown  in  Fig.  4.12.  Severe  interaction  was  observed  and  the  interaction  zone  was 
measured  to  be  about  10  pm  at  the  as-hot  pressed  condition.  The  reaction  products, 
however,  could  not  be  determined  because  of  the  fine  microstructure  (smaller  than  1 
pm)  in  the  interaction  zone  as  shown  in  Fig.  4.12.  Nevertheless,  the  result  clearly 
shows  that  mullite  could  not  be  used  as  a diffusion-barrier  coating  in  MoSi2/Nb 
composite. 

4.1.4  Effectiveness  of  the  Inert  Coating  as  a Diffusion  Barrier 

The  microstructure  and  composition  profiles  of  the  ZrC>2  uncoated  Nb/MoSi2 
interface  which  was  vacuum  hot  pressed  at  1700°C  for  40  minutes  are  shown  in  Fig. 

4. 13.  It  is  noted  that  the  thickness  of  the  interaction  zone  between  Nb  and  MoSi2  has 
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Distance  (microns) 


Fig.  4.10  Microstructure  and  composition  profile  across  an  mullite/MoSi2  interface  hot 
pressed  at  17()(PC  for  40  min. 


Fig.  4.1 1 Microstructure  and  composition  profile  across  an  Nb/Al203  interface  hot  pressed 
at  1700°C  for  40  min. 
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Table  4.3  The  standard  free  energy  changes  for  some  possible  reactions  among  the 
coating  materials,  matrix  and  reinforcements. 


Reactions 

AG1400K 

AG2000k 

(kcal) 

(kcal) 

MoSi2  + 5/6  Zr02  = 5/6  Zr  + 2/6  Mo3Si  + 10/6  SiO(g) 

102.338 

62.021 

MoSi2  + 7/5  Zr02  = 7/5  Zr  + 1/5  Mo5Si3  + 7/5  Si02 

76.964 

74.266 

MoSi2  + 4 Zr02  = 2 Zr2Si  + Mo02  + 3 02(g) 

646.815 

569.275 

Nb  + 1/2  Zr02  = NbO  + 1/2  Zr 

29.28 

27.799 

Nb  + 5/4  Zr02  = 1/2  Nb205  + 5/4  Zr 

94.664 

89.586 

MoSi2  + 5/9  A1203  = 10/9  A1  + 1/3  Mo3Si  + 5/3  SiO(g) 

99.396 

54.509 

MoSi2  + 7/5A1203  = 14/5A10(g)  + l/5Mo5Si3  4-  7/5SiO(g) 

341.994 

237.263 

MoSi2  + 7/5A1203  = 7/5A120(g)  + l/5Mo5Si3  + 7/5Si02 

125.433 

86.435 

Nb  + 5/6  A1203  = 5/3  A1  + 1/2  Nb205 

90.096 

78.318 

Nb  + 1/2  A1203  = 1/2  A120(g)  + NbO 

46.590 

32.145 
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Distance,  micrometer 


Fig.  4.12  Microstructure  and  composition  profile  across  an  Nb/mullite  interface  hot 
pressed  at  1700°C  for  1 hour. 


Fig.  4.13  Microstructure  and  composition  profile  across  the  MoSi2/Zr02-coated  Nb 
interface  hot  pressed  at  1700°C  for  40  min. 
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been  reduced  from  about  30  pm  for  the  uncoated  Nb  (Fig.  4.6)  to  4 urn  upon  coating 
of  Nb  with  ZrC>2.  In  addition,  the  reaction  product  at  the  Nb/Zr02  interface  was 
determined  to  be  Nb5Si3,  instead  of  (Mo,Nb)5Si3  for  the  uncoated  samples. 
Furthermore,  no  appreciable  diffusion  of  Nb  and  Mo  across  the  coating  was  detected, 
as  indicated  by  the  Nb  and  Mo  composition  profiles  in  Fig.  4.13. 

The  growth  rate  of  the  interaction  layer  (Nb5Si3)  in  the  ZrC>2  coated  Nb 
filament  reinforced  composites  were  determined  by  annealing  the  samples  for  25  and 
100  hours  at  1300®C.  As  shown  in  Fig.  4.14,  the  thickness  of  Nb5Si3  was  increased  to 
6 pm  after  annealing  for  100  hours.  The  thickness  was  also  found  to  fit  into  a parabolic 
function  of  the  annealing  time.  The  results  indicate  that  the  growth  of  the  Nb5Si3  layer 
is  controlled  by  the  diffusion  of  the  elemental  species  through  the  Nb5Si3  rather  than 
the  diffusion  of  Si  through  ZrC>2.  In  other  words,  the  diffusion  of  Si  is  relatively  fast 
through  ZrC>2.  Thus,  the  ZrC>2  coating  may  not  be  an  effective  diffusion  barrier  for 
Nb/MoSi2  system  because  of  the  high  Si  diffusivity,  even  though  it  is  chemically 
compatible  with  Nb  and  MoSi2- 

Similar  observations  were  made  for  the  AI2O3  sol-gel  coated  Nb/MoSi2 
interface.  Typical  composition  profiles  and  microstructure  at  the  Nb/Al203/MoSi2 
interface  produced  by  hot  pressing  at  1400®C  for  1 h are  shown  in  Fig.  4.15.  Similar 
to  the  ZrC>2  coated  Nb,  the  formation  of  Nb5Si3  at  the  Nb  side  of  the  Nb/Al203 
interface  was  observed.  Its  thickness  was  measured  to  be  around  15  pm,  compared 
with  30  pm  of  (Mo,Nb)5Si3  formed  with  uncoated  Nb  in  composites  produced  under 
the  same  hot  pressing  conditions.  These  experimental  observations  indicate  that  the 
diffusivities  of  Si  in  ZrC>2  and  AI2O3  must  be  relatively  high  compared  with 
diffusivities  of  other  species. 

An  order  of  magnitude  calculation  of  Si  diffusivity  through  Zr(>2  and  AI2O3 
coatings  made  from  the  thickness  of  the  intermetallics  formed  indicated  that  the 
apparent  diffusivity  of  Si  was  larger  than  1.0  xlO’9  cm2/sec  in  ZrC>2  at  1300°C  and 
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Fig.  4.14  The  thickness  of  Nb5Si3  vs  the  square  root  of  the  annealing  time  in  MoSi2/Zr02- 
coated  Nb  composites. 


Distance,  micrometer 


Fig.  4.15  Microstructure  and  composition  profile  across  the  MoSi2/Al203-coated  Nb 
interface  hot  pressed  at  140(PC  for  1 hour. 
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larger  than  2.  x 10"^  cm^/sec  in  AI2O3  at  1400^C.  These  values  are  considerably 
higher  than  the  bulk  diffusion  of  Si  in  these  ceramic  materials,  suggesting  that  Si 
diffusion  might  have  occurred  via  different  diffusion  paths  than  bulk  diffusion.  The 
most  feasible  fast  diffusion  route  may  be  through  the  grain  boundary.  Since  the  grain 
size  of  the  coating  produced  by  the  sol-gel  technique  is  small,  the  diffusion  flux  from 
grain  boundary  diffusion  may  dominate  the  diffusion  process  because  of  the  large  grain 
boundary  area.  Another  potential  fast  diffusion  path  is  the  free  surface  generated  by  the 
porosities  in  the  coating.  The  porosities  may  also  provide  a probable  path  for  gaseous 
diffusion  of  SiO.  Therefore,  for  these  coating  materials  to  be  used  as  diffusional 
barriers,  grain  boundary  area  and  processing  defects  must  be  minimized.  In  the  next 
section,  it  will  be  shown  that  the  effectiveness  of  the  coating  as  a diffusion  barrier 
indeed  depends  strongly  on  the  processing  techniques  and  the  process  control 
parameters. 


4.2  Optimization  of  Coating  Processes 

From  Section  4.1,  it  is  clear  that  AI2O3  and  ZK>2  are  chemically  compatible 
with  both  MoSi2  and  Nb.  However,  Z1O2  has  a much  higher  thermal  expansion 
coefficient  than  MoSi2  and  Nb.  This  may  cause  a very  low  interfacial  fracture  energy 
owing  to  the  thermal  expansion  mismatch  between  the  coating  and  the  matrix  and 
reinforcement  [170].  In  contrast,  AI2O3  has  a thermal  expansion  coefficient  close  to 
those  of  both  MoSi2  and  Nb.  Thus,  it  has  been  chosen  as  the  coating  material  for  the 
present  composite  system.  In  the  following  sections,  three  different  coating  techniques 
for  applying  a thin  AI2O3  film  on  Nb  reinforcement  will  be  discussed.  The  techniques 
consist  of  (i)  sol-gel  coating;  (ii)  physical  vapor  deposition;  (iii)  hot  dipping  the  metal 
in  molten  Al,  followed  by  anodizing  A1  to  form  AI2O3.  It  will  be  shown  that  the 
effectiveness  of  the  coating  as  a diffusion  barrier  depends  strongly  on  the  processing 
techniques  and  the  process  control  parameters. 
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4.2.1  Coating  via  Sol-Gel  Technique 

The  microstructures  of  electrophoretically  sol-gel  coated  Nb/MoSi2  interface 
with  different  processing  conditions  are  shown  in  Fig.  4.16.  As  seen  in  the  figure,  the 
coatings  have  degraded  during  the  processing  and  have  become  discontinuous,  thus 
allowing  for  direct  contact  and  interdiffusion  between  Nb  and  MoSi2  in  the 
discontinuous  regions.  The  degradation  of  the  coating  is  likely  partially  due  to  the 
presence  of  Si02  in  the  commercial  MoSi2  powder.  As  discussed  in  Section  4.1.1,  the 
as-received  MoSi2  powder  contains  SiC>2  and  Mo5Si3,  most  of  which  end  up  at  triple 
points  and  grain  boundaries  after  consolidation.  Experimental  results  in  Section  4.1.3 
suggested  that  AI2O3  may  react  with  SiC>2  present  in  the  commercial  MoSi2  powder. 
Therefore,  the  degradation  of  the  coating  at  the  present  case  is  likely  due  to  the  reaction 
between  the  AI2O3  coating  and  Si02  particles  at  the  region  where  the  two  are  in 
contact.  Another  reason  for  the  breakdown  of  the  continuous  coating  could  be  due  to 
cracking  during  drying  and  firing  of  the  gel.  It  was  found  that  when  the  coating  was 
thick,  as  in  the  cases  of  Fig.  4.16  (b)  and  (c),  cracking  occurred  more  readily  during 
drying  and  firing  in  the  furnace  or  even  during  the  air  drying  for  the  case  of  Fig.  4. 16 
(c).  The  addition  of  glycerol  prevented  cracking  somewhat  when  the  coating  was  thin, 
but  not  when  it  was  thick.  Thus,  the  discontinuity  of  the  coating  is  primarily  attributed 
to  the  attack  of  Si02  in  the  case  of  thin  coatings  and  additional  cracking  during  the 
drying  and  firing  in  the  case  of  thick  coatings. 

Notwithstanding  the  above  mentioned  localized  attack,  the  application  of  the 
coating  reduced  the  thickness  of  the  interaction  layer  between  Nb  and  MoSi2,  e.g., 
from  ~ 30  pm  for  uncoated  samples  shown  in  Fig.  4.5  to  ~ 15  pm  for  the  coated 
samples  shown  in  Fig.  4.16.  In  addition,  it  could  be  seen  from  Fig.  4.16  that  the 
thickness  of  the  oxide  coating  increased  with  decreasing  the  ratio  of  water  to  ASB.  This 
result  suggests  that  the  higher  the  concentration  of  ASB,  the  higher  the  rate  of  the 
build-up  of  hydrated  aluminum  hydroxide,  i.e.,  AIO(OH)  which  is  positively  surface- 
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Fig.  4.16  Microstructures  of  MoSi2/sol-gel  coated  Nb  interfaces  hot  pressed  at 

140()0C,  40  MPa  for  1 hour.  The  coatings  were  produced  via  electrophoresis 
with  the  ratios  of  H2O  : ASB  and  A1(N03)3  :ASB  in  the  sols  being  (a) 

100  :1  and  0.04  :1,  (b)  50  :1  and  0.04  :!  and  (c)  25  :1  and  0.04  :1, 
respectively. 
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charged  with  a formula  of  A10(0H2)  + , on  the  niobium  surface.  Furthermore,  it  was 
found  that  if  the  coating  was  thick  and  dense  enough,  the  diffusion  of  Si  across  the 
coating  could  be  minimized  An  example  is  provided  in  Fig.  4.17  which  was  taken 
from  a region  where  there  was  a relatively  long  continuous  alumina  coating.  The 
coating  in  the  figure  was  obtained  via  electrophoresis  using  SG5  sol  (Table  3.2)  with  a 
DC  voltage  of  1 volt.  As  seen  in  the  figure,  the  thickness  of  the  interaction  region  has 
been  drastically  reduced  to  ~ 1 pm  from  ~ 30  pm  for  the  uncoated  samples.  The 
reaction  product  of  the  electrophoretically  coated  samples  could  not  be  determined 
because  the  thickness  of  the  interaction  zone  was  close  to  the  beam  size  of  the  electron 
microprobe  (1  pm).  Nevertheless,  the  result  clearly  indicates  that  if  the  coating  is  dense 
enough  and  continuous,  the  diffusion  of  Si  can  be  effectively  limited. 

In  an  effort  to  further  establish  the  effect  of  the  coating  quality  on  the  elemental 
diffusion  and  interfacial  interaction,  the  following  experiments  were  conducted.  The 
experiments  involved  making  alumina  laminae  using  tape  casting  technique,  sintering 
the  tape,  placing  the  sintered  tape  in  between  the  MoSi2  and  Nb  laminae,  and  hot 
pressing  them  at  1400®C  for  1 hr.  The  pre-sintering  of  the  alumina  was  conducted  at 
1400^C  for  1 hour  with  a heating  rate  of  3®C/min.  Fig.  4.18  shows  the  result  of  the 
experiment.  As  seen  in  the  Figure,  the  Si  diffusion  has  been  totally  eliminated, 
suggesting  that  if  the  oxide  layer  is  thick  and  dense,  the  interdiffusion  and  interphase 
formation  can  be  prevented  thoroughly. 

The  above  discussion  suggests  that  as  the  thickness  of  the  coating  increases,  the 
ability  of  the  coating  to  retard  Si  diffusion  increases.  In  addition,  for  the  present 
composite  system  which  contains  SiC>2,  a thick  and  dense  coating  is  necessary  to 
defend  the  attack  of  SiC>2.  Thus,  it  seems  that  a low  water/ASB  ratio  to  produce  a thick 
coating  is  more  favorable.  However,  there  is  a lower  limit  to  the  water/ASB  ratio 
below  which  bulk  gelation  occurs  readily  and  A10(0H)  + hardens,  not  responding  to 
the  electric  field  [171].  Therefore,  there  is  an  optimum  ratio  which  allows  for  an 
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Fig.  4.17  Microstructure  and  composition  profile  of  a MoSi2/sol-gel  coated  Nb  interface, 
hot  pressed  at  1400°C,  40  MPa  for  1 hour.  The  coating  is  obtained  using  SG5 
sol  (Table  3.2)  with  a DC  voltage  of  1 volt. 


Fig.  4.18  Interfaces  of  MoSi2/pre-sintered  AI2O3  tape/Nb  compact,  hot  pressed  at 
1400^C,  40  MPa  for  1 hour. 
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adequate  thickness  of  the  coating  without  bulk  gelation.  For  the  present  study,  this  ratio 
was  found  to  be  25:1,  indicated  as  the  solution  SG5  in  Table  3.2.  However,  the 
cracking  problem  associated  with  the  thick  coating  still  presents  a formidable 
processing  challenge  which  is  yet  to  be  solved. 

No  obvious  effect  of  the  ratio  of  aluminum  nitrate  to  ASB  on  the  quality  and 
thickness  of  the  coating  was  observed  for  the  two  ratios  investigated  in  the  present 
study.  Finally,  the  AI2O3  coating  formed  using  dip-coating  technique  was  found  to  be 
thinner  and  more  porous  than  that  produced  with  electrophoresis  technique. 

4.2.2  Coating  via  Physical  Vapor  Deposition. 

The  microstructures  of  MoSi2  plates  diffusion-bonded  to  PVD  coated  Nb  with 
different  thicknesses  of  the  alumina  coating  are  shown  in  Fig.  4.19.  Interphase 
formation  was  detected  in  all  cases.  However,  the  thickness  of  the  interaction  zone  was 
reduced  dramatically  from  ~ 30  pm  for  uncoated  samples  (Fig.  4.5)  to  15,  5 and  1 pm 
for  PVD  coated  samples  with  a coating  of  1,  2 and  3 pm  thicknesses,  respectively.  The 
result  clearly  indicates  that  the  thickness  of  the  interaction  zone  decreases  with 
increasing  thickness  of  the  coating,  a trend  also  observed  with  sol-gel  coatings. 
Detailed  examination  of  Fig.  4.19  also  reveals  that  degradation  of  the  oxide  coating 
occurs  in  all  cases,  but  is  quite  localized,  as  can  be  seen  clearly  in  Fig.  4.19(c).  Since 
the  degradation  is  not  an  overall  problem  of  reaction  between  AI2O3  and  MoSi2,  but  a 
localized  attack,  it  is  again  most  likely  due  to  the  presence  of  Si02.  It  should  be  noted 
that  the  breakdown  of  the  coating  is  not  due  to  the  damage  by  MoSi2  particles  during 
the  processing.  This  can  be  illustrated  using  Fig.  4.20  which  compares  the  interfaces 
produced  with  different  processing  conditions.  Figure  4.20(a)  was  obtained  by  hot 
pressing  the  coated  Nb  with  MoSi2  powder  directly,  and  Fig.  4.20(b)  by  hot  pressing 
the  coated  Nb  with  the  pre-sintered  and  polished  MoSi2  plates.  It  is  clear  that  the 
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Fig.  4. 19  Microstructures  of  MoSi2/PVD  coated  Nb  interfaces,  hot  pressed  at 

1400°C,  40  MPa  for  1 hour,  (a)  1 (lm  AI2O3  coating;  (b)  2 p.m  AI2O3 
coating;  (c)  3 |im  AI2O3  coating. 
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Fig.  4.20  Effect  of  the  contact  of  the  thin  oxide  coating  with  MoSi2  powder  and  plates  on 
the  integrity  of  the  coating,  (a)  the  coating  in  contact  with  the  powder;  (b)  the 
coating  in  contact  with  the  pre -polished  plate.  (Hot  pressed  at  1400°C,  40  MPa 
for  1 hour.) 


Fig.  4.21  Microstructure  and  composition  profile  of  Al/Nb  interface  in  a hot  dipped  Nb 

sample.  The  temperature  of  the  molten  bath  is  930°C  and  dipping  time  is  2 min. 
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coating  characteristics  using  MoSi2  powder  are  similar  to  those  of  using  pre-sintered 
and  polished  MoSi2  plates. 

Comparing  the  results  of  PVD  coating  and  sol-gel  coating  indicates  that  the 
former  technique  produces  a much  denser  AI2O3  on  Nb  than  the  latter  one;  in  the  case 
of  sol-gel  coating  the  interphase  thickness  was  reduced  to  ~ 1 pm  only  when  the  coating 
thickness  reached  ~ 10  pm.  For  the  PVD  coating,  the  same  reaction  zone  thickness  was 
observed  when  the  coating  thickness  was  * 3 pm.  The  reduction  in  the  zone  thickness  is 
directly  related  to  the  density  of  the  coatings  produced.  For  PVD,  the  oxide  compound 
is  transformed  to  the  vapor  state  which  then  packs  on  the  surface  of  Nb.  In  contrast, 
sol-gel  coating  involves  a large  volume  shrinkage  associated  with  the  crystallization  of 
the  gel,  leading  to  the  formation  of  pinholes,  porosity  or  cracks  in  the  coating. 

4.2.3  Coating  via  Hot  Dipping  and  Anodizing. 

A typical  microstructure  and  compositional  profiles  of  Al/Nb  interface  in  a Nb 
sample  hot  dipped  in  a molten  aluminum  bath  at  930®C  for  2 minutes  are  shown  in  Fig. 
4.21.  As  can  be  seen,  the  hot  dipping  process  produced  a sharp  Al/Nb  interface  with  a 
thickness  of  ~ 1 pm.  The  compositional  analysis  also  shows  negligible  interdiffusion 
between  A1  and  Nb  during  the  dipping  operation.  In  addition,  the  interface  is  very 
clean  and  no  oxygen  is  detected  at  the  interface  by  the  compositional  analysis,  as  shown 
in  Fig.  4.21.  The  thickness  of  the  aluminum  coating  ranged  from  25  pm  to  40  pm. 
Hot  dipping  at  870®C  for  5 minutes  showed  a similar  result  except  that  the  uniformity 
of  the  coating  was  reduced. 

Figure  4.22  shows  the  microstructure  and  composition  profile  of  the  coating 
after  being  anodized  with  a current  density  of  ~ 14  mA/cm2  for  30  min.  As  seen  in  the 
figure,  aluminum  has  been  converted  to  alumina  and  the  thickness  of  the  coating  has 
been  reduced  from  ~ 30  pm  to  ~ 10  pm  during  the  anodizing  process.  Moreover,  the 
interface  remains  clean  and  sharp  after  the  anodizing.  It  was  found  that  the  anodizing 


i 


70 


Table  4.4  Thickness  of  AI2O3  layer  as  a function  of  the  current  density  during 

the  anodizing  of  pure  A1  plates. 


Current  density 

Anodizing  time 

Thickness 

(mA/cm^) 

(minute) 

(4m) 

5-7 

30 

- 2 

8-  10 

30 

-5 

12-  14 

30 

- 10 

23-25 

30 

-25 

time  and  current  density  had  a profound  influence  on  the  thickness  of  the  AI2O3 
formed,  as  given  in  Fig.  4.23  and  Table  4.4.  Fig.  4.23  shows  the  variation  with  time 
of  the  current  density  and  thickness  of  the  alumina  formed  for  anodizing  pure 
aluminum  plates  (99.999  wt%)  under  a constant  voltage.  As  expected,  the  thickness  of 
the  alumina  increases  with  time  rapidly  at  the  early  stage  of  anodizing,  but  it  slows 
down  after  a certain  period  of  time  (about  30  min),  suggesting  that  the  alumina 
approaches  a limit  thickness  asymptotically  for  a given  voltage.  Correspondingly,  the 
current  density  decreases  rapidly  at  the  early  stage  of  the  anodizing,  and  levels  off 
when  the  alumina  becomes  thick.  The  effect  of  current  density  on  the  thickness  of  the 
oxide  film  is  shown  in  Table  4.4.  As  seen  in  the  table,  increasing  the  current  density 
tends  to  increase  the  limit  thickness  of  the  oxide  film.  The  current  density  dependence 
of  the  limit  thickness  suggests  that  the  oxide  film  formed  dissolves  in  the  electrolyte 
during  the  anodizing  process  and  the  limit  thickness  is  a result  of  the  rate  of  dissolution 
matching  the  rate  of  formation.  These  results  are  in  agreement  with  the  literature 
[150,172,173]. 


71 


Dictance,  p m 


Fig.  4.22  Microstructure  and  composition  profile  of  AI2O3/NI)  interface  after  a hot  dipped 
Nb  has  been  anodized  with  a current  density  of  14  mA/cm^  for  30  min. 
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Fig.  4.23  Current  density  and  thickness  of  the  alumina  film  vs.  time  during  the  anodizing 
with  a constant  voltage. 
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When  anodizing  the  hot  dipped  Nb,  the  situation  was  complicated  by  the  fact 
that  the  total  thickness  of  the  coating,  i.e.,  the  (A1  4-  AI2O3)  layer,  on  the  surface  of 
Nb  became  thinner  and  thinner  as  anodizing  proceeded  because  of  the  dissolution  of  the 
oxide  in  the  electrolyte  [173].  If  the  anodizing  time  was  too  long,  the  entire  coating 
would  dissolve,  exposing  interior  Nb.  If  the  anodizing  time  was  too  short,  some  A1 
remained  unconverted  which  would  reacted  severely  with  Nb  during  the  subsequent 
consolidation.  Therefore,  when  anodizing  the  hot  dipped  Nb,  it  is  necessary  to  take  into 
account  the  decrease  of  the  thickness  of  (A1  + AI2O3)  layer  in  addition  to  the 
consideration  of  the  effect  of  the  current  density  on  the  thickness  of  the  AI2O3  film. 
The  present  study  showed  that  for  a current  density  of  ~ 12  mA/cm2  and  an  electrolyte 
with  5 wt%  sulfuric  acid,  the  rate  at  which  the  thickness  of  (A1  + AI2O3)  layer 
decreased  was  ~ 0.45  pm/min.  Thus,  knowing  the  thickness  of  A1  from  the  hot  dipping 
and  the  current  density  for  the  anodizing,  it  is  possible  to  calculate  the  anodizing  time 
to  control  the  coating  quality. 

Figure  4.24  shows  the  microstructure  and  corresponding  composition  profile  for 
a MoSi2/anodized  Nb  interface.  The  processing  condition  for  the  Nb  foil  was  hot 
dipping  in  the  molten  aluminum  at  930^C  for  2 min  and  anodizing  for  45  min  with  a 
current  density  of  ~ 12  mA/cm2.  As  seen  in  the  figure,  the  thickness  of  the  coating  is  ' 
10  pm  and  a reasonably  dense  coating  is  formed.  The  thickness  of  the  interaction 
region  has  been  reduced  from  ~ 30  pm  for  the  uncoated  samples  to  ~ 4 pm  for  the 
coated  samples.  Note  that  for  sol-gel  coating  technique,  when  the  coating  was  ~ 10  pm, 
the  interaction  zone  was  limited  to  ~ 1 pm,  indicating  that  the  coating  produced  via  sol- 
gel  processing  is  denser  than  that  produced  via  hot  dipping  and  anodizing  technique. 

The  effect  of  the  anodizing  current  density  on  the  quality  of  AI2O3  coating  is 
shown  in  Fig.  4.25.  The  difference  in  processing  condition  for  the  two  micrographs  is 
that  the  dipped  Nb  was  anodized  with  ~ 12  mA/cm2  for  45  min  for  Fig.  4.25(a),  and 
with  - 20  mA/cm2  for  30  min  for  Fig.  4.25(b).  As  seen  in  Fig.  4.25,  the  final 
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Fig.  4.24  Microstructure  and  composition  profile  of  MoSi2/anodized  Nb  interface. 

Processing  condition:  Nb  foil  is  hot  dipped  at  930^C  for  2 min,  anodized  with 
~12  mA/cm2  for  45  min,  and  finally  hot  pressed  with  MoSi2  at  1400°C, 

40  MPa  for  1 hour. 


Fig.  4.25  Effect  of  the  anodizing  current  density  on  the  quality  of  AI2O3  coating,  hot 

pressed  at  MOO^C,  40  MPa  for  1 hour,  (a)  Nb  foil  is  hot  dipped  at  930^0  for 
2 min,  anodized  with  ~12  mA/cm2  for  45  min;  (b)  Nb  foil  is  hot  dipped  at 
930^C  for  2 min,  anodized  with  ~20  mA/cm2  for  30  min. 
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thicknesses  of  the  coatings  are  very  similar  ( ~ 10  mm  ) for  both  processing  conditions. 
However,  the  interphases  formed  have  different  thicknesses  with  Fig.  4.25(a)  having 
about  4 mm  and  Fig.  4.25(b)  about  12  mm,  indicating  that  the  oxide  film  formed  with 
the  condition  of  Fig.  4.25(b)  is  more  porous  than  that  formed  with  the  condition  of 
Fig.  4.25(a).  The  results  imply  that  better  coatings  are  produced  using  lower  current 
densities.  However,  lower  current  densities  tend  to  produce  thinner  oxide  layers.  Thus, 
an  optimized  anodizing  current  density  is  needed  to  produce  a coating  with  proper 
density  and  thickness.  For  the  present  MoSi2  matrix  which  contains  SiC>2,  it  seems  that 
a medium  current  density  (12  - 14  mA/cm^)  which  produces  a coating  with  proper 
combination  of  thickness  and  density  would  be  recommended. 

In  conclusion,  the  present  set  of  experiments  has  demonstrated  that  a thick  and 
dense  AI2O3  coating  is  effective  in  eliminating  the  interactions  between  MoSi2  and  Nb 
totally,  and  the  best  coating  is  provided  with  physical  vapor  deposition,  followed  by 
sol-gel  technique. 


4.3  Fracture  Energy  of  the  Matrix/Reinforcement  Interfaces 
In  this  section,  the  fracture  energies  of  MoS^/Nb  interfaces  with  and  without 
coatings  will  be  discussed.  The  relationships  between  processing  routes, 
microstructures  and  mechanical  properties  of  the  interfaces  will  be  explored. 

4.3.1  MoSh/Uncoated  Nb  Interface 

Typical  interfacial  microstructures  of  MoSi2/Nb  composites  have  been  shown  in 
Figures  4.4  and  4.6.  When  the  hot-pressing  temperature  was  1400^C,  the  interfaces 
between  the  uncoated  Nb  and  MoSi2  was  MoSi2/(Nb,Mo)Si2/Nb5Si3/Nb  (Fig.  4.4), 
while  they  became  MoSi2/(Mo,Nb)5Si3/Nb  when  the  hot-pressing  temperature  was 
1700^C  (Fig.  4.6).  Despite  the  different  interfacial  microstructures,  mechanical 
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behaviors  were  observed  to  be  almost  the  same  for  the  samples  processed  at  the  two 
different  temperatures.  Figure  4.26  shows  a typical  load-displacement  curve  of  the 
chevron  notched  short  bar  for  a MoS^/uncoated  Nb  laminated  composite  hot  pressed  at 
1 400^0 . A stable  crack  propagation  has  been  achieved,  as  indicated  by  the  fluctuation 
in  the  curve  near  the  maximum  load.  Therefore,  the  maximum  load  was  used  to 
calculate  Gmax  with  the  aid  of  eq.  (3.1).  The  corresponding  fracture  surface  at  the 
chevron  crack  is  also  presented  in  Fig.  4.27.  By  examining  the  two  fracture  surfaces  of 
a broken  short  bar  with  scanning  electron  microscopy  (SEM)  and  energy  dispersive 
spectrometer  (EDS),  the  failure  location  at  the  interface  was  determined.  The  results 
are  summarized  in  Table  4.5.  For  comparison,  the  fracture  energy  of  monolithic  MoSi2 
measured  by  this  method  is  also  included  in  the  table.  It  was  noted  that  the  failure  of 
MoS^/uncoated  Nb  specimens  took  place  inside  MoSi2  rather  than  at  the  interfaces, 
indicating  that  the  interfacial  region  has  a higher  toughness  than  the  matrix.  Indeed,  the 
value  of  Gmax  of  MoSi2  measured  in  the  present  study  was  33.7  J/m^,  and  all  the 
measurements  conducted  on  MoS^/uncoated  Nb  specimens  showed  a fracture  energy 
of  about  36.4  J/m^  because  the  crack  in  MoSi2  failed  to  follow  the  chevron-notched 
plane  strictly. 

The  above  discussion  shows  clearly  that  the  failure  of  the  notched  short  bar 
specimens  does  not  necessarily  propagate  along  an  interface.  In  the  following  sections, 
more  examples  will  be  provided  to  show  that  depending  on  the  processing  conditions, 
the  crack  could  propagate  along  an  interface,  by  kinking  between  two  interfaces,  or 
only  inside  MoSi2  matrix.  Owing  to  this  feature,  the  measured  fracture  energy  of 
interface  using  the  notched  short  bar  has  been  referred  to  as  "fracture  energy  of  the 
interfacial  region"  of  the  laminates  in  the  present  study. 


4.3.2  MoSb/AbO^  Coated  Nb  Interface 

Typical  microstructures  of  MoS^/Nb  interface  coated  with  AI2O3  using 
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Fig.  4.26  A typical  load-displacement  curve  of  chevron  notched  short  bar  for  a 

MoSi2/uncoated  Nb  laminated  composite  and  the  corresponding  fracture 
surface. 
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Table  4.5.  Fracture  energy  and  failure  location  of  the  interfacial  region  in  the  coated 

and  uncoated  Nb/MoSi2  systems 


System 

Failure  location 

Toughness, 

Gmax, 

MoSi2  alone 

33.7  ± 1.4 

MoSi2/uncoated  Nb 

Inside  the  MoSi2 

> 33.7 

MoSi2/Al2C>3  coated  Nb  via 
sol-gel  technique 

Al203/(Nb,Mo)Si2  interface 
or  Al203/Nb5Si3  interface 

16.1  ± 1.3 

MoSi2/Al2C>3  coated  Nb  via 
hot  dipping  and  anodizing 

Al203/NbSi2  interface 
or  Al203/Nb5Si3  interface 

15.5  ± 1.6 

MoSi2/Al2C>3  coated  Nb  via 
PVD 

Partially  inside  the  MoSi2  and 
partially  at  Al2C>3/(Nb,Mo)Si2 
interface 

31.7  ± 3.4 

MoSi2/Zr02  coated  Nb 

ZrC>2/Nb5Si3  interface 
or  inside  the  ZrC>2 

12.8  ± 1.0 
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different  coating  techniques  have  been  presented  in  Fig.  4.16  (sol-gel  coating),  Fig. 
4.19  (PVD),  and  Fig.  4.25  (hot-dipping-and-anodizing).  For  convenience,  these 
microstructures  are  presented  in  Fig.  4.27  again.  Included  for  comparison  is  also  a 
micrograph  of  MoS^/ uncoated  Nb  interface. 

A typical  load-displacement  curve  of  the  hot-dipping-and-anodizing  coated 
Nb/MoSi2  laminate  and  the  corresponding  fracture  surfaces  at  the  chevron  crack  are 
shown  in  Fig.  4.28.  It  is  clear  that  the  failure  of  the  laminate  occurred  mainly  along  the 
Al203/NbSi2  interface  and  a small  portion  of  the  failure  was  in  the  Nb5Si3  phase.  The 
latter  was  probably  due  to  the  discontinuity  of  the  oxide  coating  generated  during  the 
anodizing  process  and/or  the  roughness  of  the  interface.  Depending  on  the  processing 
control  parameters,  the  interphase  formed  can  be  only  one  phase,  i.  e.,  NbsSi3  (Fig. 
4.24  and  Fig.  4.25(a)),  or  two  phases,  NbSi2  and  Nb5Si3  (Fig.  4.25(b)  and  Fig. 
4.27(d)).  However,  the  measured  fracture  energies  of  the  interfacial  region  for  the  two 
cases  are  very  similar.  These  results  have  been  summarized  in  Table  4.5. 

Fracture  surface  for  the  sol-gel  coated  Nb/MoSi2  laminate  are  shown  in  Fig. 
4.29.  It  is  noted  that  most  of  the  failure  occurred  along  Al203/(Nb,Mo)Si2  interface, 
and  a small  portion  of  the  failure  occurred  in  MoSi2,  a situation  similar  to  the  hot- 
dipping-and-anodizing  case.  Correspondingly,  the  specimens  using  the  two  coating 
techniques  mentioned  above  exhibited  a similar  fracture  energy  of  the  interfacial 
region,  as  given  in  Table  4.5.  Again,  the  measured  fracture  energies  of  the  interfacial 
region  for  the  case  of  only  one  interphase  formed,  i.e.,  Nb5Si3  (Fig.  4.16(b))  and  for 
the  case  of  two  interphases  formed,  (Nb,Mo)Si2  and  Nb5Si3  (Fig.  4.16(a)),  are  nearly 
the  same,  as  observed  with  the  hot-dipping-and-anodizing  coated  specimens. 

Quite  different  fracture  surface  was  observed  for  the  PVD  coated  specimens  (1 
pm  thickness  coating).  As  shown  in  Fig.  4.30,  failure  location  in  this  case  was 
partially  in  MoSi2  and  partially  at  the  interface  of  Al203/(Nb,Mo)Si2.  This  is  believed 
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Fig.  4.27  Interfacial  microstructures  of  MoSi2/Nb  laminated  composite  with  different 

processing  conditions,  (a)uncoated  Nb,  (b)  PVD  coated  Nb,  (c)  sol-gel  coated 
Nb,  and  (d)  hot-dipping-and-anodizing  coated  Nb. 
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Fig.  4.28  A typical  load-displacement  curve  of  chevron  notched  short  bar  for  a 

MoSi2/hot-dipping-and-anodizing  coated  Nb  laminated  composite  and  the 
corresponding  fracture  surfaces. 
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Fig.  4.29  General  view  of  the  two  corresponding  fracture  surfaces  of  a broken  short  bar  at 
the  chevron  crack  for  a MoSi2/sol-gel  coated  Nb  laminated  composite. 
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Fig.  4.31  An  enlarged  fracture  surface  of  the  tip  area  in  Fig.  4.30. 
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Table  4.6  Fracture  energy  of  the  interfaces 


System  (category) 

Gmax,  J/m^ 

MoSi2 

33.7 

Al203/(Nb,Mo)Si2 

(I) 

• 16 

Al203/NbSi2 

(I) 

• 16 

Al203/Nb5Si3 

(I) 

• 16 

Al203/MoSi2 

(I) 

>16? 

MoSi2/(Nb,Mo)Si2 

(II) 

> 33.7 

Nb5Si3/(Nb,  Mo)Si2 

(II) 

> 33.7 

Nb5Si3/Nb 

(III) 

> 33.7 

to  be  due  to  the  small  and  frequent  discontinuity  of  the  oxide  coating  (Fig.  4.27(b)) 
caused  by  the  presence  of  the  Si02  in  the  commercial  MoSi2  powder,  as  discussed  in 
Section  4.4.2.  It  appears  that  because  of  such  small  and  frequent  discontinuity  of  the 
coating,  the  fracture  resistance  along  the  Al203/(Nb,Mo)Si2  interface  has  been 
increased.  Owing  to  the  roughness  of  the  interface  (Fig.  4.27(b)),  crack  is  unable  to 
propagate  along  the  Al203/(Nb,Mo)Si2  interface,  but  have  to  kink  between  the 
interface  and  MoSi2  matrix  in  order  to  obtain  crack  propagation  with  minimum  energy 
consumption. 

It  is  apparent  that  various  interphases  were  formed  in  the  present  study  due  to 
the  diffusion  of  Si,  Mo  and/or  Nb  through  the  coating.  In  order  to  analyze  the  fracture 
energy  of  the  interface,  the  interfaces  at  the  interfacial  region  have  been  divided  into 
three  categories,  as  shown  in  Table  4.6.  Category  I consists  of  an  interface  (bonding) 
between  the  oxide  coating  and  intermetallics.  Category  II  is  composed  of  a bond 
between  two  intermetallics.  Category  III  is  the  interfaces  which  have  a bond  between  a 
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metal  and  an  intermetallic.  Because  of  the  formation  of  the  intermetallics,  no  direct 
bonding  between  oxide  and  metal  was  observed.  By  examining  the  failure  locations  and 
fracture  energy  of  the  interfacia]  region  measured,  the  fracture  energy  for  the  interfaces 
present  in  this  study  can  be  estimated  as  summarized  in  Table  4.6. 

It  is  noted  that  fracture  energy  for  the  oxide  bond,  category  I,  is  lower  than  the 
bonds  between  two  intermetallics  or  between  a metal  and  an  intermetallic  (categories  II 
and  III).  The  fracture  energies  of  the  former  is  at  a level  of  16  J/m^  in  the  present 
instance,  which  is  even  lower  than  the  fracture  energy  of  the  alumina  itself  (~  20  - 40 
J/m^).  The  low  fracture  energy  of  the  oxide  bond  is  attributed  to  the  low  bond  strength 
between  alumina  and  silicides  involved.  As  concluded  by  Sutton  and  Feingold  [174], 
the  bond  strength  between  oxide  and  other  materials  is  directly  related  to  the  amount  of 
interaction  between  the  two  materials,  and  the  free  energy  of  formation  of  the  other 
materials'  oxide  is  an  important  criterion  in  determining  the  interaction.  Weak  bonding 
will  form  when  the  conditions  are  unfavorable  for  the  formation  of  other  oxides.  In  the 
present  case,  elemental  Mo,  Nb  and  Si  have  a lower  free  energy  of  formation  of  their 
oxides  in  comparison  with  AI2O3  [175].  As  such,  the  formation  of  silica,  niobium  and 
molybdenum  oxides,  directly  or  by  decomposition  of  MoSi2,  is  unfavorable  in  terms  of 
free  energy  change.  Therefore,  it  is  expected  that  the  reactivity  of  the  silicides  with 
alumina  would  be  very  low,  resulting  in  a low  bond  strength.  On  the  other  hand,  if  the 
free  energy  of  formation  is  more  favorable  for  forming  other  oxides  instead  of  alumina, 
strong  bond  may  form.  An  example  is  given  by  the  work  of  Dalgleish  et  al.  [176]. 
They  found  that  crack  always  initiated  in  the  alumina  adjacent  to  the  interface  of  a Al- 
Mg  alloy  bonded  to  an  alumina.  High  bond  strength  in  this  case  can  be  attributed  to  a 
much  more  negative  free  energy  of  formation  for  MgO  than  that  of  AI2O3. 

Moreover,  the  fracture  energy  of  category  I interface  was  found  to  be 
insensitive  to  the  change  in  the  composition  of  the  silicides,  as  shown  in  Table  4.6. 
This  is  probably  due  to  the  similarity  of  Nb  and  Mo  in  their  atomic  structures.  If  this  is 
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true,  Al203/MoSi2  interface  should  show  similar  fracture  energy.  However,  no  failure 
at  Al203/MoSi2  interface  was  observed  in  the  present  study.  As  shown  in  Fig.  4.27(c) 
and  (d),  one  major  difference  of  Al203/MoSi2  interface  from  the  Al203/NbSi2  and 
Al203/(Nb,Mo)Si2  interfaces  is  that  the  former  is  rougher  than  the  latter.  Such 
roughness  of  the  interface  can  increase  fracture  resistance  by  the  mechanism  of 
interlocking  and  crack  deflection.  Thus,  Al203/MoSi2  interface  shows  a higher 
fracture  resistance  than  Al203/NbSi2  and  Al203/(Nb,Mo)Si2  interfaces,  although  their 
bond  strengths  may  be  similar. 

Category  II  and  III  interfaces  all  exhibit  higher  fracture  resistance  than  MoSi2 
matrix,  indicating  the  bond  strength  between  two  intermetallics  or  between  a metal  and 
an  intermetallic  is  high.  In  this  case,  crack  was  observed  to  propagate  inside  the  MoSi2 
matrix  instead  of  inside  the  interphases  formed.  However,  this  does  not  necessarily 
mean  that  MoSi2  has  a lower  inherent  fracture  resistance  than  that  of  the  interphases, 
(Nb,Mo)Si2  and  Nb5Si3.  It  is  believed  that  the  low  apparent  fracture  resistance  of 
MoSi2  in  the  present  instance  is  partially  related  to  the  high  porosity  and  SiC>2  in 
MoSi2  adjacent  to  the  interface,  as  shown  in  Fig.  4.27(a).  High  porosity  is  due  to  the 
fast  diffusion  of  Si  into  Nb  and  segregation  of  vacancy  onto  the  opening  channel  of 
MoSi2  powder  compact.  In  contrast,  the  interphases  formed  are  very  dense,  as  shown 
in  Fig.  4.27.  Because  the  pores  distribute  along  the  grain  boundary  of  MoSi2,  it  is 
expected  that  crack  would  propagate  along  the  grain  boundary  and  an  intergranular 
fracture  results.  Indeed,  this  is  confirmed  by  SEM  observation.  Figure  4.31  shows  an 
enlarged  fracture  surface  of  the  tip  area  in  Fig.  4.30.  Therefore,  low  apparent  fracture 
resistance  of  MoSi2  is  partially  caused  by  the  presence  of  porosity  and  Si02. 

4.3.3  MoSb/ZrOo  Coated  Nb  Interface 

The  interfacial  microstructure  of  MoSi2/Zr02  coated  Nb  has  been  presented  in 
Fig.  4.13  which  shows  that  the  interfacial  region  is  MoS^/ZrCtyNbsS^/Nb  for 
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MoSi2/ZrC>2  coated  Nb  systems.  A typical  load-displacement  curve  of  the  chevron- 
notched  short  bar  for  a MoS^/ZrCb  coated  Nb  laminate  and  the  corresponding  fracture 
surface  are  shown  in  Fig.  4.32.  Again,  the  fluctuation  in  the  curve  near  the  maximum 
load  is  taken  as  an  indication  of  the  stable  propagation  of  the  crack.  The  failure  location 
was  identified  as  at  the  Zr02/Nb5Si3  interface  or  inside  the  ZrC>2,  and  the  fracture 
energy  of  the  interface  was  measured  to  be  12.8  J/m^.  These  results  have  been  included 
in  Table  4.5  for  comparison.  It  is  noted  that  ZrC>2  coated  Nb  system  has  the  lowest 
interfacial  fracture  energy  among  the  systems  studied.  This  is  likely  due  to  the 
existence  of  residual  tensile  stresses  in  the  zirconia  coating  caused  by  a higher 
coefficient  of  thermal  expansion  of  the  zirconia  than  that  of  the  matrix  and 
reinforcement. 

In  conclusion,  the  experiments  in  this  section  have  demonstrated  that  fracture 
energy  of  an  interfacial  region  depends  on  the  interfacial  bond  strength,  roughness  of 
interface  and  microstructure  of  the  two  component  materials  at  the  interface.  Weak 
bond  strength  leads  to  a low  interfacial  fracture  energy.  Roughness  of  the  interface 
increases  the  interfacial  fracture  energy.  Porosity  in  one  of  the  two  component 
materials  at  the  interface  could  lead  to  a low  fracture  energy  of  the  interfacial  region 
due  to  the  crack  propagation  along  the  weak  path  in  one  of  the  component  materials 
instead  of  along  the  interface. 

4.4  Interactions  between  Cracks  and  Ductile  Reinforcements 

In  this  section,  the  crack  propagation  and  crack-ductile  reinforcement 
interactions  at  various  stages  of  the  loading  will  be  discussed.  The  technique  used  was 
by  unloading  the  chevron  notched  samples  at  various  loads  during  four-point  bending 
tests,  followed  by  sectioning  the  unloaded  samples  and  investigating  the  crack 
geometry.  The  tensile  test  specimens  shown  in  Fig.  3.2  were  also  used  to  aid  the 
investigation  of  the  characteristics  of  the  crack  interaction  with  ductile  reinforcements. 
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Fig.  4.32  A typical  load-displacement  curve  of  chevron  notched  short  bar  for  a 

MoSi2/Zr02  coated  Nb  laminated  composite  and  the  corresponding  fracture 
surface. 
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4.4.1  Relation  between  the  External  Load  and  Crack  Propagation  in 

MoSb/Nb  Composites 

Typical  load-displacement  curves  of  chevron  notched  specimens  for  a 
monolithic  MoSi2  and  an  uncoated-Nb-foil-reinforced  composite  are  shown  in  Fig. 
4.33.  The  load  levels  at  which  the  specimens  were  unloaded  are  indicated  by  points  (a) 
through  (e)  in  the  figure.  The  unloaded  specimens  were  sectioned  parallel  to  the  crack 
propagation  direction  but  perpendicular  to  the  direction  of  the  foil  alignment,  and  then 
examined  for  the  extent  of  crack  propagation.  Figures  4.34(a)-(e)  show  the  nature  of 
crack  propagation  in  the  samples  at  the  corresponding  load  levels  indicated  in  Fig. 
4.33.  For  the  sample  loaded  up  to  15  lbs  (point  (a)  in  Fig.  4.33)  which  is 
approximately  10  % of  the  peak  load  of  the  composite,  a crack  was  observed  at  the  tip 
of  the  chevron  notch,  which  had  propagated  in  a stable  manner  as  shown  in  Fig. 
4.34(a).  When  the  load  was  increased  to  35  lbs  (point  (b)  in  Fig.  4.33)  which  is  about 
23%  of  the  peak  load,  the  result  indicated  that  the  crack  had  already  propagated 
throughout  the  entire  thickness  of  the  matrix.  This  load  is  approximately  80  % of  the 
peak  load  of  the  monolithic  matrix.  Beyond  point  (b),  the  matrix  is  split  in  two 
sections,  and  many  more  secondary  cracks  were  also  observed. 

The  extensive  cracking  of  the  matrix  at  a load  level  considerably  below  the  peak 
load  of  the  composite  may  be  explained  using  the  principle  of  strain  compatibility  in 
composite  materials.  During  the  initial  loading  and  prior  to  cracking  of  the  matrix, 
most  of  the  applied  load  is  carried  by  the  matrix  because  elastic  modulus  of  MoSi2 
matrix  (380  GPa)  is  much  higher  than  that  of  Nb  foils  (105  GPa).  For  the  Nb  foils  to 
carry  the  same  stress  as  the  matrix  does,  it  needs  to  be  elongated  three  times  more  than 
the  matrix  (when  there  is  no  constraint  for  deformation  from  the  matrix),  which  is  not 
possible  if  the  matrix  is  not  cracked.  Therefore,  the  reinforcement  with  a lower  elastic 
modulus  will  carry  a minor  part  of  the  load  until  a matrix  crack  has  opened  appreciably 
to  allow  for  straining  of  the  reinforcement.  This  has  been  schematically  shown  in  Fig. 
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Fig.  4.33  Typical  load-displacement  curves  of  chevron-notched  monolithic  MoSi2  matrix  and  uncoated  niobium-foil- 
reinforced  composite  specimens.  The  letters  on  the  composite  curve  indicate  the  approximate  unloading  load 
levels  for  the  microstructures  shown  in  Fig.  4.34. 
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Fig.  4.34  Caption  overleaf. 
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Fig.  4.34  Caption  overleaf. 
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Fig.  4.34  The  propagation  of  the  cracks  in  uncoated  niobium-foil-reinforced  MoSi2  composites  at  various  load  levels 

shown  in  Fig.  4.33. 
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2.5,  where  the  reinforcements  which  are  farther  behind  the  crack  tip  are  strained  much 
more  than  the  ones  closer  to  the  tip  or  ahead  of  the  crack.  Obviously,  the  required 
level  of  the  crack  opening  for  the  additional  straining  of  the  reinforcement  will  depend 
on  the  length  of  debonded  section  of  the  reinforcement;  the  larger  the  debonding 
length,  the  larger  is  the  reinforcement  section  which  is  being  strained,  and  the  larger  is 
the  required  crack  opening.  Moreover,  the  latter  will  also  depend  on  the  relative  values 
of  the  elastic  moduli  of  the  matrix  and  the  reinforcement. 

Closer  examination  of  Fig.  4.34  reveals  interesting  aspects  of  the  interaction 
between  cracks  and  Nb  foils  in  the  composites  reinforced  by  the  uncoated  Nb  foils.  At 
the  early  stages  of  loading,  there  is  only  one  primary  crack  propagating  in  the  matrix. 
This  primary  crack  propagates  through  the  matrix  in  a step-wise  manner,  instead  of 
propagation  on  one  crack  plane  (Fig.  4.34(b)).  The  step-wise  crack  propagation  in  the 
matrix  may  be  explained  through  geometrical  effects  of  the  Nb  reinforcement  foils.  In 
these  laminated  composites,  as  a crack  ends  on  one  side  of  Nb  foils,  it  needs  to 
reinitiate  in  the  other  side,  making  the  crack  propagation  discontinuous  in  nature.  The 
reinitiation  plane,  however,  may  not  necessary  be  on  the  same  plane  of  the  crack  on  the 
other  side  of  the  foil  because  of  the  existence  and  statistical  distribution  of  defects  in 
MoSi2,  which  can  initiate  the  crack.  The  net  result  is  the  step-wise  crack  propagation. 

When  the  geometry  of  the  Nb  reinforcements  is  changed  from  foils  to  filaments, 
the  crack  propagation  in  the  matrix  is  also  expected  to  change.  For  filament  reinforced 
composites,  the  crack  can  propagate  continuously  throughout  the  matrix  without  the 
need  for  reinitiation.  Therefore,  the  crack  will  propagate  on  the  plane  of  the  maximum 
stress  intensity  and  the  secondary  cracks  are  not  expected  to  form.  Fig.  4.35  shows  the 
fractured  surface  of  a uncoated  Nb  filament  reinforced  composite.  As  seen  in  the 
figure,  there  is  no  appreciable  debonding  at  the  reinforcement/matrix  interface  and  the 
filaments  has  failed  in  a brittle  cleavage  manner.  The  brittle  fracture  of  the  filaments  is 
thought  to  occur  because  of  the  lack  of  decohesion  at  the  reinforcement/ matrix  interface 
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due  to  strong  bonding  and  because  of  the  absence  of  the  virtual  decohesion  by 
secondary  cracks  as  in  the  laminated  composites.  As  a result,  the  filaments  are 
constrained  by  the  matrix,  creating  a tri-axial  stress  state,  which  leads  to  the  brittle 
fracture  of  the  filaments. 

When  Nb  foils  were  coated  with  AI2O3,  the  cracks  in  the  laminated  composites 
propagated  in  a manner  similar  to  the  uncoated  Nb  foil  reinforced  composites:  step- 
wise propagation  of  the  primary  crack,  extensive  secondary  cracking  and  ductile  failure 
of  Nb  foils.  However,  when  Nb  foils  were  coated  with  ZrC>2,  the  mode  of  crack 
propagation  was  found  to  be  dependent  of  the  processing  temperature.  If  the  hot 
pressing  temperature  was  1400^C,  the  crack  was  observed  to  propagate  in  a manner 
similar  to  the  uncoated  Nb  foil  reinforced  composites.  However,  when  the  hot  pressing 
temperature  was  raised  to  1700^C,  a different  propagation  mode  was  observed,  as 
shown  in  Fig.  4.36.  In  this  composite,  there  was  no  extensive  secondary  cracking  and 
Nb  foils  failed  in  a brittle  manner,  even  though  the  primary  crack  propagated  in  the 
step-wise  manner.  The  brittle  fracture  of  Nb  foil  is  believed  to  be  due  to  the 
embrittlement  of  Nb  foils,  as  will  be  discussed  further  in  the  following  section. 

4.4.2  Behavior  of  Cracks  at  the  Interface 

The  scenario  of  how  cracks  behave  near  the  interface  will  be  explored  in  this 
section.  Details  of  cracks  impinging  on  Nb  foils  for  uncoated  and  coated  composites  are 
shown  in  Fig.  4.37.  The  micrographs  were  taken  from  the  cross  sections  of  laminated 
composites  unloaded  at  about  20%  of  the  peak  load  of  the  composites  during  the  four- 
point  bending  tests.  It  is  noted  from  the  figures  that  debonding  at  the  interface  does  not 
occur  when  a crack  approaches  the  interface  or  just  impinges  on  it  for  both  uncoated 
and  coated  reinforcement  systems.  This  is  not  surprising  since  all  the  fracture  energies 
of  the  interfaces  measured  in  the  present  study  are  higher  than  1/5  of  the  matrix 
fracture  energy.  The  observation  is  consistent  with  an  estimation  made  by  Cook  and 
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Fig.  4.35  Fracture  surface  of  a uncoated  niobium-filament-reinforced  MoSi2  composite, 
showing  limited  decohesion  at  the  interface  and  brittle  cleavage  fracture  of  the 
niobium  filament:  (a)  a fractured  chevron-notched  specimen;  (b)  a fractured 
niobium  filament. 
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failed  in  a brittle  manner. 
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Gordon  [177].  They  found  that,  for  an  elliptical  crack,  an  interfacial  fracture  energy  of 
1/5  or  less  of  the  matrix  fracture  energy  would  cause  interfacial  debonding  in  advance 
of  the  crack. 

Figure  4.38  shows  a typical  microstructure  of  an  uncoated  Nb  tensile  test 
specimen.  A notch  tip  can  be  seen  at  the  left  side  of  Fig.  4.38(a),  which  shows  that  a 
crack  initiates  at  the  notch  tip  and  ends  at  the  Nb  foil.  Fig.  4.38(b)  is  a close-up  of  one 
of  the  crack  tips  in  Fig.  4.38(a).  As  seen  in  the  figure,  the  impingement  of  the  crack  on 
the  interface  causes  local  dislocation  slip  of  the  reinforcement  instead  of  interfacial 
failure,  leading  to  the  release  of  the  stress  concentration.  In  addition,  there  are  also 
multiple  cracks  at  the  interface  which  add  to  relaxation  of  tip  stress  intensity  and 
effectively  enlarge  the  initial  plastic  zone  in  the  ductile  phase.  The  result  shows  that 
interfacial  failure  is  not  the  only  mechanism  of  blunting  cracks  in  the  case  of  ductile 
reinforcements.  Instead,  the  slip  capability  of  ductile  reinforcement  can  play  an 
additional  and  important  role. 

As  load  continues  to  increase  during  the  tensile  tests  after  the  cracks  have 
impinged  on  the  Nb  foil,  debonding  at  the  interfaces  and/or  multiple  fracture  of  the 
matrix  near  the  interfaces  occur  due  to  a relatively  large  lateral  deformation  of  the  Nb 
reinforcement  compared  to  the  matrix  and  load  transfer  from  the  matrix  to  the 
reinforcement.  Typical  features  of  debonding  at  the  interface  and  multiple  fracture  of 
the  matrix  near  the  interface  for  the  coated  and  uncoated  specimens  are  shown  in  Fig. 
4.39.  It  is  noted  that  interfacial  debonding  has  occurred  at  the  oxide  coated  systems 
(Fig.  4.39(a)),  as  contrasted  with  multiple  matrix  fracture  near  the  interface  for  the 
uncoated  systems  (Fig.  4.39(b)).  This  is  in  agreement  with  the  interface  fracture  energy 
measurement  since  the  interfaces  in  the  uncoated  composites  have  a higher  fracture 
energy  than  the  matrix,  as  such  the  matrix  is  expected  to  fail  more  easily  than  the 
interfaces  in  these  composites.  On  the  other  hand,  for  the  oxide  coated  composites  the 
fracture  energy  of  the  interfaces  is  lower  than  the  fracture  energy  of  the  matrix.  Thus, 


98 


Fig.  4.37  Optical  microstructures  of  composite  laminates  showing  that  cracks  end  at  the  front  of  the  reinforcement 
without  causing  interfacial  debonding,  (a)  Uncoated  Nb,  (b)  AI2O3  coated  Nb,  and  (c)  Z1O2  coated  Nb 
reinforced  composites. 
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Fig.  4.38  An  edge  view  of  a unloaded  tensile  test  specimen,  (a)  A general  view  of  the 
cracks  and  (b)  a close-up  of  the  crack  tips,  showing  dislocation  slip  of  the  Nb 
at  the  crack  tips. 


Fig.  4.39.  An  edge  view  of  tensile  tested  specimens,  (a)  Debonding  at  the  interface  in  a 
specimen  with  AI2O3  coated  Nb.  (b)  Multiple  matrix  fracture  near  the  interface 
in  a specimen  with  uncoated  Nb. 
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interfacial  debonding  prevailed  in  these  composites.  Similar  results  were  also  observed 
on  bending  tests  of  chevron  notched  composite  laminates. 

Both  interfacial  debonding  and  multiple  matrix  fracture  create  a "gauge  length" 
at  the  matrix/reinforcement  interface  which  is  a region  virtually  free  from  constraints  of 
the  matrix.  This  region  is  called  decohesion  length  in  this  paper.  Measurement  of  the 
decohesion  length  was  conducted  for  laminated  composites  and  the  results  are  shown  in 
Table  4.7.  Since  the  decohesion  length  of  the  laminated  composites  varies  with  load 
and  position  of  the  reinforcement,  the  values  reported  in  Table  4.7  were  measured  at 
the  peak  load  and  for  the  second  foil  from  the  notch  tip.  Microhardness  of  the  Nb 
reinforcements  after  hot  pressing  is  also  included  in  Table  4.7  to  show  how  hardness  of 
the  Nb  foils  was  affected  by  the  processing  conditions.  Since  little  change  in 
microhardness  across  the  Nb  foil  cross  section  was  observed,  the  microhardness  was 
taken  as  a constant  for  each  specific  condition.  The  results  for  the  specimens  hot 
pressed  at  1400°C  show  that  the  lower  the  fracture  energy  of  the  interface,  the  longer 
the  decohesion  length.  This  is  expected  since  lower  interfacial  fracture  energy  means 
lower  resistance  to  interfacial  decohesion.  However,  in  the  case  of  Zr02  coated 
systems,  when  hot  pressing  temperature  is  increased  to  1700°C,  the  decohesion  length 
becomes  negligible  in  contrast  with  the  general  trend  exhibited  by  the  composites 
processed  at  lower  temperatures.  Fig.  4.40  shows  two  typical  load-displacement  curves 
of  bending  tests  on  Zr02  coated  Nb/MoSi2  laminated  composites  hot  pressed  at 
1400°C  and  1700°C.  It  is  noted  that  fracture  energy  of  the  composites  is  reduced  and 
the  failure  becomes  catastrophic  when  hot  pressing  temperature  changes  from  1400°C 
to  1700°C.  Examination  of  fracture  surfaces  also  shows  a change  of  fracture  modes 
from  quasi-cleavage  to  cleavage,  as  shown  in  Fig.  4.41.  As  indicated  in  Table  4.7, 
microhardness  of  the  Nb  foils  increases  from  131  to  236  as  the  processing  temperature 
is  increased,  indicating  an  increase  in  the  slip  resistance  of  the  Nb  foils.  The 
embrittlement  of  Nb  is  probably  due  to  the  diffusion  of  interstitial  oxygen  into  the  foils 


Table  4.7  Decohesion  length  and  microhardness  measured  from  the  laminated  composites  reinforced  with 
20  vol.%  of  Nb  foils  with  a thickness  of  0.25  mm* 
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at  higher  processing  temperature  from  the  decomposition  of  Z1O2.  Loss  of  the  ductility 
of  the  foil  leads  to  little  lateral  deformation  of  the  Nb  and  resultant  transverse  stresses, 
which  results  in  negligible  decohesion  length.  To  support  this  inference,  tensile  tests  on 
single  constrained  reinforcement  specimens  were  conducted  and  the  results  showed  that 
whenever  Vickers  hardness  of  uncoated  Nb  foils  increased  to  about  200  by  controlling 
hot  pressing  conditions,  the  fracture  mode  was  always  by  cleavage  even  if  the  notches 
were  cut  directly  into  the  Nb  foils.  This  result  shows  clearly  that  brittle  fracture  of  the 
Nb  foils  is  due  to  embrittlement  of  the  foils. 

From  the  above  discussion,  it  can  be  summarized  that  as  a crack  approaches  the 
ductile  reinforcements  in  these  systems,  there  exist  three  competitive  mechanisms: 
interfacial  debonding,  multiple  matrix  fracture  and  direct  crack  propagation  through  the 
ductile  reinforcement.  Figure  4.42  shows  schematically  these  interactions.  For  the  case 
of  high  ductility  reinforcements,  the  crack  is  blunted  by  local  dislocation  slip  of  the 
ductile  phase.  The  present  experiments  showed  that  extensive  debonding  did  not  occur 
at  the  early  stage  of  the  crack/ductile  phase  interaction,  so  that  the  interaction  could  be 
approximated  as  the  case  of  no  debonding  at  this  stage.  Thus,  to  fail  the  ductile 
reinforcement,  a much  higher  tensile  stress  has  to  be  applied  which  increases  shear 
stresses  at  the  interface  due  to  load  transfer  and  transverse  stresses  caused  by  the 
difference  between  lateral  displacements  of  the  matrix  and  reinforcement.  Both  shear 
and  transverse  stresses  enhance  interfacial  debonding  and/or  multiple  matrix  fracture. 
Thus,  whether  interfacial  debonding  or  multiple  matrix  fracture  predominates  depends 
on  the  values  of  the  fracture  energies  of  the  matrix  and  interface.  On  the  other  hand, 
for  the  case  of  low  ductility  reinforcements,  the  crack  can  relatively  easily  propagate 
through  the  reinforcement  before  occurrence  of  the  interfacial  debonding  or  multiple 
matrix  fracture.  The  decohesion  length,  which  results  from  the  competition  of  the  three 
above  mentioned  mechanisms,  is  controlled  by  a combination  of  three  material 
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Fig.  4.40  Typical  load-displacement  curves  of  bending  tests  on  chevron  notched 
composite  laminates  reinforced  by  Z1O2  coated  Nb  foils. 


Fig.  4.41  Fracture  surfaces  of  the  Nb  foils  in  composite  laminates  reinforced  by  Z1O2 

coated  Nb,  (a)  hot  pressed  at  1400°C  for  1 hour  and  (b)  hot  pressed  at  1700°C 
for  40  min. 
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(a)  Interfacial  debonding  (b)  Multiple  matrix  fracture 


(c)  Direct  crack  propagation 


Fig.  4.42  Schematic  of  the  three  competitive  mechanisms  at  the  interface.  Conditions 

for  the  process  to  dominate  shown  in  (a)  are  high  toughness  reinforcement  and 
weak  interface,  (b)  are  high  toughness  reinforcement  and  strong  interface,  (c) 
are  low  toughness  reinforcement  and  strong  or  weak  interface. 
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properties:  interface  fracture  energy,  toughness  of  the  matrix  and  slip  capability  of  the 
reinforcement. 

4.4.3  Fracture  Surface  Observation  on  MoSn  Matrix 

An  overall  view  of  the  fracture  surface  of  chevron  notched  specimens  for  the 
monolithic  MoSi2  is  presented  in  Fig.  4.43.  Higher-magnification  views  of  the 
fracture  surfaces  of  the  monolithic  MoSi2  at  different  positions  relative  to  the  notch  are 
shown  in  Fig.  4.44  where  figures  (a),  (b),  and  (c)  are  taken  from  the  positions  a,  b, 
and  c shown  in  Fig.  4.43,  respectively.  As  seen  in  the  figures,  intergranular  fracture 
prevails  at  the  tip  of  the  notch,  while  transgranular  fracture  dominates  at  the  bottom  of 
the  notch.  Between  these  two  positions,  a transition  of  fracture  mode  is  observed.  Most 
of  the  visible  second  phases  on  the  intergranular  fracture  surface  (Fig.  4.44(a))  have 
been  identified  as  amorphous  silica,  and  the  rest  as  carbon-stabilized  Mo5Si3  particles 
[158].  There  also  exist  some  pores  along  the  grain  boundary,  as  seen  on  the 
intergranular  fracture  surface.  The  presence  of  the  silica  glass  and  pores  weakens  the 
grain  boundary.  The  change  in  fracture  mode  with  the  position  is  ascribed  to  the 
different  natures  of  crack  propagation  at  different  positions.  Owing  to  the  geometry  of 
the  V-notch,  propagation  of  the  crack  is  stable  near  the  notch  tip,  unstable  at  the 
bottom  of  the  notch  and  changes  from  stable  to  unstable  at  the  intermediate  position 
[178,179].  Because  of  these  characteristics,  crack  propagation  near  the  notch  tip  is 
controlled  by  the  crosshead  speed  of  the  testing  machine,  and  it  becomes  fast  at  the 
bottom  of  the  notch.  Thus,  at  the  notch  tip  the  crack  has  more  time  to  choose  the 
energy  path  which  is  grain  boundary  for  the  present  matrix  due  to  the  presence  of  silica 
glass  and  pores,  while  at  the  bottom  of  the  notch  the  crack  propagates  too  fast  to 
choose  the  low  energy  path,  leading  to  transgranular  fracture. 

A general  view  of  a typical  fracture  surface  at  the  notch  plane  of  a laminate 
containing  20  vol.  % of  Nb  laminae  with  a thickness  of  0.25  mm  is  presented  in  Fig. 
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Fig.  4.43  Fracture  surface  of  a MoSi2  chevron  notched  specimen. 


Fig.  4.44  High  magnification  views  of  the  fracture  surfaces  of  MoSi2.  Figures  (a),  (b)  and 
(c)  are  taken  from  the  positions  a,  b and  c shown  in  Fig.  4.43,  respectively. 
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4.45.  The  a,  b,  c,  and  d in  Fig.  4.45  denote  the  positions  where  the  fractographs  of 
Fig.  4.46(a),  (b),  (c)  and  (d)  are  taken,  respectively.  Note  that  intergranular  fracture 
again  dominates  at  the  tip  of  the  notch,  while  transgranular  fracture  prevails  at  the 
bottom  of  the  notch,  a general  trend  also  observed  in  the  monolithic  MoSi2.  The 
presence  of  ductile  laminae  makes  crack  propagation  in  the  matrix  become 
discontinuous,  i.  e.,  the  crack  has  to  stop  at  one  side  of  the  Nb  lamina  and  reinitiates 
on  the  other  side  to  continue  the  propagation.  However,  the  position  dependence  of 
fracture  mode  of  MoSi2  in  the  composite  laminate  shows  little  difference  from  that  of 
the  monolithic  MoSi2.  Such  phenomena  could  be  explained  in  terms  of  the  discussion 
made  in  Section  4.4.1  where  it  was  shown  that  a primary  crack  has  already 
propagated  through  the  entire  MoSi2  matrix  at  a load  level  of  about  20%  of  the  peak 
load  in  the  bend  test  of  chevron  notched  specimens.  Therefore,  the  presence  of  Nb 
laminae  does  not  have  much  influence  on  the  fracture  mode  of  MoSi2  which  fractures 
at  the  early  stage  of  the  bend  test.  The  presence  of  Nb  laminae,  however,  improves 
the  damage  tolerance  of  the  brittle  materials,  as  will  be  discussed  in  Sections  4.5  and 
4.6. 


4.5  Flow  Behavior  of  Constrained  Niobium 
As  discussed  in  Section  2.3,  the  primary  toughening  mechanism  of  ductile 
reinforcement  is  due  to  the  bridging  of  ductile  ligaments  [20-26,105-107],  The  ductile 
bridging  ligaments  tend  to  hold  the  crack  faces  together  and  reduce  the  stress  intensity 
seen  at  the  crack  tip.  Thus,  the  key  to  predict  the  increased  fracture  toughness  is  to 
calculate  the  clamping  forces  of  the  bridging  ligaments,  ct(u),  as  a function  of  crack 
opening.  In  this  section,  the  clamping  forces  of  the  bridging  ligaments  as  a function 
of  crack  opening,  i.e.,  the  flow  behavior  of  the  constrained  ductile  phases,  will  be 
presented.  The  experimental  procedure  used  has  been  presented  in  Section  3.4. 
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Fig.  4.45  An  overall  view  of  fracture  surface  of  a laminated  composite  containing 
20  vol.%  of  Nb  laminae  with  a thickness  of  0.25  mm.  The  a,  b,  c and  d 
denote  the  positions  where  the  fractographs  of  Fig.  4.46(a),  (b),  (c)  and  (d) 
are  taken,  respectively. 


109 


Fig.  4.46  Fracture  surfaces  of  MoSi2  matrix  in  a laminated  composite.  Figures  (a),  (b),  (c)  and  (d)  are  taken  from  the  positions 

a,  b,  c and  d shown  in  Fig.  4.45,  respectively. 
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4.5.1  Intrinsic  Properties  of  Unconstrained  Niobium 

Mechanical  properties  measured  from  the  as-received  Nb  have  been 
summarized  in  Table  3.1.  However,  these  properties  cannot  represent  the  intrinsic 
properties  of  niobium  reinforcements  embedded  in  the  MoSi2  matrix  (called  as  the 
processed  Nb  in  the  text)  because  the  embedded  niobium  have  undergone  the  hot 
compaction  process.  Thus,  the  properties  of  the  processed  Nb  have  to  be  obtained  if  a 
more  accurate  relation  between  the  flow  behaviors  of  the  constrained  and 
unconstrained  ductile  phases  is  to  be  obtained. 

To  obtain  the  processed  Nb,  the  following  procedures  has  been  taken:  first, 
one  layer  of  ZrC>2  coated  Nb  foil  was  hot  pressed  with  MoSi2',  second,  the  hot 
pressed  discs  were  cut  into  rectangular  bars;  and  finally,  the  MoSi2  on  both  sides  of 
the  Nb  foil  was  broken  off  by  mechanical  means.  The  last  step  was  not  very  difficult 
due  to  a relative  weak  bonding  between  Nb  and  MoSi2  at  the  presence  of  a ZrC>2 
coating  layer. 

The  properties  measured  from  the  processed  Nb  are  summarized  in  Table  4.8. 
Comparison  between  Tables  3.1  and  4.8  reveals  that  strength  of  Nb  has  increased 
after  hot  pressing  except  for  the  1.0  mm  Nb  foils.  The  rate  of  strain  hardening  has 
also  increased  for  all  the  niobium.  Such  changes  are  believed  to  be  due  to  a 
combination  of  two  competitive  processes  during  hot  pressing:  (1)  grain  growth  of 
Nb,  as  shown  in  Fig.  4.47;  (2)  solid  solution  strengthening  because  of  its  contact  with 
the  silicides  formed  at  the  matrix/reinforcement  interface,  as  shown  in  Fig.  4.13.  For 
1.0  mm  Nb  foils,  the  grain  size  has  changed  from  15  to  ~500  pm,  which  leads  to  a 
large  decrease  in  strength  and  overshadows  the  increase  due  to  the  solid  solution 
hardening.  For  0.5  and  0.25  mm  Nb  foils,  the  grain  size  has  changed  from  50  to 
~370  pm  and  from  48  to  '300  pm,  respectively.  However,  shorter  diffusion  distance 
to  the  center  of  the  niobium  has  made  solid  solution  hardening  dominate  over  the 


Table  4.8.  Mechanical  properties  of  the  processed  Nb 
(Zr02  coated,  hot  pressed  at  1400°C,  40  MPa  for  1 hour) 
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Fig.  4.47  Optical  microstructures  of  Nb  foils,  (a)  as-received,  1.0  mm  thick;  (b)  processed,  1.0  mm  thick;  (c)  as-received, 
0.5  mm  thick  and  (d)  processed,  0.5  mm  thick.  Note  that  different  magnifications  are  used. 
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softening  due  to  the  grain  growth.  Therefore,  0.5  and  0.25  mm  Nb  foils  show  an 
increased  strength  after  hot  pressing. 

Examination  of  the  true  stress-strain  data  of  the  processed  Nb  indicates  that  as 
the  true  strain  exceeds  a certain  value,  the  true  stress-strain  follows  a power  law  quite 
well.  For  example,  the  true  strain  at  which  the  true  stress-strain  starts  to  follow  a 
power  law  is  * 0.34%  for  0.25  mm  processed  Nb,  and  is  ~ 1.0%  for  1.0  mm 
processed  Nb.  Thus,  the  true  stress-strain  data  after  the  above  mentioned  strains 
(e.g.,  0.34%  for  0.25  mm  Nb)  have  been  used  to  calculate  both  the  strength 
coefficient  and  the  strength  hardening  coefficient.  The  results  are  included  in  Table 
4.8  and  also  presented  below  for  convenience. 


ot  = 356  et  °-165  (MPa) 

for  1.0  mm  Nb 

ot  = 363  et  °-085  (MPa) 

for  0.5  mm  Nb 

(4.8) 

at  = 386  et  0-076  (MPa) 

for  0.25  mm  Nb 

The  above  true  stress-strain  relations  have  been  utilized  to  describe  the  intrinsic 
properties  of  the  unconstrained  Zr02  coated  Nb  as  well  as  the  unconstrained  AI2O3 
coated  and  uncoated  Nb  because  of  the  following  reasons.  As  shown  in  Tables  3.1 
and  4.8,  although  the  tensile  strength  of  the  niobium  varies  with  the  thickness,  the 
tensile  strength  is  propotional  to  its  hardness,  i.e.,  au/HV  equals  to  2.48  ± 0.04  and 
2.19  ± 0.08  for  the  as-received  and  processed  Nb,  respectively.  Such  relation  has 
also  been  reported  for  other  materials,  especially  for  those  that  show  a necking  at  the 
maximum  tensile  load  [180].  Ratio  of  the  yield  strength  to  the  Vicker's  hardness  is 
also  approximately  a constant,  as  shown  in  Tables  3.1  and  4.8.  The  results  indicate 
that  it  is  reasonable  to  infer  the  strength  of  the  niobium  from  the  hardness  measured. 
Hardnesses  of  Nb  hot  pressed  with  different  conditions  are  summarized  in  Table  4.9. 
The  table  shows  that  the  hardnesses  are  almost  the  same  for  each  group  of  Nb  based 
on  the  size.  Due  to  such  hardness  results,  the  stress-strain  relations  for  the 
unconstrained  uncoated  and  AI2O3  coated  Nb  can  be  assumed  to  be  the  same  as  that 
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of  the  processed  Nb  which  is  the  Zr02  uncoated  Nb  before  MoSi2  on  both  sides  of 
the  Nb  foils  has  been  broken  away. 

It  is  clear,  from  the  above  discussion,  that  the  intrinsic  mechanical  properties 
of  the  Nb  reinforcement  in  the  composites  are  different  from  those  of  the  Nb  before 
the  processing.  It  is  believed  that  such  change  in  properties  due  to  the  high 
temperatures  during  hot  compaction  processing  is  a common  phenomenon  for  ductile- 
phase-reinforced  brittle  matrix  composites.  Even  in  chemically  compatible 
composites,  the  change  of  grain  size  will  cause  property  change  of  the  ductile 
reinforcements.  Such  change  imposes  difficulty  on  the  prediction  of  the  composite 
properties. 

4.5.2  Failure  Mechanism  and  Related  Observation  on  Constrained  Niobium 

An  edge  view  of  the  unconstrained  and  constrained  uncoated  Nb  foils  after 
tensile  test  is  shown  in  Fig.  4.48.  As  seen  in  the  figure,  Nb  foils  fail  by  drawing 
down  to  a wedge  for  both  the  constrained  and  unconstrained  conditions.  All  three 
different  sizes  of  Nb  foils  show  the  same  failure  mechanism.  Furthermore,  similar 
necking  was  also  observed  for  all  the  coated  Nb  hot  pressed  at  1400^C. 

The  decohesion  lengths  (defined  in  Section  4.4.2)  measured  for  the  tensile  test 
specimens  are  summarized  in  Table  4.10.  The  data  show  that  within  each  size  group, 
decohesion  length  increases  from  the  uncoated  to  AI2O3  coated  and  then  to  Zr02 
coated  Nb/MoSi2  composites.  Such  results  are  consistent  with  the  measurement  of  the 
interfacial  fracture  energy,  as  discussed  in  Section  4.3  which  shows  that  the  uncoated 
Nb/MoSi2  has  the  highest  interfacial  fracture  energy,  followed  by  the  AI2O3  coated 
and  then  ZrC>2  coated  Nb/MoSi2  composites. 

The  increase  of  the  decohesion  length  with  the  thickness  of  Nb  foils,  as  shown 
in  Table  4.10,  is  believed  to  be  mainly  related  to  the  necking  of  the  Nb.  Figure  4.49 
shows  schematically  the  effect  of  the  necking  on  the  decohesion  length.  The  length  of 
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Fig.  4.48  An  edge  view  of  fractured  specimens,  showing  contour  of  the  neck,  (a)  is  a unconstrained,  processed  Nb  foil  with 
a thickness  of  1.0  mm,  (b),  (c)  and  (d)  are  constrained  Nb  with  thicknesses  of  1.0,  0.5  and  0.25  mm,  respectively. 
Note  that  different  magnifications  are  used. 
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the  Nb  region  affected  by  necking,  2H,  is  proportional  to  the  thickness  of  the  Nb,  as 
shown  in  Fig.  4.48.  Approximately,  H is  1.3  times  of  the  thickness  of  Nb,  measured 
from  Fig.  4.48.  Large  lateral  displacement  of  Nb  in  the  necking  affected  region  gives 
rise  to  large  transverse  stresses  which  lead  to  the  interfacial  debonding  and/or  matrix 
fracture.  Thus,  the  thicker  the  niobium,  the  longer  the  necking  affected  region,  and 
therefore  the  longer  the  decohesion  length.  If  the  interface  bonding  is  relatively  weak, 
such  as  the  AI2O3  and  Zr02  coated  interfaces,  the  transverse  stresses  aroused  by 
lateral  displacement  of  the  Nb  during  uniform  strain  can  be  large  enough  to  cause  the 
interfacial  debonding  all  the  way  to  the  test  grips.  During  the  uniform  strain,  the 
lateral  displacement  of  Nb  right  at  the  interface  is  proportional  to  thickness  of  the  Nb. 
Therefore,  decohesion  by  this  mechanism  is  easier  to  occur  in  thick  Nb  composites 
than  in  thin  Nb  counterparts,  as  indicated  in  Table  4.10. 

It  has  been  shown  in  Section  4.4  that  the  failure  of  the  constrained  Nb  laminae 
in  MoSi2  matrix  involves  the  following  scenario.  When  an  advancing  crack  just 
impinges  a Nb  reinforcement,  no  debonding  at  the  matrix/reinforcement  interface 
occurs.  As  external  load  continues  to  increase,  debonding  at  the  interface  and/or 
multiple  fracture  of  the  matrix  near  the  interface  occurs  due  to  a relatively  large 
lateral  deformation  of  the  Nb  reinforcement  compared  to  the  matrix  and  load  transfer 
from  the  matrix  to  the  reinforcement.  Then,  the  necking  and  fracture  of  the  niobium 
begin  as  the  load  continues  to  increase.  Based  on  this  observation,  the  flow  behavior 
of  the  constrained  Nb  can  be  simplified  into  three  stages:  (1)  elastic  deformation;  (2) 
plastic  deformation  within  the  decohesion  region;  (3)  localized  plastic  deformation 
within  the  necking  region.  Three  stages  are  illustrated  in  Fig.  4.50.  The  crack 
opening  of  the  matrix  can  be  directly  related  to  the  deformation  of  Nb  in  the 
decohesion  zone.  At  the  early  stage  of  crack  opening,  Nb  only  undergoes  elastic 
deformation,  as  shown  in  Fig.  4.50(a).  As  crack  opening  increases,  Nb  enters  the 
stage  of  plastic  deformation  (Fig.  4.50(b)).  At  last,  the  work  hardening  rate  of  the 
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Fig.4.49  A schematic  showing  effect  of  necking  on  the  decohesion  length. 
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necking  affected  region,  2a  is  the  remaining  thickness  of  the  reinfrocement,  and  u is 
displacement. 
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niobium  at  the  central  portion  of  the  decohesion  region  is  lower  than  the  stress 
increasing  rate  due  to  decrease  in  the  cross  section  area  of  the  reinforcement,  the 
localized  deformation  begins  (Fig.  4.50(c)),  which  leads  to  the  final  failure  of  the  Nb 
lamina.  Such  simplified  deformation  processes  will  be  used  to  model  stress- 
displacement  behavior  of  the  constrained  Nb,  as  presented  in  Section  4.5.4. 

4.5.3  Stress-Displacement  Curves  of  Constrained  Niobium 

Figure  4.51  shows  typical  stress-displacement  curves  for  the  uncoated  Nb  foils 
with  different  thicknesses  constrained  in  MoSi2  matrix.  The  parameters  measured  are 
summarized  in  Table  4.11.  The  work  of  rupture  normalized  by  the  yield  strength,  o0; 
and  half  thickness  of  ductile  phase,  t , is  given  by 

_ ru'  g(u)  d u 
‘ " J°  CT0  t 

The  values  of  Et,  called  the  normalized  work  of  rupture  in  the  text  for  convenience, 
are  also  included  in  Table  4.11.  As  seen  in  the  figure  and  table,  due  to  the  different 
intrinsic  properties  of  the  niobium,  the  maximum  stresses  reached  by  the  constrained 
Nb,  amax,  are  different  for  different  sizes  of  Nb.  However,  if  the  omax  is 
normalized  by  its  own  yield  strength,  the  normalized  maximum  stress,  amax/CTo, 
becomes  independent  of  the  thickness  of  Nb.  This  result  suggests  that  the  maximum 
stress  reached  is  not  a function  of  size  of  the  ductile  phase.  On  the  other  hand,  it  is 
noted  that  the  work  of  rupture  increases  with  size  of  Nb,  indicating  that  large  size  of 
ductile  reinforcement  is  more  effective  in  improving  toughness  of  the  brittle  matrix 
composites.  It  is  also  noted  that  the  normalized  work  of  rupture,  Et,  increases  with 
increasing  size  of  the  niobium.  The  size  dependency  of  the  normalized  work  of 
rupture  is  attributed  to  the  size  dependency  of  the  decohesion  length,  as  will  be 
discussed  in  Section  4.5.5.  The  result  suggests  that  the  data  of  the  normalized  work 
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Displacement,  mm 


Fig.  4.51  Typical  stress-displacement  curves  of  constrained,  uncoated  Nb  laminae  with 
different  thicknesses. 


Fig.  4.52  Effect  of  the  coatings  on  the  stress-displacement  curve,  measured  from 
composite  laminates  with  0.5  mm  thick  Nb  lamina. 


Table  4.1 1 Summary  of  the  parameters  calculated  from  the  model  and  measured  for  the  uncoated  Nb/MoSi2  laminates 

Measured  Calculated 
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Eight  specimens  for  each  condition  were  tested. 
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of  rupture  obtained  from  the  tests  on  large  size  of  ductile  reinforcements  cannot  be 
extended  to  small  size  of  the  reinforcements  directly. 

Effect  of  coatings  on  the  stress-displacement  curves  is  shown  in  Fig.  4.52. 
The  parameters  measured  are  summarized  in  Table  4.12.  The  data  show  that  as  the 
decohesion  length  increases,  i.e.,  as  the  constraint  on  the  reinforcement  decreases, 
the  work  of  rupture  increases.  This  is  not  surprising  since  the  longer  the  decohesion 
length,  the  more  ductile  material  participates  in  the  plastic  deformation,  and  therefore 
the  more  energy  is  consumed  before  fracture.  It  is  noted  that  as  the  decohesion  length 
increases,  the  normalized  maximum  stress  decreases,  indicating  that  increasing 
constraints  on  the  ductile  phase  increases  the  omax-  Table  4.12  also  shows  that  the 
normalized  work  of  rupture  is  dependent  on  the  interfacial  bonding  strength.  This  is 
believed  to  be  due  to  the  increase  of  work  of  rupture  with  increasing  decohesion 
length.  The  above  results  indicate  that  limited  decohesion  enhances  the  energy 
consumed  to  fracture  the  constrained  ductile  reinforcement,  and  therefore  improves 
toughness  of  the  composites. 

4.5.4  Modeling  of  Flow  Behavior  of  the  Constrained  Niobium 

As  discussed  in  Section  4.5.2,  the  flow  behavior  of  the  constrained  Nb  in  the 
present  study  can  been  simplified  into  three  stages:  (1)  elastic  deformation;  (2)  plastic 
deformation  within  the  decohesion  region;  (3)  localized  plastic  deformation  within  the 
necking  region.  To  model  the  stress-displacement  relation  of  constrained  ductile 
phase,  ct-u,  some  assumptions  are  made.  First,  displacement  of  the  stress- 
displacement  curve  is  assumed  to  only  come  from  the  deformation  of  Nb  lamina 
inside  the  decohesion  region  because  the  Nb  outside  the  decohesion  region  is  bonded 
to  the  matrix  and  only  undergoes  elastic  deformation.  Second,  contour  of  the  Nb 
lamina  in  the  decohesion  region  is  assumed  to  be  part  of  the  outside  surface  of  a 
cylinder  with  a varied  radius  of  R tor  all  the  stages  of  the  deformation,  as  shown  in 
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Fig.  4.50.  Such  an  assumption  is  a good  approximation  to  the  real  contour  of  the  Nb 
observed  in  the  experiments  and  makes  the  calculation  possible.  Third,  plane  strain  is 
assumed  for  the  computation  in  all  the  stages.  In  addition,  in  the  calculation  the  value 
of  u is  measured  as  the  axial  displacement  of  the  matrix  point  at  the  boundary  of  the 
decohesion  region,  while  a is  computed  from  the  total  force  exerted  on  the  remaining 
cross  section  of  the  midplane  at  the  neck,  normalized  by  the  original  midplane  area. 

Stage  of  elastic  deformation.  To  compute  o-u  data  for  elastic  deformation,  two 
more  assumptions  are  made:  (1)  the  volume-conserving  nature  of  large  scale  plastic 
deformation  is  also  applicable  to  this  stage;  (2)  the  effective  stress  across  the 
midplane  of  the  neck  is  constant.  Although  the  assumed  contour  and  constant 
effective  stress  do  not  represent  the  real  situation  of  the  decohesion  region,  they 
introduce  negligible  error  to  the  work  of  rupture  (see  the  Appendix). 

Based  on  the  assumption  (1)  and  the  assumed  contour,  the  displacement,  u, 
can  be  related  to  the  radius  of  the  cylinder,  R,  by  geometry  as  follows. 


- ( d + u )3  , 4 t u 

12  t u 3 ( d + u ) 


(4.9) 


From  the  assumption  of  plane  strain,  we  have 

Gy  = |i  (o2  + ox)  (4.10) 

where  p is  the  Poisson's  ratio  of  Nb  and  equals  0.39  [181]. 

By  definition,  the  effective  stress,  a , is  written  as  [182] 
a 2 = 1/2  { ( ox  - Cy  )2  + (ay  - oz  )2  + (az  - ax  )2  } (4.11) 

where  ax,  Gy,  and  oz  are  principal  stresses  because  the  shear  stress  components 
vanish  on  the  midplane  of  the  neck  by  symmetry  [183].  Substituting  eq.  (4.10)  into 
eq.  (4.11),  we  have 


o2  = 0.762  ax2  + 0.762  oz2  - 1.4758  Gxoz  (4.12) 

From  the  boundary  condition  at  the  neck,  i.e.,  ax  = 0,  eq.  (4.12)  reduces  to 


o 2 = 0.762  oza2 


(4.13) 
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where  aza  is  the  stress  in  the  z-direction  on  the  free  surface  of  the  midplane  of  the 
neck.  Substituting  eq.  (4.13)  back  into  (4.12)  and  solving  for  az,  we  have 


<JZ  = 0.968  ax  + tJo^2  - 0.062  trx 


= 0.968  ox  ± oza 


(4.14) 

because  0.062  ctx^  < < (as  shown  in  eq.  (4.16)).  The  negative  sign  in  eq. 
(4.14)  can  be  discarded  because  of  the  boundary  condition  at  the  neck  that  ax  = 0 
and  oz  = o ^ when  x = a. 

Based  on  the  force  balance  in  the  x-direction  acting  on  an  element  at  the 
midplane  of  the  neck,  Bridgman  has  derived  a formula  for  a necking  plate  [183] 

d or 


(“)  °z  =(-)<*, 

a a 


d x 


(4.15) 


where  x and  a have  been  defined  in  Fig.  4.50.  Introduction  of  eq.  (4.14)  into  (4.15) 
results  in  a linear  first  order  differential  equation  for  ax  , i.e., 


d x 


0.032  x 


a2  x2 

aR  + V (1  ' ^ ) 

2 a 


The  solution  for  eq.  (4.16)  is 


■=r  ct_  = - •=■ 


aR  + — ( 1 
2 


x2 
— ) 
a2  ) 


=r  cr. 


(4.16) 


<7  = <7 

x za 


31.25  - 31.25 


0.032 


(-) 

a 


( 


R + 0.5  - 0.5  \)0032 


(4.17) 


a a‘ 

The  above  can  be  related  to  the  average  elastic  strain,  ex  , which  can  be 
measured  from  the  dimension  of  the  neck  using  the  relation  of  ix  = In  (a/t)  and 
related  to  the  elastic  strain  in  the  x-direction  at  the  position  of  x on  the  midplane  of 
the  neck,  Sx(^x>^y>^z)>  ^ follows 
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cjz  ) dx 

ex  = L 

a 

(4.18) 

where 

ex  can 

be  found  from  the  plane  strain  condition  using  the  equations  of 

cy  = p.  ( ox  + oz  ) 

(4.19) 

and 

£x  = 1/E  { gx  - p (o 

y °z)  } 

(4.20) 

where 

E is  Young's  modulus  of  Nb. 

Substituting  eqns.  (4.14),  (4.17)  and  (4. 

19)  into 

(4.20) 

yields 

£ = 


9.55  - 10.1 


^0  032 

a 


R 

a 


(-  + 0.5  - 0.5  \)0032 


(4.21) 


Substituting  eq.  (4.21)  into  (4.18)  results  in 


£ x = 


R - 0.032 


aE 


9.55  a - 10.1  f“ 

Jo 


(-) 

a 


(-  + 0.5  - 0.5  ^-)0032 


dx 


(4.22) 


a a* 

Thus,  by  numerical  integration  of  eq.  (4.22),  ct^  can  be  found  from  ex  for 
any  specific  R and  a.  With  oza  in  hand,  az  can  be  found  from  eqns.  (4.14)  and 
(4.17). 


Stage  of  plastic  deformation.  The  Von  Mises  yield  criterion  is  used  to  monitor 
the  initial  yielding.  Thus,  when  the  effective  stress  reaches  the  yield  strength  of  the 
Nb,  plastic  deformation  begins.  Now,  Bridgman's  formula  for  a necking  plate  on  the 
condition  of  plane  strain  [183] 


a,  = a. 


1 + In  [ 1 + 0.5  — 
R 


(4.23) 


can  be  used  directly  to  compute  the  nominal  axial  stress.  By  applying  boundary 
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conditions  at  the  neck,  a relation  between  az a and  the  effective  stress,  a , can  be 
found  as  follows.  The  plane  strain  condition  at  the  stage  of  plastic  deformation  gives 
cy  = 0.5  ( cx  + oz  ) (4.24) 

Substituting  eq.  (4.24)  into  (4.11)  yields 

a 2 = 0.75  (ax-az)2  (4.25) 

Combining  eq.  (4.25)  with  the  boundary  condition  that  ox  = 0 on  the  free  surface  of 
the  midplane  of  the  neck,  we  have 

Oza  = 1.1547  a (4.26) 

where  o is  determined  by  the  effective  strain  £.  The  effective  strain  can  be 
calculated  from  the  plastic  strain  sXp  in  the  x-direction  using  the  plasticity  equation 
which  can  be  written  as  [184] 

d£xp  = iJ-  [ crx  - 0.5  ( ay  + at  ) ] (4.27) 

Substituting  eqns.  (4.24)  - (4.26)  into  (4.27)  results  in 

d exp  = - 0.866  de  (4.28) 

Integration  of  eqn.  (4.28)  yields 

£ = - 1.1547  exp  (4.29) 

where  eXp  is  constant  over  the  plane  of  the  neck  because  £ is  constant  over  the  whole 
plane  of  the  neck  [183].  In  the  computation,  the  effective  strain  at  the  initial  yield  is 
assigned  as  zero  (i.e.,  eXp  is  equal  to  the  total  strain  ln(a/t)  deducted  by  the  elastic 
strain  at  the  yield),  while  the  effective  stress  has  a value  of  the  yield  strength  of  the 
unconstrained  Nb.  After  the  initial  yield,  the  niobium  is  assumed  to  have  perfect 
isotropy  of  strain  hardening,  and  a power  law,  <7  = C £ m,  is  used  to  compute  the 
effective  stress.  Thus,  ctz  can  be  found  by  substituting  the  effective  stress  into  eqns. 
(4.23)  and  (4.26). 

Stage  of  localized  plastic  deformation.  The  localized  plastic  deformation  is 
assumed  to  begin  when  the  nominal  axial  stress  is  about  to  decrease.  The  length  of 
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the  region  affected  by  necking,  2H,  shown  in  Fig.  4.50.,  is  chosen  as  2t.  After 
necking  down  to  a point  at  the  neck,  the  region  with  such  dimension  gives  two 
triangles  with  the  height  equal  to  the  base.  The  base  is  assumed  to  have  no  more 
deformation  when  the  localized  deformation  begins.  The  assumptions  give  rise  to  a 
contour  as  shown  in  Fig.  4.50(c),  and  such  contour  is  an  approximation  for  the 
contour  of  Nb  in  the  decohesion  region,  as  shown  in  Fig.  4.48.  The  displacement  in 
the  stage  of  the  localized  deformation  is  assumed  to  come  only  from  the  deformation 
of  the  necking  affected  region  and  the  computation  of  the  nominal  axial  stress  is 
carried  out  using  eqs.  (4.23)  and  (4.26). 

In  summary,  with  the  above  approximate  model  the  effect  of  bonding 
conditions  at  the  interface  on  the  calculated  cr(u)  ~ u curve  is  introduced  by  properly 
selecting  the  decohesion  length  at  the  interface,  since  the  decohesion  length  is 
dependent  on  the  interfacial  conditions  as  revealed  in  the  experiments  (Sections  4.4.2 
and  4.5.2).  The  effects  of  the  intrinsic  mechanical  properties  of  ductile  phase  are 
addressed  by  utilizing  the  elastic  modulus,  yield  strength  and  true  stress-strain  relation 
of  the  ductile  phase  throughout  the  computation  of  the  o(u)  ~ u curve.  Finally,  the 
size  of  the  ductile  phase  comes  into  play  through  the  selection  of  the  decohesion 
length  as  well  as  the  direct  entry  of  the  half  thickness  of  the  ductile  phase  in  the 
model. 

4.5.5  Comparison  between  the  Mo  lei  and  Experiments 

Effect  of  the  decohesion  length.  Dependence  of  the  stress-displacement  curve 
on  the  decohesion  length  is  shown  in  Fig.  4.53  which  is  generated  using  the  true 
stress-strain  data  of  0.5  mm  Nb.  The  results  indicate  that  the  maximum  stress  reached 
by  the  constrained  Nb,  omax>  increases  with  increasing  constraints  on  the  ductile 
reinforcement.  However,  increasing  constraints  decreases  the  work  of  rupture.  Such 
results  suggest  that  a relatively  weak  bond  at  the  interface  enhances  the  work  of 
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rupture  and  therefore  is  conducive  lo  toughening  of  the  brittle  matrix  composites. 
Similar  results  are  obtained  by  other  investigators  using  different  models  [23,24,26]. 

The  calculated  parameters  from  the  model  based  on  the  experimental  data  of 
0.5  mm  thick  Nb  lamina  in  Table  4.10  are  presented  in  Table  4.12.  Ratios  of  the 
calculated  and  measured  normalized  work  of  rupture,  Et(cal)/Et(exp),  is  also  included 
in  the  table.  A comparison  between  the  calculated  and  measured  data  in  Table  4.12 
reveals  that  the  calculated  work  of  rupture  deviates  from  the  measured  counterparts 
by  about  20%.  However,  the  general  trend  that  work  of  rupture  decreases  with 
increasing  interfacial  fracture  energy  observed  in  the  experiments  has  been  maintained 
in  the  prediction  of  the  model.  The  discrepancy  in  the  values  between  the  calculated 
and  measured  work  of  rupture  lies  in  the  assumptions  involved  in  the  model, 
especially  the  one  relating  to  the  contour  of  the  Nb  lamina  in  the  decohesion  region. 
Although  the  assumed  contour  of  the  neck  makes  the  computation  feasible,  it  does  not 
exactly  reflect  the  real  contour  of  the  Nb  which  actually  displays  a convex  shape  at 
the  region  away  from  the  neck,  as  shown  in  Fig.  4.48.  On  the  one  hand,  the  assumed 
contour  (a  concave  shape)  leads  to  a larger  strain  to  failure,  when  compared  with  the 
convex  shape  of  the  real  contour;  on  the  other  hand,  it  also  results  in  a larger  radius 
of  curvature  at  the  neck  than  does  the  real  contour.  As  such,  a less  severe  state  of 
triaxial  stresses  results  and  a lower  stress  in  axial  direction  (Z-direction)  is  predicted 
by  the  model.  All  of  these  contribute  to  the  discrepancy  in  the  work  of  rupture. 
Nevertheless,  the  predicted  values  are  reasonable  approximations  to  the  measured 
counterparts. 

It  is  also  noted  from  Table  4.12  that  amax  predicted  by  the  model  is 
insensitive  to  the  interfacial  bonding  conditions  in  the  range  of  the  measured 
decohesion  lengths.  The  results  suggest  that  the  present  coated  and  uncoated 
Nb/MoSi2  laminates  are  in  the  range  of  low  constraints  because  in  the  high  constraint 
range  (i.e. , short  decohesion  length),  amax  is  very  sensitive  to  the  decohesion  length, 
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Fig.  4.53  Effect  of  decohesion  length  on  the  stress-displacement  curve. 


Displacement,  rr  ^> 


Fig.  4.54  Effect  of  reinforcement  size  on  the  stress-displacement  curve.  Al,  B1  and  Cl 
(solid  lines)  are  experimental  curves  of  the  composites  with  0.25,  0.5  and  1.0 
mm  thick  Nb  laminae,  respectively;  A2,  B2  and  C2  (dashed  lines)  are  results 
from  the  model  for  composites  with  0.25,  0.5  and  1.0  mm  thick  Nb  laminae, 
respectively. 
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as  shown  in  Fig.  4.53.  An  experimental  evidence  indicating  that  the  present  laminates 
are  in  the  range  of  low  constraints  is  that  the  maximum  stress  reached  by  the  laminate 
composites  with  1.0  mm  Nb  lamina  are  almost  the  same  for  different  decohesion 
lengths  ranging  from  ~ 5 to  ~ 15  pm. 

Effect  of  reinforcement  size.  A comparison  of  the  size  effect  between  the 
model  and  experiments  for  the  uncoated  Nb/MoSi2  laminates  is  shown  in  Fig.  4.54. 
The  input  data  for  the  model  are  from  eq.  (4.8),  Tables  4.8  and  4.10  with  each  size 
of  Nb  having  its  own  measured  parameters  and  properties.  The  input  of  Young's 
modulus  is  from  the  slope  of  the  stress-strain  curve  of  the  unconstrained  Nb  in  a 
simple  tensile  test.  As  shown  in  Fig.  4.54,  the  model  fits  the  experimental  data 
reasonably  well  except  the  stage  of  the  localized  plastic  deformation.  The  deviation  in 
this  stage  is  believed  to  be  due  to  the  difference  between  the  real  and  assumed 
contours  of  the  neck,  as  discussed  earlier.  The  parameters  calculated  from  the  model 
are  summarized  in  Table  4.11.  It  is  noted  that  the  calculated  maximum  stress,  the 
work  of  rupture  and  the  normalized  work  of  rupture  are  all  close  to  the  measured 
counterparts,  clearly  indicating  that  these  can  be  estimated  from  the  model. 

The  reinforcements  used  in  composites  for  structural  applications  may  be 
smaller  in  size  than  those  tested  in  the  present  study.  To  predict  the  work  of  rupture 
for  small  size  of  reinforcements,  it  is  necessary  to  know  the  decohesion  length  in 
advance  for  the  use  of  the  model.  A plot  of  omax  vs.  decohesion  length  computed 
from  the  model  for  the  composites  with  0.25  mm  Nb  lamina  is  shown  in  Fig.  4.55.  It 
is  noted  that  there  exists  a characteristic  decohesion  length,  dc,  below  which  amax 
becomes  very  sensitive  to  the  decohesion  length.  Computation  of  the  model  shows 
that  the  characteristic  decohesion  length  is  equal  to  twice  of  the  thickness  of  the 
niobium  and  the  normalized  work  of  rupture  with  the  characteristic  decohesion 
length,  called  as  the  characteristic  normalized  work  of  rupture,  is  constant  regardless 
of  size  of  the  reinforcement. 
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Fig.  4.55  A plot  of  Gmax  vs  decohesion  length.  Arrow  A shows  characteristic 

decohesion  length,  and  arrows  B,  C and  D show  the  measured  decohesion 
lengths  of  the  composites  containing  uncoated,  AI2O3  coated  and  Zrt>2  coated 
Nb  laminae,  respectively. 


Fig.  4.56  A plot  of  ratio  of  the  measured  to  characteristic  decohesion  lengths  as  a 
function  of  reinforcement  size. 


Fig.  4.57  Normalized  work  of  rupture  as  a function  of  reinforcement  size. 
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The  measured  decohesion  lengths,  dm,  in  the  present  study  are  all  larger  than 
dc.  However,  the  ratio  of  dm  to  dc  is  found  to  be  a linear  function  of  size  of  the 
reinforcement,  as  shown  in  Fig.  4.56.  The  relations  found  are 
<Wdc  = 0.66  + 10.72  t ( for  AI2O3  coated  laminates) 

dm/dc  = 0.66  + 8.70 1 ( for  uncoated  laminates ) (4.30) 

It  is  noted  that  dm  becomes  closer  to  dc  as  size  of  the  reinforcement 
decreases,  suggesting  that  the  normalized  work  of  rupture  becomes  smaller  and  closer 
to  the  characteristic  normalized  work  of  rupture.  The  oxide  coating  changes  the  slope 
of  the  line,  indicating  that  the  coating  enhances  decohesion  and  such  effect  becomes 
larger  as  the  reinforcement  size  increases.  With  the  aid  of  eq.  (4.30),  the  decohesion 
length  of  small  size  of  the  reinforcement  can  be  predicted  and  the  work  of  rupture  be 
calculated  from  the  model.  A plot  generated  in  this  way  is  shown  in  Fig.  4.57  which 
shows  that  Et  decreases  with  decreasing  the  reinforcement  size,  as  found  in  the 
experiments.  Thus,  the  model  can  be  used  to  estimate  the  contribution  of  ductile 
reinforcements  to  the  toughness  of  a brittle  matrix  composite. 

Effect  of  yield  strength  and  work  hardening.  Effects  of  yield  strength  and 
work  hardening  on  the  work  of  rupture  are  evaluated  by  substituting  the  properties  of 
several  different  materials  into  the  model.  The  properties  are  obtained  from  Reference 
[185]  and  the  true  stress-strain  curves  for  some  materials  are  shown  in  Fig.  4.58.  It  is 
noted  that  both  4340  steel  and  70/30  brass  have  a higher  work  hardening  rate  than  the 
Nb  used  in  the  present  study.  The  computed  stress-displacement  curves  are  shown  in 
Fig.  4.59,  and  the  input  and  output  data  are  summarized  in  Table  4. 13.  It  is  clear  that 
work  of  rupture  increases  with  increasing  work  hardening  rate.  Although  70/30  brass 
has  the  lowest  yield  strength,  its  high  work  hardening  rate  enables  it  to  have  a much 
higher  work  of  rupture  than  both  4340  steel  and  Nb  both  of  which  have  higher  yield 
strengths  than  the  brass.  The  reason  for  such  phenomenon  is  that  a high  work 
hardening  rate  reinforcement  has  a higher  increase  rate  in  the  load-carrying  ability, 


Table  4.13  The  input  and  output  data  for  computing  the  stress-dsiplacement  curves  shown  in 
Fig.  4.59  (the  thickness  of  all  materials  is  assumed  as  0.5  mm) 
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Fig.  4.58  The  true  stress-strain  curves  of  the  some  materials  used  for  evaluating 
effects  of  yield  strength  and  work  hardening. 
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displacement,  mm 


Fig.  4.59  Computed  stress-displacement  curves  using  the  data  in  Fig.  4.58  and 
Table  4.13. 
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and  instability  (i.e.,  the  localized  plastic  deformation  within  the  necking  region), 
comes  relatively  later  than  a low  work  hardening  rate  reinforcement.  Thus,  high  work 
hardening  rate  reinforcements  have  a relatively  longer  stage  of  plastic  deformation  in 
the  whole  decohesion  region,  leading  to  more  energy  dissipation  and  higher  work  of 
rupture.  It  is  noteworthy  that  the  brass  has  a very  high  normalized  work  of  rupture, 
which  suggests  that  the  normalized  work  of  rupture  could  probably  be  used  as  an 
indicator  of  the  efficiency  with  which  the  work  hardening  contributes  to  enhancing 
the  work  of  rupture. 

When  yield  strength  is  very  high,  such  as  0.6  % C steel  in  Table  4.13,  the 
work  of  rupture  would  be  also  very  high.  As  shown  in  Table  4.13,  0.6%  C steel  has 
a much  higher  work  of  rupture  than  70/30  brass  although  it  has  a lower  work 
hardening  rate  than  the  brass.  Contribution  of  yield  strength  to  work  of  rupture  is 
mainly  to  increase  the  maximum  stress  reached  by  the  constrained  reinforcement. 
From  the  above  discussion,  it  is  concluded  that  high  work  hardening  rate  and  high 
yield  strength  are  both  beneficial  to  the  work  of  rupture.  Work  hardening,  however, 
is  more  effective  in  enhancing  the  work  of  rupture  than  yield  strength  because  the 
former  delays  the  instability  of  deformation  for  constrained  ductile  reinforcements, 
therefore  more  material  of  the  reinforcements  participates  in  plastic  deformation  and 
more  energy  has  to  be  consumed. 

In  summary,  based  on  the  small-scale  bridging  concept,  the  experiments  and 
the  model  indicate  that  work  of  rupture  is  enhanced  by  a relatively  weak  bond  at  the 
matrix/reinforcement  interface,  by  large  size  of  reinforcement,  and  by  high  yield 
strength  and  high  work  hardening  r:;te.  High  work  hardening  rate  is  more  effective  in 
enhancing  work  of  rupture  than  high  yield  strength. 
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4.6  Effects  of  Interface.  Size  and  Orientation  of  Ductile  Phases  on  Toughness 
In  this  section  the  effects  of  interface,  size  and  orientation  of  ductile  phases  on 
the  toughness  of  composites  will  be  discussed.  For  the  procedures  of  measuring  the 
toughness  of  the  composites,  please  refer  to  Section  3.6. 

4.6.1  Dual-Effect  of  the  Matrix/Reinforcement  Interface  on  the  Toughness 
of  the  Composites 

Typical  load-displacement  curves  of  bending  tests  on  chevron  notched 
composite  laminates  hot  pressed  at  1400°C  are  shown  in  Fig.  4.60.  There  are  two 
prominent  features  in  the  figure:  (1)  uncoated  Nb  reinforced  composites  exhibit  the 
highest  peak  load,  (2)  mechanical  behavior  of  Zr02  and  AI2O3  coated  systems  are 
similar,  both  of  them  showing  an  increase  in  the  carried  load  in  the  last  part  of  the 
displacement  curves.  The  latter  observation  is  caused  by  the  extensive  delamination  at 
the  interface,  and  correspondingly  more  niobium  participates  in  deformation  and 
deforms  under  much  less  constrained  condition.  Because  of  this  feature,  the  total 
energy  consumed  to  break  a specimen  (area  under  the  curve)  for  the  coated  composites 
is  larger  than  the  uncoated  ones.  The  work  of  fracture,  defined  as  the  total  energy 
normalized  with  respect  to  the  generated  crack  area  [186],  for  the  various  composites 
are  presented  in  Table  4.14.  As  seen  in  the  table,  the  coated  composites  show  the 
higher  work  of  fracture,  indicating  that  the  low  interfacial  fracture  energy,  and  thus 
long  decohesion  length,  is  beneficial  to  improving  the  toughness  of  the  composite.  This 
result  is  consistent  with  predictions  based  on  other  numerical  and  analytical  studies 
[23,24,26]. 

However,  the  damage  tolerance  determined  from  the  peak  load  of  the  chevron 
notched  specimens  shows  an  opposite  trend.  The  damage  tolerance  determined  in  this 
way  is  also  included  in  Table  4.14.  The  data  show  that  uncoated  Nb  reinforced 
composites  exhibit  the  highest  damage  tolerance,  followed  by  AI2O3  coated  and  then 


Table  4.14  Measured  damage  tolerance,  steady-state  toughness  and  other  relevant  parameters  from  the 

tensile  tests  and  bend  tests. 
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a.  From  bend  tests  on  chevron  notched  composite  laminates  reinforced  with  20  vol.%  of  Nb  laminae  (0.25  mm  thick). 

b.  For  details,  see  Tables  4.5  and  4.7. 

c.  From  tensile  tests  on  a single  constrained  Nb  (0.25  mm  thick). 

#.  Four  specimens  for  each  condition  have  been  tested  except  AI2O3  coated  Nb  system  for  which  eight  specimens  have 
been  tested. 
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Fig.  4.60  Typical  load  displacement  curves  of  chevron  notched  specimens  for  laminated 
composites  containing  20  vol%  of  Nb  laminae  with  different  interfacial 
conditions,  hot  pressed  at  1400®C  for  1 hour.  The  origin  for  each  curve  has 
been  shifted  for  the  convenience  of  observation. 


Fig.  4.61  Typical  load-displacement  curves  of  chevron  notched  specimens  for  monolithic 
MoSi2  and  laminated  composites  containing  different  thicknesses  of  uncoated 

Nb  laminae,  hot  pressed  at  140fPC  for  1 hour.  The  origin  for  each  curve  has 
been  shifted  for  the  convenience  of  observation. 
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Zr02  coated  composites.  This  sequence  is  parallel  to  the  increase  in  the  interfacial 
fracture  energy  of  the  composites,  that  is,  the  higher  the  interfacial  fracture  energy,  the 
higher  the  damage  tolerance  of  the  composites. 

The  data  in  Table  4. 14  also  indicate  that  loss  of  ductility  of  the  ZrC>2  coated  Nb 
foils  hot  pressed  at  1700°C  results  in  a large  drop  in  the  toughness  of  the  laminated 
composites.  Therefore,  the  damage  tolerance  data  suggest  that  a high  degree  of 
constraint  is  conducive  to  toughness  as  long  as  the  ductile  reinforcement  has  a high  slip 
capability.  When  the  reinforcement  loses  its  slip  capability,  as  in  the  case  of  Zr02 
coated  composites  processed  at  1700°C,  decreasing  the  constraint,  i.e.,  increasing  the 
gauge  length  by  decreasing  the  interfacial  bonding,  is  more  beneficial  to  improving  the 
toughness. 

In  summary,  the  results  from  the  bend  test  of  chevron  notched  specimens 
indicate  that  the  role  of  the  interface  for  ductile  reinforcement  and  debonding  in  the 
toughness  depends  on  the  criterion  used  to  describe  the  toughness  of  the  composites.  If 
the  peak  load  of  the  chevron  notched  specimen  is  used  as  an  indicator  of  the  toughness, 
a strong  bonding,  and  therefore  a high  degree  of  constraints,  would  be  desirable.  On 
the  other  hand,  if  the  total  energy  consumed  to  break  a specimen  is  used  as  indicative 
of  its  toughness,  then  a relatively  weak  bonding  is  required. 

The  aforementioned  dual-effects  of  the  interfacial  properties  on  the  toughness  of 
the  composites  can  also  be  related  to  the  data  obtained  from  the  simple  tensile  tests  on  a 
single  constrained  Nb  foil  which  has  been  discussed  in  Section  4.5.  For  convenience, 
the  properties  measured  from  the  tensile  tests  on  0.25  mm  Nb  reinforced  composites 
have  also  been  included  in  Table  4.14.  The  steady-state  toughness  of  the  composites  in 
the  table  is  calculated  from  eq.  (2.6)  and  is  proportional  to  the  work  of  rupture  of  the 
constrained  ductile  phase.  It  can  be  seen  from  Table  4.14  that  increasing  decohesion 
length  is  conducive  to  improving  steady-state  toughness  of  the  composites,  a trend  also 
shown  by  the  work  of  fracture  measured  from  chevron-notched  specimens.  Table  4. 14, 
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at  the  same  time,  also  reveals  that  the  peak  stress  reached  by  the  constrained  Nb 
decreases  with  increasing  decohesion  length,  a trend  similar  to  that  exhibited  by  the 
damage  tolerance  measured  from  the  chevron  notched  specimens.  This  apparently  dual 
role  of  the  interface  can  be  explained  by  relating  the  maximum  stress  reached  by  the 
constrained  ductile  phases  to  the  crack  propagation  resistance  in  the  case  of  large-scale 
bridging  (i.e.,  crack  length  is  at  the  same  order  of  magnitude  as  the  bridging  length)  as 
observed  in  the  present  composites.  In  Section  4.7,  the  large-scale  bridging  will  be 
treated,  and  the  results  indicate  that  the  maximum  crack  propagation  resistance  of  the 
composites  increases  with  increasing  decohesion  length.  However,  the  maximum  crack 
propagation  resistance  is  achieved  only  after  the  crack  has  propagated  extensively. 
When  crack  size  is  small,  the  composites  with  less  decohesion  show  a higher  crack 
propagation  resistance  than  the  counterparts  with  more  decohesion.  Therefore,  it  seems 
that  the  choice  of  weak  or  strong  interface  in  describing  the  fracture  behavior  of 
ductile-phase-toughened  composites  depends  on  the  criterion  used  to  describe  the 
toughness  of  the  composites. 

4.6.2  Effects  of  Niobium  Size  on  the  Toughness 

Representative  load-displacement  curves  of  chevron  notched  specimens  for 
monolithic  MoSi2  and  laminated  composites  containing  different  thicknesses  of  Nb 
laminae  are  shown  in  Fig.  4.61.  The  origin  for  each  curve  in  this  figure  has  been 
shifted  for  the  convenience  of  observation.  It  is  noted  that  the  catastrophic  failure  of 
MoSi2  has  changed  to  graceful  failure,  and  the  area  under  the  load-displacement  curve 
as  well  as  the  peak  load  have  increased  by  the  incorporation  of  ductile  Nb  laminae  into 
the  matrix.  A comparison  between  laminated  composites  reveals  that  the  area  under  the 
load-displacement  curve  increases  with  increasing  thickness  of  Nb  laminae.  The  peak 
load  reached  by  the  composites  also  exhibits  a similar  trend,  i.e.,  as  thickness  of  Nb 


142 


laminae  increases,  the  peak  load  increases,  with  an  exception  of  0.5  mm  Nb  laminae 
which  shows  a similar  peak  load  as  0.25  mm  Nb  laminae. 

The  measured  work  of  fracture  and  damage  tolerance  are  summarized  in  Table 
4.15.  It  is  clear  from  the  table  that  the  work  of  fracture  of  the  composites  increases 
with  increasing  thickness  of  Nb  laminae.  Similar  size  effects  have  also  been  observed  in 
Al203/Ni  [130],  WC-Co/Mo  [131]  and  NiAl/Nb  [132]  composite  systems.  The 
damage  tolerance,  however,  does  not  exhibit  a monotonic  variation  with  lamina 
thickness.  It  is  believed  that  this  is  due  to  the  combined  effects  of  size  as  well  as 
intrinsic  properties  of  Nb  laminae.  As  mentioned  in  Section  4.4,  the  peak  load  of  the 
bend  test  is  mainly  carried  by  Nb  laminae.  The  damage  tolerance,  therefore,  is  actually 
an  indicator  of  the  bridging  capability  of  the  ductile  phase  and  should  be  related  to  size, 
intrinsic  properties  and  fracture  behaviors  of  Nb  laminae. 

Overall  views  of  representative  fracture  surfaces  for  composites  containing 
different  size  Nb  laminae  are  provided  in  Fig.  4.62.  Higher-magnification  views  of  the 
fracture  surfaces  of  different  thicknesses  of  Nb  laminae  are  presented  in  Fig.  4.63. 
The  crack  propagation  is  from  the  top  to  bottom  in  Fig.  4.63.  The  positions  from 
where  the  fractographs  of  Fig.  4.63  are  taken  are  indicated  in  Fig.  4.62  with  the  a,  b,  c 
and  d in  Fig.  4.62  corresponding  to  the  figure  designation  (a),  (b),  (c)  and  (d)  in  Fig. 
4.63.  It  is  noted  that  the  fracture  mode  of  Nb  laminae  changes  from  elongated  dimples 
(Fig.  4.63(a))  to  quasi-cleavage  (Fig.  4.63(c)  and  (d))  as  the  thickness  of  Nb  laminae 
decreases.  The  change  in  the  fracture  mode  with  the  thickness  of  the  niobium  laminae 
can  be  related  to  the  change  of  Nb  ductility.  The  microhardness  and  yield  strength  of 
hot  pressed  Nb  laminae  have  been  included  in  Table  4.15.  As  shown  in  Table  4.15,  the 
microhardness  and  yield  strength  of  the  Nb  laminae  increases  with  decreasing 
thickness.  Thus,  the  thinner  the  Nb  lamina,  the  less  ductile  it  is. 

Based  on  the  above  observations,  it  is  clear  that  the  damage  tolerance  of  the 
composites  in  the  present  study  is  a reflection  of  the  combined  effects  of  the  thickness 


Table  4.15.  Mechanical  properties  measured  from  MoSi2/Nb  composite  laminates  showing  the  effects  of  size  and 

orientation  of  the  ductile  phase  (niobium) 
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* 3,  3,  9,  5 and  3 specimens  have  been  tested  for  laminates  containing  0.125  (T-orientation),  0.125  (N-orientation),  0.25, 
0.5  and  1.0  mm  Nb  laminae,  respectively. 

# Three  specimens  for  each  size  of  Nb  lamina  have  been  tested. 

@ Go  is  not  measured  directly,  but  obtained  by  assuming  G(/HV  = 1.80. 
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Fiu  4.62  Overall  views  of  typical  fracture  surfaces  of  the  laminated  composites,  (a)  reinforced  with  1.0  mm  thick  Nb  laminae, 
' (b)  reinforced  with  0.5  mm  thick  Nb  laminae;  (c)  reinforced  with  0.25  mm  thick  Nb  laminae;  (d)  reinforced  with  0.1 

mm  thick  Nb  laminae.  The  a,  b,  c and  d in  the  figures  indicate  the  positions  where  Fig.  4.63  has  been  taken. 
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and  the  intrinsic  properties  of  Nb  laminae.  An  attempt  has  been  made  to  separate  these 
two  combined  effects,  and  the  results  will  be  discussed  in  Section  4.6.5. 

4.6.3  Effects  of  Orientation  of  Niobium  Lamina  on  the  Toughness 

The  two  different  orientations  of  chevron  notches  are  presented  schematically  in 
Fig.  4.64.  For  the  sample  shown  in  Fig.  4.64(a),  crack  propagation  in  the  matrix  is 
normal  to  the  niobium  laminae.  Therefore,  this  orientation  is  called  N-orientation.  For 
the  sample  shown  in  Fig.  4.64  (b),  the  crack  propagation  is  transverse  to  the  niobium 
laminae.  Thus,  it  is  called  T-orientation.  Typical  load-displacement  curves  of  chevron- 
notched  specimens  for  two  different  orientations  of  the  composites  containing  0.127 
mm  thick  Nb  laminae  are  shown  in  Fig.  4.65.  Again,  the  origins  of  the  curves  in  the 
figure  has  been  shifted  to  facilitate  the  comparison.  The  corresponding  damage 
tolerance  and  work  of  rupture  for  the  two  orientations  have  been  summarized  in  Table 
4.15.  It  is  clear  form  the  table  that  T-orientation  not  only  has  a larger  work  of  fracture 
but  also  has  a larger  damage  tolerance  than  N-orientation. 

A general  view  of  the  fracture  surfaces  of  chevron  notched  specimens  for  the 
two  orientations  is  presented  in  Fig.  4.66.  A relatively  higher  magnification  of  the 
fracture  surface  is  shown  in  Fig.  4.67  to  reveal  the  general  characteristics  of  crack 
propagation  for  both  orientations.  The  crack  propagation  in  Fig.  4.67  is  from  the  top  to 
bottom.  For  the  N-orientation,  the  crack  propagation  in  the  matrix  is  discontinuous  in 
nature,  as  mentioned  before.  In  contrast,  for  the  T-orientation  the  crack  propagates 
continuously  once  initiated,  making  less  matrix  fracture  and  less  tortuosity  on  the 
fracture  surface  of  the  matrix  (Fig.  4.67).  It  is  expected  that  the  discontinuous  crack 
propagation  is  more  energy  consuming  process  than  the  continuous  crack  propagation. 
However,  T-orientation  (continuous  crack  propagation)  actually  exhibits  a higher 
damage  tolerance  and  work  of  fracture  than  N-orientation,  as  shown  in  Table  4.15.  The 
apparent  contradiction  could  be  explained  in  terms  of  the  observation  that  for  the 
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(a) 


(b) 


Fie  4 64  (a)  The  notched  plane  of  N-orientation  where  crack  propagates  normal  to  the 
^ lamina  plane;  (b)  The  notched  plane  of  T-onentation  where  crack  propagates 

transverse  to  the  lamina  plane.  Crack  propagates  from  the  bottom  to  top  in  the 

figures. 


displacement,  mm 


Fie  4 65  Representative  load-displacement  curves  of  chevron  notched  specimens 
for  two  different  orientations  of  the  composites  containing  0.1-/  mm 
thick  Nb  laminae.  The  origins  of  the  curves  have  been  shifted  for  the 
convenience  of  comparison. 
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Fig.  4.66  Fracture  surfaces  of  the  laminated  composites,  (a)  Crack  propagates  normal  to 
the  lamina  plane  (N-orientation);  (b)  Crack  propagates  transverse  to  the  lamina 
plane  (T-orientation). 
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Fig.  4.67  A relatively  higher  magnification  view  of  the  fracture  surfaces  of  the  laminated 
composites,  crack  propagating  from  the  top  to  bottom,  (a)  N-orientation;  (b) 
T-orientation. 
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Fig.  4.68  Fracture  surfaces  of  the  Nb  laminae  in  laminated  composites.  The  direction  of 
the  overall  crack  propagation  is  from  the  top  to  bottom,  (a)  N-orientation;  (b) 
T-orientation. 
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present  specimens  the  crack  propagation  in  the  matrix  actually  occurs  at  the  low  level 
of  load,  and  the  peak  load  and  the  area  after  the  peak  load  is  controlled  by  the 
deformation  and  fracture  of  Nb  reinforcements.  Thus,  the  crack  propagation 
characteristics  of  the  matrix  is  not  a dominant  factor  in  determining  the  toughness  of 
the  present  composites. 

Fracture  surfaces  of  Nb  laminae  for  the  N-  and  T-orientations  are  shown  in  Fig. 
4.68.  The  fracture  surfaces  of  Nb  laminae  reveals  that  fracture  modes  for  the  two 
orientations  are  very  similar,  both  exhibiting  quasi-cleavage  fracture.  The  dependence 
of  the  measured  toughness  on  the  lamina  orientation,  therefore,  might  lie  on  the 
differences  in  the  bending  contributions  of  Nb  laminae  in  the  different  orientations, 
which  will  be  further  discussed  in  the  following  sections. 

4.6.4  Stress  Analysis  of  the  Chevron  Notched  Specimens  of  the  MoSn/Nb 

Laminated  Composites 

The  discussion  in  Section  4.4  has  clearly  shown  that  the  extensive  crack 
propagation  in  the  matrix  has  occurred  at  the  low  level  of  load  and  the  peak  load  is 
mainly  carried  by  Nb  laminae  in  the  bend  test.  Based  on  such  observation,  a simple  and 
approximate  model  has  been  proposed  to  analyze  the  present  specimens.  The  basic 
purpose  of  the  model  is  to  relate  the  bridging  capability  of  the  Nb  laminae  to  the  peak 
load  of  the  bend  test.  Assumptions  involved  in  the  present  mechanics-of-material 
approach  are  (1)  one-dimensional  stress  is  assumed  for  chevron  notched  bend  bar;  (2) 
all  sections  that  are  far  from  the  notch  plane  and  are  initially  plane  and  perpendicular  to 
the  axis  of  the  beam  remain  plane  and  perpendicular  to  it  after  bending;  (3)  all 
longitudinal  elements  bend  into  concentric  circular  arcs;  (4)  near  the  peak  load  the 
matrix  in  the  tensile  portion  has  already  failed  and  the  load  is  carried  only  by  the 
bridging  ductile  phases;  (5)  the  deformation  in  the  tensile  portion  only  comes  from  the 
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Fig.  4.69  Schematic  of  deformed  bend  bar,  showing  the  relevant  parameters  used  in  the 
analysis. 
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ductile  phases  which  are  inside  the  "decohesion  region"  because  other  ductile  phases 
are  bonded  to  the  matrix  and  only  undergo  elastic  deformation. 

Figure  4.69  shows  the  relevant  parameters  used  in  the  model.  The  hatched  area 
in  the  tensile  portion  of  the  beam  denotes  the  bridging  region.  The  dashed  line  indicates 
the  real  neutral  axis  which  changes  gradually  near  the  notch  plane  from  the  midplane  of 
the  beam  to  the  position  close  to  the  bottom  of  the  notch.  The  position  of  the  neutral 
axis  right  at  the  notch  plane  varies  during  the  bending,  because  in  the  tensile  portion  of 
the  notch  plane  load  is  carried  only  by  Nb  ligaments,  and  the  ligaments  enter  the  stage 
of  plastic  deformation,  while  the  compression  portion  (containing  MoSi2  matrix  and 
Nb  laminae)  still  undergoes  elastic  deformation.  To  simplify  the  calculation,  the  neutral 
axis  right  at  the  notch  plane  has  been  chosen  as  the  neutral  axis  for  the  beam  and  is 
indicated  as  the  dotted  line  in  Fig.  4.69.  This  is  a reasonable  approximation  since  the 
moment  and  force  balance  are  calculated  right  at  the  notch  plane.  Thus,  the  position  of 
the  neutral  axis  (N.A.)  is  determined  from  force  balance  at  the  x-axis  right  at  the  notch 
plane,  i.e., 

J £ Cr«m)  dy  dz  + J j°(iyi)  <T  <tel)  dydz  =0  (4.31) 

where  crx(com)  is  the  stress  component  at  x-direction  in  the  compressive  portion, 
CTx(ten)  is  the  nominal  stress  on  the  bridging  ductile  phases  at  x-direction  in  the  tensile 
portion,  h and  y'  are  defined  in  Fig.  4.69,  and  z direction  is  perpendicular  to  the  x-y 
plane. 

Based  on  the  geometry  and  the  assumptions  mentioned  above,  it  can  be  shown 
that  the  engineering  strain  in  the  compressive  portion,  ec,  is 

£.  = - | (4-32) 

and  the  nominal  strain  inside  the  decohesion  region  in  the  tensile  portion,  e^,  is 


(4.33) 
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where  R is  the  radius  of  curvature  of  the  neutral  axis,  L the  outer  span  of  the  bend  test 


setup,  and  d0  the  decohesion  length  at  the  matrix/reinforcement  interface.  ax(com)  in 
eq.  (4.31)  is  determined  from  the  corresponding  compressive  strain  in  eq.  (4.32) 
according  to  the  elastic  stress-strain  relations  for  MoSi2  and  Nb.  ax(ten)  is  calculated 
from  the  corresponding  strain  in  eq.  (4.33).  The  relation  between  ax(ten)  and  st  is 
calculated  based  on  the  equations  derived  in  Section  4.5.  Since  at  the  beginning  of  the 
calculation  the  position  of  the  neutral  axis  is  not  known  for  a specific  value  of  R,  the 
values  of  Y in  eqns.  (4.32)  and  (4.33)  are  unknown.  An  iteration  approach  has  been 
used  to  calculate  the  position  of  the  neutral  axis  by  first  assuming  the  midplane  of  the 
beam  as  the  neutral  axis,  then  calculating  the  corresponding  ec  and  et,  and  substituting 
them  into  eq.  (4.31)  to  check  whether  or  not  eq.  (4.31)  is  satisfied.  If  eq.  (4.31)  is  not 
satisfied,  a new  neutral  axis  closer  to  the  bottom  of  the  notch  than  the  former  neutral 
axis  is  assumed,  and  the  new  ec  and  are  calculated,  and  eq.  (4.31)  checked  again. 
The  procedure  is  repeated  until  the  absolute  values  of  the  two  terms  on  the  left  hand 
side  of  eq.  (4.31)  are  within  ± 5%  of  each  other. 

Once  the  neutral  axis  for  one  specific  value  of  R has  been  found,  then  it  is 
straightforward  to  find  the  bending  force,  P,  using  the  following  equation 


Note  that  for  each  R value,  there  is  only  one  corresponding  bending  force.  By 
decreasing  R gradually,  therefore,  a maximum  bending  force  can  be  found  which  is 
regarded  to  be  equivalent  to  the  peak  load  of  the  load-displacement  curve  of  chevron 
notched  specimens.  Comparison  between  this  simple  and  approximate  model  and 
experiments  is  made  in  the  following  section. 


(4.34) 
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4.6.5  Comparison  between  the  Experiments  and  Analytical  Results 

Table  4.16  compares  the  measured  and  calculated  damage  tolerances,  designated 
as  Kmax  in  the  table,  for  N-orientation  laminated  composites  reinforced  with  0.5  mm 
or  0.25  mm  thick  Nb  laminae.  Since  the  samples  tested  are  not  exactly  the  same  in 
dimension,  the  calculation  is  carried  out  on  the  basis  of  the  dimensions  of  each 
specimen.  As  mentioned  before,  different  thicknesses  of  Nb  laminae  have  different 
intrinsic  mechanical  properties  (eq.  (4.8)).  Thus,  the  relations  between  ox(ten)  and  st 
are  different  for  0.5  and  0.25  mm  Nb  laminae.  As  can  be  seen  from  the  table,  the 
calculated  damage  tolerances  are  about  15%  lower  than  the  measured  data.  The 
discrepancy  probably  comes  from  the  assumption  (1)  in  Section  4.6.4  because  the  state 
of  stress  near  the  notch  plane  is  3-dimensional.  The  position  dependence  of  decohesion 
length  at  the  interface  may  also  contribute  to  the  deviation  of  the  calculated  damage 
tolerance  from  the  measured  data.  Nevertheless,  the  simple  model  reflects  the  main 
feature  of  the  ductile-phase-reinforced  composites,  that  is,  the  peak  load  of  the 
composites  is  related  to  the  bridging  capability  and  the  size  of  the  ductile  ligaments. 

By  assuming  that  different  sizes  of  Nb  laminae  have  the  same  intrinsic 
mechanical  properties,  the  model  can  be  used  to  predict  the  effect  of  size  of  Nb  laminae 
on  the  toughness  of  the  composites.  One  of  the  results  calculated  in  this  way  has  been 
presented  in  Fig.  4.70  in  which  it  is  assumed  that  all  different  sizes  of  Nb  laminae  have 
the  same  properties  as  that  measured  from  0.25  mm  Nb  laminae.  It  is  clear  from  Fig. 
4.70  that  the  damage  tolerance  of  the  composites  increases  with  increasing  size  of  Nb 
laminae.  Thus,  It  can  be  concluded  that  in  the  thickness  range  studied  (from  0.1  to  1.0 
mm),  as  Nb  lamina  thickness  increases,  both  the  damage  tolerance  and  the  work  of 
fracture  of  the  composites  increase. 

Since  the  nominal  stress  on  the  bridging  Nb  as  a function  of  the  crack  opening 
in  the  T-orientation  is  not  known,  the  present  model  cannot  be  used  directly  to  compare 
with  the  experiments  at  the  T-orientation.  However,  as  discussed  above,  the  damage 
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Thickness  of  Nb  laminae,  mm 


Fig.  4.70  Effect  of  Nb  lamina  thickness  on  toughness  of  the  composites,  calculated  from 
the  proposed  model  assuming  that  different  thicknesses  Nb  laminae  have  the 
same  intrinsic  properties. 


Table  4.17  Measured  toughness  of  the  composites  containing  Nb  laminae  with  a thickness  of  0.127  mm  and 
moment  of  inertia  of  Nb  cross  section  about  the  neutral  axis  at  the  notch  plane. 
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tolerance  of  the  composites  is  related  to  bridging  capability  of  the  ductile  laminae. 
Therefore,  it  is  believed  that  the  observed  orientation  dependency  of  the  damage 
tolerance  is  due  to  the  different  bending  contributions  of  Nb  laminae  in  the  different 
orientations.  Since  the  moment  of  inertia  of  Nb  cross  section  about  the  neutral  axis  at 
the  notch  plane  reflects  the  contribution  of  Nb  laminae  to  resisting  the  external  bending 
moment  at  the  stage  of  elastic  deformation,  the  moment  of  inertia  of  Nb  cross  section  is 
calculated  and  listed  in  Table  4.17.  Note  that  both  the  ratio  of  the  damage  tolerances  in 
the  N-orientation  to  the  T-orientation  and  the  ratio  of  the  work  of  fracture  in  the  N- 
orientation  to  the  T-orientation  are  very  close  to  the  ratio  of  the  moment  of  inertia  of 
Nb  cross  section  at  the  N-orientation  to  the  T-orientation,  indicating  that  the  orientation 
dependence  of  damage  tolerance  is  indeed  due  to  the  different  bending  contributions  of 
Nb  laminae  in  the  different  orientations. 

4.7  Modeling  of  Ductile  Phase  Toughening  — Large  Scale  Bridging 

As  discussed  in  Section  2.3,  ductile  phase  bridging  has  been  analyzed 
theoretically  by  many  investigators.  The  primary  mechanism  responsible  for  the 
toughening  has  been  ascribed  to  the  bridging  of  ductile  ligaments  behind  the  advancing 
crack  tip  [20,21,23-26],  In  these  studies,  the  toughening  effect  has  been  evaluated  on 
the  basis  of  the  assumption  of  the  small-scale  bridging  (SSB)  for  which  the  bridge- 
length  is  small  relative  to  the  crack  length  and  the  specimen  size.  However,  in  most  of 
the  ductile-phase-toughened  composites,  the  bridging  length  is  often  at  the  same  scale 
as  the  crack  length  and  specimen  size  [116,133].  Our  present  research  results  also 
demonstrated  that  MoS^/Nb  composites  is  large-scale  bridging  (LSB). 

Ravichandran  [134,135]  and  Cox  [136,137]  have  recently  proposed  models  for 
describing  large  scale  bridging  to  account  for  the  effects  of  the  specimen  geometry  and 
crack  length.  Their  approaches  were  intended  to  give  a better  estimation  of  the  ductile 
phase  toughening.  However,  the  effect  of  size  of  ductile  phases  has  not  been  addressed 
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in  their  studies  and  the  ductile  phase  constitutive  relations  are  assumed  artificially  [134- 
137],  In  this  section,  therefore,  not  only  the  size  effect  will  be  addressed  but  also  the 
constitutive  relation  of  the  ductile  reinforcement  based  on  the  measurement  from 
MoSi2/Nb  composites  will  be  included  to  give  insight  into  the  direct  relation  between 
the  constituents'  properties  (fracture  energy  of  the  matrix/reinforcement  interface, 
decohesion  at  the  interface,  size  of  the  ductile  phase,  etc.)  and  the  overall  performance 
of  the  composites. 

4.7.1  Basic  Model 

The  problem  is  defined  in  Fig.  4.71  which  shows  a crack  of  length  2a  with 
ductile  phase  bridging  in  an  infinite  body.  The  hatched  areas  represent  the  bridging 
ligaments  and  the  crack  is  assumed  to  grow  symmetrically  along  the  x-axis  at  the  both 
positive  and  negative  directions,  a*  and  Oj  designate  the  externally  applied  stress  and 
bridging  stressess,  respectively.  For  the  present  study,  the  effective  stress  intensity 
acting  at  the  crack  tip  is  treated  as  a superimposition  of  an  opening  stress  intensity 
caused  by  the  externally  applied  stress  and  a stress  intensity  of  opposite  sign  caused  by 
the  clamping  forces  of  the  bridging  ligaments.  Due  to  the  nature  of  the  bridging 
analysis  and  for  the  sake  of  simplicity,  it  is  assumed  that  during  the  propagation,  the 
crack  tip  does  not  stop  inside  the  ductile  ligaments  nor  at  the  interface,  but  in  the 
matrix.  Thus,  the  effective  stress  intensity  acting  at  the  crack  tip,  Ktjp,  can  be  written 
as 

Ktip  — Kext  - Kbj-jdgg  (4.35) 

where  Kext  and  K^^ge  are  the  stress  intensity  factors  induced  by  the  external  stress 
and  bridging  tractions,  respectively.  Rearranging  eqn.  (4.35),  we  have 

Kext  = K^p  + Kjjj-j^gg  (4.36) 

The  crack  propagates  when  Ktip  = Km  or  Kext  = Kg  where  Km  and  Kg  are  the 
fracture  toughness  of  the  matrix  and  composite,  respectively.  For  a given  matrix,  Km 
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is  a constant.  Thus  the  crack-growth  resistance  of  the  composite  varies  with  K^dge 
(eqn.  4.36).  To  obtain  the  effective  stress  intensity  factor,  a stress  function  proposed  by 
Sedov  [187]  satisfying  the  present  boundary  conditions  is  utilized.  The  function  has 
been  given  as: 


<t> ' = 


n yf 


2 2 

z - a 


f 

•'-a 


geo 

z - Z 


(4.37) 


where  a is  defined  in  Fig.  4.71,  g(£)  = a*  - (cq  + 02  + ...  + cn)  and  z = x + iy. 
Following  Westergaard  [188],  the  Sedov's  stress  function  can  be  related  to  the  effective 
stress  intensity  factor  (mode  I)  and  the  equilibrium  crack  opening  field  along  the  y-axis, 
v(x,y),  by  the  following  equations 

oy(x,y)  = Re  <J>'  + y Im  0"  (4.38) 


K 


Up 


= Lim  a (r,  6)  V 2;rr  , when  0 = 0 


and 


2Gv(x,y)  = K + 1 Im  0 - y Re  0’ 


(4.39) 


(4.40) 


where  r and  0 are  the  polar  coordinate  with  the  origin  at  the  crack  tip,  k = 3 - 4v,  v 


2 j 1' 

the  Poisson's  ratio,  G the  shear  modulus,  0 = [ 0 dz  and  0 = — — 

* d z 


K*  =Kext 


and 


Substituting  eqns.  (4.37)  and  (4.38)  into  (4.39)  and  (4.37)  into  (4.40),  we  have 

(4.41) 


=T  I.  ( ai  + °2  + - + ^ 


V^a 


a - £ 


v(x)  = v“<x)  ' Im  iff.  < a'  +a>  + - + ff. ) ; 7- ! 

^ v x - a (x  - 5) 


(4.42) 


where  and  vext(x)  are  the  contributions  of  the  external  stress  to  the  effective  stress 
intensity  factor  and  the  equilibrium  crack  face  displacement,  respectively,  and 
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Fig.  4.7 1 A schematic  of  an  internal  bridged  crack. 


Crack  Extension,  mm 


Fig.  4.72  Typical  R-curve  behavior  of  ductile-phase-toughened  composites.  Calculation 
parameters:  The  fracture  toughness  of  the  matrix  is  3.3  MPa.m^2 ; Thickness  of 
the  Nb  lamina  20  (i  m;  Decohesion  length  at  the  interface  5 pm;  Volume  fraction 
of  Nb  20  vol.%.  The  unit  of  the  Ktip  is  MPa.m1/2  . The  detailed  crack  opening 
profile  and  bridging  tractions  at  points  A,  B and  C in  the  figure  are  presented  in 
Fig.  4.73. 
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VextW  = 


k + 1 V a2  - x2 


4 G a 

The  second  terms  in  the  right  hand  side  of  eqns.  (4.41)  and  (4.42)  are  the  contributions 
of  the  bridging  ligaments  to  the  effective  stress  intensity  factor  (designated  as  K^^dge) 
and  to  the  equilibrium  crack  face  displacement  (designated  as  vj5r[dge(x))>  respectively. 
From  Fig.  4.71,  it  is  clear  that  the  bridging  stress  is  a function  of  the  bridging 
ligament's  position,  x.  Thus,  by  carrying  out  the  integration  in  eqns.  (4.41)  and  (4.42), 
the  bridging  contributions  to  the  effective  stress  intensity  factor  and  the  equilibrium 
crack  face  displacement  can  be  rewritten  as 


K. 


(4.43) 


(4.44) 
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and 


^ bridge  (X)  = ~~~  ^ " X2  ^ O (X; ) jsin^  ) . Sm'(  )} 

4G  w JT,  I a a J 

(4.46) 

where  xj  is  the  position  of  the  center  of  the  ith  bridging  ligament,  and  t is  the  half 
thickness  or  radius  of  the  bridging  ligament  for  the  ductile  plate  and  fiber,  respectively. 
Substituting  eqns.  (4.43)  - (4.46)  into  (4.41)  and  (4.42),  the  effective  stress  intensity 
and  the  crack  face  displacement  profile  can  be  calculated  for  any  combination  of  the 
externally  applied  stress  and  bridging  tractions. 


4.7.2  Computation  of  the  Crack  Growth  Resistance 

To  utilize  eqns.  (4.41)  and  (4.42),  the  bridging  stresses  across  the  crack  face 
must  be  determined  from  the  equilibrium  crack  face  displacement  profile  which  in  turn 
depends  on  the  bridging  and  external  stresses.  Thus,  an  approach  similar  to  the  iterative 
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approach  [134,135,189]  has  been  taken  to  obtain  the  self-consistent  crack  face 

displacement  profile.  For  the  propose  of  illustration,  the  procedure  to  find  the  crack- 

growth  resistance  for  a given  crack  of  length  2a  are  outlined  below. 

1.  Choose  an  external  stress,  a*1,  slightly  larger  than  the  fracture  stress  for  the 

unreinforced  matrix  with  a crack  of  length  2a. 

2.  Compute  the  crack  opening  profile,  vext1(x),due  to  the  external  stress  a*1  using 
eqns.  (4.43)  and  (4.44). 

3.  Assuming  an  equilibrium  crack  opening  profile  vJ(x)  using  eqn.  (4.44)  with  a stress 
oJ  much  smaller  than  a*1,  calculate  the  corresponding  bridging  stress  distribution 
°bridgeKx)  using  the  ductile  phase  constitutive  relations  derived  in  Section  4.5. 

4.  Calculate  the  crack  opening  profile  vbjjdgeKx)  due  to  the  % ridge*  (x)  using  eqn. 
(4.46). 

5.  Check  eqn.  (4.42).  If  vJ(x)  = vext^(x)  - V5rj<jge*(x)>  then  the  assumed  vJ(x)  is  the 
equilibrium  crack  opening  profile.  If  not,  repeat  step  3,  4 and  5 but  with  a new 
a)  + 1 slightly  larger  than  cd  until  eqn.  (4.42)  is  satisfied  or 

v^(x) " { vext*(x) _ vbridge*(x)  } ^ 0.1%  vJ(x). 

6.  After  the  equilibrium  crack  opening  profile  is  found,  Kqp  is  calculated  using  eqns. 
(4.41),  (4.43)  and  (4.45).  The  aj  in  eqn.  (4.45)  is  determined  from  v(x)  using  the 
ductile  phase  constitutive  relations  derived  in  Section  4.5. 

7.  If  ( Km  - Ktjp  ) < 1%  Km,  then  the  crack-growth  resistance  Kc  is  calculated  from 
eqn.  (4.36)  or  (4.43).  If  not,  repeat  steps  1 to  7 with  a new  external  stress,  + 1, 
slightly  larger  than  cr*i  until  the  inequality  is  satisfied.  For  this  study,  Km  = 3.3 
MPa.  m 1/2  js  USed  which  is  the  fracture  toughness  of  MoSi2  matrix  measured  for  the 
present  MoSi2/Nb  system. 


4.7.3  Typical  R-Curve  Behavior  of  Ductile-Phase-Reinforced  Composites 

A calculated  curve  of  the  crack-growth  resistance  vs.  crack  extension  has  been 
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shown  in  Fig.  4.72.  For  the  present  computation,  laminated  composites  are  assumed 
and  the  crack  propagation  is  supposed  to  be  normal  to  the  ductile  laminae.  In  addition, 
the  initial  crack  length  is  assumed  to  be  the  distance  between  two  consecutive  ductile 
laminae.  Then  the  crack-growth  resistance  is  calculated  correspondingly  (see  Section 
4.7.2).  This  calculation  rule  is  also  followed  in  the  rest  of  this  chapter,  unless 
otherwise  mentioned.  As  expected,  the  crack-growth  resistance  increases  initially  as  the 
crack  extends  (from  the  original  crack  size  to  point  B),  because  more  and  more  ductile 
laminae  participate  in  bridging.  However,  as  the  crack  extends  more,  this  general  trend 
may  change  because  another  competitive  factor  becomes  pronounced,  that  is,  the 
central  ductile  ligaments  neck  down  or  even  fracture  due  to  the  large  crack  face 
displacement  at  the  center  of  the  crack.  At  a certain  crack  length  (point  B),  the 
computation  exhibits  a bifurcation:  one  shows  increasing  Kc  and  the  other  decreasing 
Kc  as  the  crack  extends.  Note  that  the  initial  curve  before  point  B has  Ktjp  = 3.3 
MPa.ml/^  while  after  point  B the  value  of  Ktjp  for  the  upper  and  lower  curves  are  not 
equal.  The  reason  for  this  discrepancy  will  be  provided  below  from  the  examination  of 
the  equilibrium  crack  opening  profile  and  bridging  stress  distribution  as  a function  of 
the  crack  length. 

The  equilibrium  crack  opening  profile  and  bridging  stresses  across  crack  face 
for  different  crack  lengths  are  shown  in  Fig.  4.73  where  (a),  (b)  and  (c)  correspond  to 
the  conditions  marked  as  A,  B and  C in  Fig.  4.72.  As  can  be  seen,  for  all  conditions, 
the  crack  opening  becomes  much  smaller  with  the  presence  of  bridging  ligaments  than 
the  unbridged  crack.  At  small  crack  length  (Fig.  4.73(a)),  the  maximum  bridging  stress 
is  almost  at  the  center  of  the  crack.  As  the  crack  length  increases,  the  maximum 
bridging  stress  moves  toward  the  vicinity  of  the  crack  tip  (Fig.  4.73(b)  and  (c)).  When 
the  crack  is  long,  it  is  found  that  the  profiles  of  the  equilibrium  crack  opening  and 
bridging  stresses  become  very  sensitive  to  the  external  stress  at  a certain  level  of 
loading  (to  be  discussed  more  in  Section  4.7.8),  as  shown  in  Fig.  4.73(c)  (note  KeXt 
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Fig.  4.73  The  equilibrium  crack  opening  profile  and  the  corresponding  bridging  tractions.  Crack  length,  2a,  is  2.1  mm  in  (a), 
12.1  mm  in  (b)  and  24.1  mm  in  (c).  Figures  (a),  (b)  and  (c)  correspond  to  A,  B and  C in  Fig.  4.72,  respectively. 
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a0).  Correspondingly,  the  effective  stress  intensity  Ktjp  changes  sharply  from  0.9  to 
12.0  MPa.ml/2  when  Kext  increases  only  by  0.3  MPa.ml/2  (about  1%  increase  in 
Kext).  The  sharp  increase  in  the  effective  stress  intensity  is  due  to  the  dramatic 
decrease  in  the  bridging  stresses  because  of  the  necking  of  the  central  ligaments,  as 
indicated  in  Fig.  4.73(c).  The  result  implies  that  = 24.6  MPa.ml/2  at  point  C 
may  not  be  dependable  because  a very  small  disturbance  in  the  external  stress  could 
make  Ktjp  jump  from  0.9  to  12  MPa.  m 1/2,  a value  much  higher  than  the  fracture 
toughness  of  the  matrix  (3.3  MPa.  m 1/2).  Furthermore,  if  the  opening  of  the  crack 
center  is  examined  as  a function  of  crack  length  as  presented  in  Fig.  4.74,  it  is  found 
that  the  opening  of  the  crack  center  increases  gradually  as  the  crack  extends  from  the 
original  crack  size  to  point  B.  However,  a bifurcation  occurs  at  point  B.  After  point  B, 
the  upper  curve  where  Ktjp  > > Km  (3.3  MPa. m 1/2)  exhibits  a fast  increase  in  the 
opening  of  the  crack  center,  indicating  that  the  ligament  necking  becomes  dominant 
after  a certain  crack  extension.  In  contrast,  the  lower  curve  where  Ktjp  < < 3.3 
MPa. m 1/2  has  a opening  smaller  than  that  at  point  B,  suggesting  an  impossible 
physical  state  because  it  implies  that  the  crack  closure  happens  when  the  crack  extends. 
Therefore,  the  upper  curve  in  Fig.  4.74  is  taken  as  the  real  behavior  of  the  composites. 
Correspondingly,  the  lower  curve  after  point  B in  Fig.  4.72  computed  from  the  larger 
crack  opening  (Ktjp  > > 3.3  MPa. m 1/2)  is  considered  as  the  representative  crack- 
growth  resistance  of  the  composite.  Thus,  point  B in  Fig.  4.72  is  regarded  as  the  point 
where  the  unstable  crack  propagation  begins  and  the  crack-growth  resistance  at  point  B 
is  considered  as  the  maximum  crack-growth  resistance,  Kc(max),  available  for  the 
composite.  The  other  two  important  parameters  describing  the  characteristics  of  the  R- 
curve  are  the  initial  slope  of  the  curve  and  the  crack  length,  ac(max),  corresponding  to 
the  Kc(max).  The  former  dominates  the  crack-growth  resistance  when  the  crack  is 
small,  while  the  latter  indicates  the  sensitivity  of  ductile  phase  toughening  to  the  crack 
length. 
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Fig.  4.74  The  equilibrium  displacement  of  crack  center  as  a function  of  crack  extension. 

The  A,  B and  C correpond  to  those  in  Fig.  4.72  and  the  calculation  parameters 
are  the  same  as  Fig.  4.72. 


Fig.  4.75  R curve  behaviors  of  composites  with  different  decohesion  lengths  at  the 
interface.  (20  vol%  of  Nb) 
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4.7.4  Prediction  of  the  Effect  of  Matrix/Reinforcement  Interface  on  Toughness 

The  crack-growth  resistance  curves  of  the  composites  with  different  interfacial 
decohesion  length  are  presented  in  Fig.  4.75.  It  can  be  seen  that  the  maximum  crack- 
growth  resistance  increases  with  increasing  decohesion  length,  while  the  initial  slope  of 
the  crack-growth  resistance  is  sharper  for  the  short  decohesion  composites  than  for  the 
long  decohesion  counterparts.  Because  of  the  sharper  slope  of  the  R-curve,  the  short- 
decohesion  composite  has  a larger  crack-growth  resistance  than  the  long-decohesion 
counterpart  when  crack  is  small.  Decohesion  length  at  the  interface  has  been  shown  to 
be  related  to  the  fracture  energy  of  the  interface,  i.e.,  the  higher  the  fracture  energy  of 
the  interface,  the  shorter  the  decohesion  length  (Section  4.5.2).  Furthermore,  it  has 
been  shown  in  Section  4.6.1  that  the  fracture  energy  of  the  chevron  notched  MoSi2/Nb 
composites  increased,  while  the  peak  load  to  fracture  the  same  specimens  decreased 
with  decreasing  the  fracture  energy  of  the  interface,  a paradox  in  choosing  the  bonding 
strength  of  the  interface.  The  present  theoretical  analysis  also  suggests  a trade-off  in  the 
selection  of  the  interfacial  properties,  that  is,  if  short  cracks  are  the  main  concern  of  the 
composite,  a relatively  strong  interface  should  be  tailored.  On  the  other  hand,  if  the 
composite  is  to  sustain  large  cracks,  a relatively  weak  interface  should  be  designed. 
Finally,  it  should  be  noted  that  a very  weak  interface  (decohesion  length  = 0.04  mm) 
may  not  be  desirable  because  the  large  toughening  is  reached  only  after  extensive  crack 
propagation. 

4.7.5  Prediction  of  the  Effect  of  Size  of  Ductile  Phase  on  Toughness 

If  different  sizes  of  Nb  laminae  are  assumed  to  have  the  same  intrinsic 
mechanical  properties,  then  the  corresponding  R-curves  of  the  composites  are  those 
presented  in  Fig.  4.76.  It  is  noted  that  although  the  composites  contain  different  sizes 
of  ductile  phases,  the  overall  R-curve  behaviors  are  very  similar,  i.e.,  the  crack-growth 
resistances  follow  a similar  locus  as  crack  extends.  However,  the  maximum  crack- 
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growth  resistance  increases  with  increasing  size  of  the  ductile  phase,  indicating  that 
large  size  of  ductile  phase  is  more  capable  of  toughening  when  crack  is  large.  This 
analysis  is  consistent  with  the  experiments  on  the  chevron  notched  composites  where 
the  initial  crack  length  is  about  1 mm,  as  shown  in  Section  4.6  and  in  Reference  132. 

For  engineering  applications,  the  crack  size  is  usually  smaller  than  few 
millimeters.  Thus,  if  we  concentrate  on  the  region  of  small  crack  size  (2a  < 6 mm) 
and  enlarge  that  portion,  we  have  Fig.  4.77.  Because  the  calculation  is  carried  out  on 
the  basis  of  constant  volume  fraction  of  ductile  phase,  the  initial  crack  size  (assumed  to 
be  equal  to  the  distance  between  two  consecutive  ductile  reinforcements)  increases  with 
increasing  size  of  ductile  phase.  Before  the  crack  grows  to  the  initial  crack  size,  the 
crack-growth  resistance  of  the  composite  is  regarded  to  be  equal  to  the  matrix 
toughness  unless  there  is  a residual  compressive  stress  developed  in  the  matrix  and/or  a 
strong  toughening  from  the  ductile  phase  shielding  is  present.  The  latter  is  usually  small 
[23].  For  the  case  of  no  residual  stresses,  therefore,  the  toughening  comes  earlier  for 
small  size  of  ductile  reinforcement  than  for  large  size  of  the  reinforcement,  as  shown  in 
Fig.  4.77.  This  is  better  appreciated  if  the  R-curve  is  converted  to  the  fracture  stress  of 
the  composites  using  the  equation  of 


The  result  is  presented  in  Fig.  4.78.  As  can  be  seen,  at  the  region  of  small  crack  size 
the  composite  containing  small  size  of  ductile  phases  has  a much  higher  fracture  stress 
than  those  containing  larger  ductile  phases.  For  example,  if  the  initial  crack  size  is 
limited  to  80  pm,  then  the  composite  containing  20  pm  size  of  ductile  phase  has  a 
fracture  stress  of  about  460  MPa,  while  the  counterpart  with  125  pm  size  ductile  phase 
has  only  a fracture  stress  of  about  280  MPa.  Therefore,  analogous  to  the  effect  of  the 
interface,  the  selection  of  the  size  of  ductile  phases  depends  on  the  major  crack  size  to 
be  dealt  with. 
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Fig.  4.76  R-curve  behaviors  of  composites  with  different  sizes  of  Nb  laminae 
( 20  vol%  of  Nb). 


Fig.  4.77  An  enlarged  view  of  the  R-curves  at  the  region  of  the  small  crack  size 
in  Fig.  4.76. 
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Fig.  4.78  Fracture  stresses  of  the  composites  with  different  sizes  of  Nb  laminae,  calculated 
from  the  data  in  Fig.  4.77. 


Fig.  4.79  R-curves  of  the  composites  containing  different  ducide  laminae.Volume  fraction 
of  the  ductile  laminae  are  20  % for  all  the  composites. 
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4.7.6  Prediction  of  the  Effect  of  the  Intrinsic  Properties  of  Ductile  phase  on  Toughness 

The  intrinsic  properties  of  the  ductile  phases  used  in  the  calculation  are  obtained 
from  Reference  185  and  have  been  listed  in  Table  4.13.  The  corresponding  R-curves 
are  presented  in  Fig.  4.79.  It  is  noted  that  the  best  toughening  is  achieved  by  0.6%  C 
Steel  which  has  the  highest  elastic  modulus  and  yield  strength  among  the  ductile  phases 
considered.  4340  steel  has  a similar  elastic  modulus  and  work  hardening  rate  to  0.6%  C 
Steel,  but  it  exhibits  a much  lower  toughening  than  0.6%  C Steel,  indicating  a strong 
influence  of  the  yield  strength.  Comparison  between  4340  steel  and  Nb  suggests  that 
higher  elastic  modulus  has  better  toughening,  because  4340  steel  and  Nb  have  similar 
yield  strengths  and  work  hardening  rates  but  different  elastic  moduli.  Effect  of  work 
hardening  rate  can  be  obtained  by  comparing  70/30  brass  with  Nb  both  of  which  have 
similar  elastic  moduli.  Although  70/30  brass  has  a higher  work  hardening  rate,  it  shows 
a lower  toughening  than  Nb  at  small  crack  size  but  a large  toughening  after  extensive 
crack  propagation,  indicating  the  yield  strength  has  a better  toughening  than  work 
hardening  rate  at  small  crack  size. 

Based  on  the  above  discussion,  it  can  be  summarized  that  the  maximum  crack- 
growth  resistance  increases  as  the  yield  strength,  elastic  modulus  and  work  hardening 
rate  of  ductile  phases  increase.  However,  the  yield  strength  and  elastic  modulus  are 
more  efficient  in  enhancing  toughening  than  work  hardening  rate  in  the  sense  that  the 
formers  can  provide  high  crack-growth  resistance  at  a relatively  short  crack  length 
while  the  latter  provides  large  toughening  only  after  extensive  crack  propagation. 
Similar  results  were  also  found  by  Ravichandran  [135]  and  Cox  [136]. 

4-7.7  Prediction  of  the  Effect  of  Volume  Fraction  of  Ductile  Phase  on  Toughness 

The  R-curves  of  the  composites  with  different  volume  fraction  of  ductile  phase 
are  presented  in  Fig.  4.80.  As  expected,  high  volume  fraction  of  ductile  reinforcement 
gives  rise  to  not  only  a high  maximum  crack-growth  resistance  but  also  a sharp  initial 
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Fig.  4.80  Effect  of  volume  fraction  of  ductile  phase  on  the  R-curve  behavior  of  the 
composites  containing  Nb  laminae  with  a thickness  of  0.02  mm  and  a 
decohesion  length  of  0.005  mm  at  the  interface. 


Table  4.18  Comparison  of  Composite  Toughness  Calculated  by  Different  Methods 
(SSB  and  LSB) 
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0.05  and  0.02  mm  Nb,  respectively. 

(b)  Thickness  of  Nb  = 0.02  mm;  Vol.  % of  Nb  = 20%. 

(c)  Thickness  of  Nb  = 0.02  mm;  Decohesion  length  = 0.005  mm. 


177 


slope  of  the  R-curve.  However,  the  maximum  crack-growth  resistance  occurs  at  a 
smaller  crack  length  when  the  volume  fraction  increases,  indicating  that  the  toughening 
becomes  more  sensitive  to  the  crack  length  as  the  volume  fraction  increases. 

At  this  point,  a comparison  can  be  made  between  the  calculations  based  on  the 
assumption  of  the  small-scale  bridging  and  the  large  scale  bridging.  The  conclusions 
made  in  Section  4.5  are  based  on  the  assumption  of  small-scale  bridging,  and  the 
fundamental  equations  in  SSB  have  been  shown  in  Section  2.3  (eqns.  (2.6)  and  (2.7)). 
The  results  calculated  using  eqns.  (2.6)  and  (2.7)  are  listed  in  Table  4.18  and  are 
compared  to  the  present  large-scale-bridging  calculations.  Since  the  crack-growth 
resistance  varies  with  the  crack  extension  in  the  present  large-scale-bridging 
calculations,  the  maximum  crack-growth  resistance  Kc(max)  is  taken  for  the 
comparison.  As  expected,  the  two  methods  lead  to  different  toughnesses,  a result  also 
found  by  Ravichandran  [135].  For  the  comparison  listed  in  Table  4.18,  the  large-scale- 
bridging  exhibits  a higher  toughness.  However,  if  the  crack-growth  resistance  after 
passing  the  Kc(max)  is  used  to  compare  with  the  steady-state  toughness,  it  is  found  that 
kc(SSB)  shows  a higher  toughness.  For  example,  a composite  containing  0.02  mm 
thick  Nb  laminae  exhibits  a crack-growth  resistance  of  about  9 MPa.  m 1/2  after  passing 
Kc(max)  (Fig-  4.76),  in  contrast  with  12.36  MPa. ml/2  calculated  from  the  small-scale- 
bridging method.  The  discrepancy  lies  obviously  in  the  different  assumptions  for  the 
two  approaches.  However,  since  large-scale  bridging  is  more  appropriate  for  the 
ductile-phase-toughened  composites,  it  is  expected  to  provide  a better  estimation. 

4.7.8  Prediction  of  the  Effect  of  Properties  of  Matrix  on  Toughness 

Effect  of  the  matrix  toughness  on  ductile  phase  toughening  can  be  evaluated 
from  Fig.  4.81  which  shows  how  various  stress  intensity  factors  change  with  the 
externally  applied  stress  for  a given  crack  length.  If  a certain  point  of  the  external  stress 
in  Fig.  4.81  is  taken  as  the  stress  to  propagate  the  crack,  then  the  corresponding 
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~ 89  113  137  162  186  210  234 


External  Stress,  MPa 


Fig.  4.81  Various  stress  intensity  factors  vs  external  stress  at  given  crack  lengths:  (a)  a = 
2.06  mm;  (b)  a = 5.06  mm;  and  (c)  a = 10.06  mm.  The  composite  contains  20 
vol.%  of  Nb  laminae  with  a thickness  of  0.02  mm  and  a decohesion  length  of 
0.005  mm  at  the  interface. 
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values  of  Kext,  Ktjp  and  K^^dge  can  be  regarded  as  the  crack-growth  resistance  of  the 
composite,  the  fracture  toughness  of  the  matrix  and  the  bridging  contribution  at  the 
time  of  the  propagation,  respectively.  The  stress  intensity  numbers  listed  in  the  figure 
indicate  the  maximum  bridging  contribution  available  for  the  given  conditions  and  the 
corresponding  matrix  toughness.  Thus,  from  the  Fig.  4.81  it  is  noted  that  the  bridging 
contribution  depends  on  the  fracture  toughness  of  the  matrix.  If  the  toughness  of  the 
matrix  is  low,  the  bridging  contribution  will  account  for  most  of  the  crack-growth 
resistance  of  the  composite.  However,  when  the  toughness  of  the  matrix  is  high  (Km  > 
10  MPa.ml/2),  the  bridging  contribution  becomes  not  only  smaller  in  magnitude  but  a 
small  portion  of  the  crack-growth  resistance  as  well.  This  result  is  believed  to  be 
related  to  the  following  scenario.  When  the  matrix  toughness  is  low,  the  crack  opening 
must  be  small  and  most  of  the  external  load  will  be  carried  by  the  bridging  ligaments 
before  the  crack  propagates,  leading  to  a large  contribution  from  ductile  phase.  In 
contrast,  when  the  matrix  toughness  is  high,  the  crack  opening  has  to  be  increased  to 
produce  enough  stress  intensity  at  the  crack  tip  before  the  crack  can  propagates.  As 
such,  the  bridging  ligaments  at  the  central  portion  of  the  crack  will  neck  down, 
resulting  in  less  contribution  from  the  bridging  ligaments. 

The  above  discussion  indicates  that  when  the  matrix  has  a high  toughness  (Km 
> 10  MPa.ml/2  ),  ductile  phase  toughening  is  not  efficient  unless  a ductile  phase  with 
a high  work  hardening  rate  and  a weak  interface  has  been  utilized.  Both  of  these 
material  properties  prevent  ductile  ligaments  from  necking.  In  reality,  however,  almost 
all  the  brittle  materials  have  toughness  below  10  MPa.m^  therefore  a more 
important  aspect  is  what  material  properties  affect  the  transition  of  the  bridging 
contribution  from  high  percentage  to  low  percentage  of  the  crack-growth  resistance, 
and  how  to  control  these  properties  for  optimal  ductile  phase  toughening.  Form  Fig. 
4.81,  it  is  obvious  that  the  transition  is  affected  by  the  crack  length,  i.e.,  the  longer  the 
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crack  length,  the  lower  the  matrix  toughness  at  which  the  transition  occurs.  Besides, 
the  transition  also  becomes  sharper  when  the  crack  length  increases. 

The  transition  is  also  affected  by  the  decohesion  length  at  the  interface. 
Comparison  between  Fig.  4.81  and  4.82  manifests  that  if  the  decohesion  length  is 
large,  the  transition  becomes  much  gentler  or  might  even  shows  no  decrease  in  the 
magnitude  of  the  bridging  contribution  (Fig.  4.82(a)).  It  is  also  noted  that  the  transition 
occurs  at  a higher  matrix  toughness  when  the  decohesion  is  long.  For  example,  the 
matrix  toughness  at  which  the  bridging  contribution  begins  to  decrease  changes  from 
2.2  to  5.1  MPa.ml/2  when  the  decohesion  length  changes  from  0.005  to  0.04  mm  for 
a crack  with  a = 5.06  mm  (Fig.  4.81(b)  and  4.82(b)).  The  results  suggest  that  when 
the  matrix  toughness  is  high,  a relatively  weak  interface  (thus  a long  decohesion  length) 
would  be  desirable  to  enhance  bridging  contribution,  otherwise  the  bridging 
contribution  would  be  too  sensitive  to  the  size  of  the  existing  crack.  On  the  other  hand, 
if  the  matrix  toughness  is  low,  a relatively  strong  interface  may  be  beneficial  to  limiting 
the  bridging  length  and  producing  a higher  crack-growth  resistance  at  a shorter  crack 
length. 

Effect  of  the  volume  fraction  of  the  ductile  phase  on  the  transition  is  presented 
in  Fig.  4.83.  Comparison  between  Fig.  4.81  and  4.83  reveals  that  the  transition 
becomes  sharper  with  increasing  volume  fraction  of  ductile  phase,  or  in  other  words,  it 
becomes  more  sensitive  to  the  crack  length.  However,  a much  higher  bridging 
contribution  at  the  onset  of  the  transition  is  reached  by  the  high  volume  fraction  of  the 
ductile  phase.  For  instance,  23.5  and  14.3  MPa.m^^  0f  the  bridging  contributions  is 
obtained  at  the  onset  of  the  transition  for  40  and  20  vol.%  of  the  ductile  phase, 
respectively  (Fig.  4.81(b)  and  4.83(b)). 

The  effects  of  various  material  properties  on  the  transition  can  also  be  deduced 
from  the  R-curves  (Fig.  4.75,  4.76,  4.79  and  4.80).  Increasing  the  matrix  toughness 
will  shift  the  R-curve  upward  and  increase  the  initial  slope  of  the  curve,  but  move 
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Fig.  4.82  Various  stress  intensity  factors  vs  external  stress  at  given  crack  lengths:  (a)  a = 
1 .0  mm;  (b)  a = 5.06  mm;  and  (c)  a = 11 .06  mm.  The  composite  contains  20 
vol.%  of  Nb  laminae  with  a thickness  of  0.02  mm  and  a decohesion  length  of 
0.04  mm  at  the  interface. 
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External  Stress,  MPa 


External  Stress,  MPa 


Fig.  4.83  Various  stress  intensity  factors  vs  external  stress  at  given  crack  lengths:  (a)  a = 
3.035  mm  and  (b)  a = 4.035  mm;  and  (c)  a = 11 .06  mm.  The  composite 
contains  40  vol.%  of  Nb  laminae  with  a thickness  of  0.02  mm  and  a decohesion 
length  of  0.005  mm  at  the  interface. 
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^(max)  t0  a smaller  crack  size.  Based  on  such  reasoning,  the  selection  of 
material  properties  can  be  made  by  referring  to  the  R-curves.  For  example,  if  a matrix 
with  a higher  toughness  than  that  shown  in  Fig.  4.75  is  to  be  ductile-phase- toughened, 
the  decohesion  length  of  0.001  mm  is  obviously  not  a good  choice  because  the 
composite  would  be  even  more  sensitive  to  the  size  of  the  existing  crack.  From  Fig. 
4.76,  the  transition  is  expected  to  be  insensitive  to  size  of  the  ductile  phase,  because 
different  sizes  of  ductile  phases  follow  a similar  locus  of  the  crack-growth  resistance. 
Fig.  4.79  suggests  that  a high  work  hardening  rate  material  such  as  70/30  brass  would 
be  less  sensitive  to  the  matrix  toughness  and  the  crack  size.  In  a word,  the  effects  on 
the  ductile  phase  toughening  of  the  matrix  toughness,  the  bonding  strength  of  the 
interface,  and  the  intrinsic  properties,  size  and  volume  fraction  of  the  ductile  phase  are 
interdependent  with  each  other.  The  selection  of  these  material  properties  should  be 
based  on  the  crack  size  being  concerned  and  the  matrix  toughness. 

In  summary,  from  the  analysis  of  the  large-scale  bridging,  the  following 
conclusions  can  be  made: 

1.  Ductile-phase-toughened  composites  will  exhibit  a rising  crack-growth 
resistance  (R-curve)  which  can  be  characterized  by  three  parameters:  the  maximum 
crack-growth  resistance  Kc(max),  the  corresponding  crack  length  2ac(max),  and  the 
initial  slope  of  the  R-curve.  Crack-growth  resistance  is  mainly  determined  by  the  initial 
slope  of  the  R-curve  when  crack  is  small  and  is  dominated  by  the  maximum  crack- 
growth  resistance  when  crack  is  large. 

2.  The  maximum  crack-growth  resistance  increases  with  decreasing  the  fracture 

energy  of  the  matrix/reinforcement  interface.  However,  2ac(max)  increases  as  well, 
suggesting  that  is  reached  only  after  extensive  crack  propagation.  On  the  other 

hand,  the  interface  with  a high  fracture  energy  has  a high  initial  slope  of  the  R-curve. 
The  duality  of  the  interface  suggests  that  a weak  interface  would  be  desirable  if  the 
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cracks  concerned  are  long,  and  a strong  interface  would  be  beneficial  if  the  cracks  to  be 
dealt  with  are  short. 

3.  The  overall  R-curve  behaviors  of  composites  with  different  sizes  of  ductile 
phases  are  very  similar;  they  follow  a similar  locus  as  crack  extends.  However, 
Kc(max)  increases  with  increasing  size  of  ductile  phase,  indicating  that  large  size  of 
ductile  phases  are  more  capable  of  toughening  when  cracks  are  large.  When  the  cracks 
are  small  (2a  < 200  pm)  which  is  of  great  interest  to  the  engineering  applications, 
composites  containing  small  size  of  ductile  phases  show  a higher  fracture  stress  than  the 
counterparts  with  large  size  of  ductile  phases.  This  paradox  is  very  similar  to  the  choice 
of  the  interface,  i.e.,  selecting  these  properties  depends  on  the  crack  size  concerned. 

4.  Kc(max)  increases  with  increasing  yield  strength,  elastic  modulus  and  work 
hardening  rate  of  ductile  phases.  However,  the  yield  strength  and  elastic  modulus  are 
more  efficient  in  enhancing  toughening  than  work  hardening  rate  in  the  sense  that  the 
formers  can  provide  high  crack-growth  resistance  at  a relatively  short  crack  length, 
while  the  latter  provides  large  toughening  only  after  extensive  crack  propagation. 

5.  High  volume  fraction  of  ductile  reinforcement  gives  rise  to  not  only  a high 
Kc(max)  but  a sharp  initial  slope  of  the  R-curve.  But  2ac(max)  also  becomes 
smaller  when  the  volume  fraction  increases,  indicating  that  the  toughening  becomes 
more  sensitive  to  the  crack  length  as  the  volume  fraction  increases. 

6.  When  the  matrix  toughness  is  low,  a ductile  phase  with  high  yield  strength 
and  elastic  modulus  and  a relatively  strong  interface  would  be  desirable  to  produce  a 
high  crack-growth  resistance  at  a short  crack  length.  If  the  matrix  toughness  is  high, 
other  material  properties,  such  as  weak  interface  and  high-work-hardening-rate 
materials  which  prevent  bridging  ligaments  from  necking,  should  be  chosen  to  enhance 
bridging  contribution. 


CHAPTER  FIVE 
CONCLUSIONS 


Ductile  phase  toughening  in  brittle  matrices  for  high  temperature  structural 
applications  has  been  investigated.  Owing  to  the  projected  use  temperatures,  the  system 
selected  for  the  study  was  MoSi2  matrix  and  Nb  ductile  reinforcements.  The  main 
conclusions  from  this  study  are  as  follows. 

1.  For  controlling  the  matrix/reinforcement  interaction,  chemical  compatibilities 
of  various  potential  coating  materials  (mullite,  AI2O3,  and  Zr02)  with  MoSi2  and  Nb 
were  studied  via  microstructural  and  compositional  analyses.  Based  on  these  analyses, 
AI2O3  was  chosen  as  the  final  diffusion  barrier  coating  because  it  is  chemically 
compatible  with  both  MoSi2  and  Nb,  and  it  also  has  a thermal  expansion  coefficient 
close  to  those  of  both  the  matrix  and  reinforcement. 

2.  The  techniques  investigated  for  producing  AI2O3  coating  on  Nb  surface 
included  (1)  sol-gel  processing;  (2)  physical  vapor  deposition;  (3)  hot  dipping  the 
reinforcement  in  molten  aluminum,  followed  by  anodizing  A1  to  form  AI2O3.  The 
results  indicated  that  a thick  and  dense  AI2O3  coating  is  necessary  for  effective 
elimination  of  the  interactions  between  Nb  and  the  matrix  of  commercially  pure  MoSi2. 
The  results  also  showed  that  the  best  coating  is  provided  with  physical  vapor 
deposition,  followed  by  the  sol-gel  technique. 

3.  Effects  of  the  matrix/reinforcement  interface,  the  mechanical  properties  and 
size  of  the  ductile  phase  on  the  flow  behavior  of  constrained  ductile  phases  were 
evaluated  via  tensile  tests  on  a single  Nb  lamina  imbedded  in  MoSi2  matrix.  An 
approximate  model  was  proposed  which,  according  to  the  small-scale-bridging 
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assumption,  gave  the  following  conclusions:  the  toughness  of  the  composites  is 
enhanced  by  a relatively  weak  bond  at  the  matrix/reinforcement  interface,  by  a large 
size  of  ductile  reinforcements,  and  by  a ductile  phase  with  a high  yield  strength  and 
high  work  hardening  rate. 

4.  The  role  of  the  matrix/reinforcement  interface  in  the  fracture  toughness  of 
MoSi2/Nb  composites  was  assessed  via  four-point  bend  tests  on  chevron  notched 
specimens.  It  has  been  established  that  whether  or  not  a strong  interfacial  bonding  is 
conducive  to  toughness  depended  on  the  criterion  used  to  describe  the  toughness  of  the 
composites.  If  the  peak  load  of  the  chevron  notched  specimen  is  used  as  an  indicator  of 
the  toughness,  a strong  bonding  would  be  desirable.  On  the  other  hand,  if  the  total 
energy  consumed  to  break  a specimen  is  used  as  indicative  of  its  toughness,  then  a 
relatively  weak  bonding  is  preferable. 

5.  Four-point  bend  tests  were  used  to  evaluate  the  effects  of  the  size  and 
orientation  of  ductile  laminae  on  the  toughness  of  the  composites.  The  results  indicated 
that  toughness  of  the  composites  increased  with  increasing  size  of  Nb  laminae,  and  that 
ductile  laminae  offered  two  dimensional  toughening.  The  fracture  toughness  in  the  two 
orientations  (N-  and  T-orientations)  has  been  increased  to  about  11-14  MPa.  m 1/2 
from  3 MPa. ml/2  for  unreinforced  matrix. 

6.  All  the  toughness  evaluations  and  the  examination  of  the  interaction  of  cracks 
with  ductile  reinforcements  in  the  present  study  indicated  that  (1)  ductile  phase 
toughening  involved  large-scale  bridging,  i.e.,  the  bridging-length  was  at  the  same 
magnitude  as  the  crack  length,  the  specimen  size  or  the  distance  from  the  crack  to  the 
specimen  boundaries,  and  (2)  the  toughening  of  brittle  matrix  with  ductile  phase  could 
be  described  by  the  equilibrium  stress  distribution  across  the  crack  face. 

7.  According  to  the  aforementioned  observations,  a general  approach  has  been 
proposed  to  compute  the  ductile  phase  toughening  in  the  case  of  large-scale-bridging. 
The  approach,  which  allowed  for  the  prediction  of  the  typical  R-curve  behaviors  of 
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ductile-phase-toughened  composites,  encompassed  the  effects  of  the  bonding  strength  of 
the  matrix/ reinforcement  interface,  the  size  of  ductile  phases,  and  the  intrinsic 
mechanical  properties  of  the  ductile  phases  and  matrix.  From  the  analysis  of  the  large- 
scale  bridging,  the  following  conclusions  can  be  made:  (1)  Ductile-phase-toughened 
composites  will  exhibit  a rising  crack-growth  resistance  (R-curve)  which  can  be 
characterized  by  three  parameters,  i.e.,  the  maximum  crack-growth  resistance 
Kc(max)>  the  corresponding  crack  length  2ac(max),  and  the  initial  slope  of  the  R-curve. 
Crack-growth  resistance  is  mainly  determined  by  the  initial  slope  of  the  R-curve  when 
crack  is  small  and  by  the  maximum  crack-growth  resistance  when  crack  is  large.  (2) 
The  interface  with  a short  decohesion  length  gives  rise  to  a steeper  initial  slope  of  R- 
curve,  while  Kc(max)  increases  with  the  decohesion  length  at  the  interface.  The  duality 
of  the  interface  suggests  that  a weak  interface  would  be  desirable  if  the  crack  concerned 
are  long,  but  a strong  interface  would  be  beneficial  if  the  cracks  to  be  dealt  with  are 
short.  (3)  increases  with  increasing  size  of  the  ductile  phase,  indicating  that 

large  size  of  ductile  phases  are  more  capable  of  toughening  when  cracks  are  large. 
When  the  cracks  concerned  are  small  (2a  < 200  pm)  which  is  of  great  interest  to  the 
engineering  applications,  composites  containing  small  size  of  ductile  phases  show  a 
higher  fracture  stress  than  the  counterparts  with  large  size  of  ductile  phases.  This 
paradox  is  very  similar  to  the  choice  of  the  interface,  i.e.,  selecting  these  properties 
depends  on  the  crack  size  concerned.  (4)  Kc(max)  increases  with  increasing  yield 
strength,  elastic  modulus  and  work  hardening  rate  of  ductile  phases.  (5)  High  volume 
fraction  of  ductile  reinforcement  gives  rise  to  not  only  a high  but  also  to  a 

sharp  initial  slope  of  the  R-curve.  (6)  When  the  matrix  toughness  is  low,  a ductile 
phase  with  high  yield  strength  and  elastic  modulus  and  a relatively  strong  interface 
would  be  desirable  to  produce  a high  crack-growth  resistance  at  a not  very  long  crack 
length.  If  the  matrix  toughness  is  high,  other  material  properties,  such  as  weak 
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interface  and  high-work-hardening-rate  materials  which  prevent  bridging  ligaments 
from  necking,  should  be  chosen  to  enhance  bridging  contribution. 


APPENDIX 

FINITE  ELEMENT  ANALYSIS  OF  STRESS  DISTRIBUTION  ACROSS  THE 
NECK  OF  THE  CONSTRAINED  DUCTILE  PHASES 


Finite  element  analysis  (FEA)  of  the  stress  distribution  across  the  neck  and  the 
contour  of  the  decohesion  region  during  the  elastic  deformation  was  implemented 
using  the  finite  element  software  package  ANSYS.  Due  to  the  symmetry,  only  one 
quarter  of  the  specimen  was  analyzed  using  two-dimensional  4-node  isoparametric 
elements  with  an  assumption  of  plane  strain  deformation.  A typical  finite  element 
model  for  the  constrained  Nb  is  shown  in  Fig.  A.  The  displacement  at  the  midplane 
of  the  neck  was  allowed  only  in  the  x-direction,  and  the  bonded  boundary  of  Nb  to 
the  matrix  and  the  axial  central  plane  were  allowed  to  displace  only  in  z-direction,  as 
shown  in  Fig.  A. 

Results  from  the  FEA  for  the  case  of  constrained  1.0  mm  Nb  lamina  with  2.28 
mm  decohesion  length  under  a nominal  axial  stress  of  150  MPa  are  presented  in 
Figures  B-D.  Figure  B compares  the  assumed  contour  of  the  decohesion  region  with  the 
FEA  result.  The  data  show  that  the  assumed  contour  deviates  from  the  FEA  result, 
leading  to  a smaller  crack  opening  displacement  compared  to  the  finite  element 
analysis.  However,  the  difference  is  so  small  that  it  introduces  negligible  errors  to  the 
computation  of  the  work  of  rupture.  Evaluation  of  effective  stress  distribution  across 
the  midplane  of  the  neck  is  shown  in  Fig.  C.  Clearly,  the  effective  stress  is  not 
constant  across  the  midplane,  but  increases  from  the  free  surface  to  the  center  of  the 
neck.  Although  the  assumed  constant  effective  stress  does  not  represent  the  real  stress 
distribution,  the  mean  axial  stresses  across  the  midplane  are  almost  the  same  for  the 
model  and  FEA,  as  shown  in  Fig.  D.  Thus,  it  is  concluded  that  the  assumed  constant 
effective  stress  and  contour  of  the  decohesion  region  impose  negligible  errors  on  the 
calculation  of  a-u  curve  and  the  work  of  rupture. 


189 


190 


Fig.  A Mesh  and  boundary  conditions  used  for  the  analysis  of  constrained  Nb. 


Z,  mm 


Fig.  B A comparison  between  the  assumed  contour  and  finite  element  analysis  of 
the  decohesion  region  during  the  stage  of  elastic  deformation.  X and  Z are 
defined  in  Fig.  A. 
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Fig.  C Effective  stress  distribution  at  the  midplane  of  neck. 


Fig.  D A comparison  of  the  axial  stresses  at  the  midplane  of  neck  between  the  model 
and  FEA.  Note  that  mean  axial  stresses,  omean,  316  almost  the  same. 
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